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Capitolo 1 

Interazioni Nucleari 



1.1 Interazioni di scambio 

Abbiamo gia visto che il potenziale di interazione tra due nucleoni si puo scrivere come 

Vnn = Vi(r) + V 2 {r)a 1 ■ a 2 + V 3 (r)Si 2 (1.1) 

dove 

3(ai -r)(<7 2 -r) 

^12 = ' y 1 - 2 ) 

Tale potenziale deve essere corretto con termini che tengano conto dei possibili scambi di 
particelle nei processi di scattering. Si vogliono quindi descrivere tutti i possibili scambi in 



modo da quantificare il termine correttivo della (1.1). Vediamo innanzitutto quali sono le 
possibili operazioni di scambio. Supponiamo quindi di prendere la funzione ip(fi, a%, r*2, a 2 ) 
che rappresenta la funzione d'onda del sistema di due particelle. Definiamo l'operatore P r 
che agisce come P r ifj{fi, a\, r 2 , cr 2 ) = tp(r 2 , 0*1, ?i, 0i) ■ tale operatore e detto di Majorana. 
Analogamente si introduce un operatore di scambio P a che agisce come P a ip = ip(ri,cr 2 ,f 2 ,ai) 
: tale operatore e detto di Bartlett. Infine si introduce un'operatore , detto di Heisenberg , 
P ra - che agisce come P ra ip = ip(^2, cr 2 , fi, o\). Ovviamente vale 

P — P P ■ P 2 — P 2 — P 2 — IT 

quindi gli autovalori di questi operatori sono dbl. L'autovalore dell'operatore P r si puo scrivere 
in termini del moto relativo , infatti P r = (— 1) . La stessa cosa si puo fare per l'operatore 
P a , visto che nei casi di singoletto (5 = 0), tripletto (S = 1) la funzione d'onda risulta 
rispettivamente antisimmetrica e simmetrica. Dunque P a = (— l) 5+1 .Combinando queste due 
relazioni si ricavano anche gli autovalori dell'operatore di Heisenberg : P ra = (— l) i+s+1 . 
Analogamente si introducono 3 tipi di forze nucleari di scambio. 

• VAi( r )Pr — > forze di Majorana 

• V B {R)Pa — > forze di Bartlett 

• VH(r)P ra — > forze di Heisenberg 

A partire da queste forze, di tipo centrale ,si possono costruire anche delle forze di tipo 
tensoriale. Queste si ottengono dalle prime , moltiplicando per l'operatore tensore definito 



dalla (1.2). Si vede pero che esistono forze tensoriali solo di tipo Majorana, visto che le altre 
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si possono ricondurre a questa. Infatti proviamo ad applicare l'operatore tensoriale di spin 



alle forze di Bartlett. Visto che l'operatore della (1.2) si puo riscrivere come 



S12 



2 

¥ 



{S-f 



^2 



s 2 



allora 



Analogamente 



5*12 Pa 



S 12 Pr Pa 



se S = 
5i2 se S = 1 

'0 se S = 
S^Pr se S = 



1 



(1.3) 



(1.4) 



(1.5) 



Quindi un eventuale forza tensoriale di scambio di tipo Heisenberg, cosf come una di tipo 
Bartlett, si riduce alia forza di tipo Majorana. Dunque l'unico tipo di forza tensoriale puo 
essere solo quella di tipo Majorana ed e data da Vryii^SuPr- In particolare, se si considera 
un sistema di due nucleoni si puo anche considerare anche l'operatore di isospin per il quale 
viene definito, analogamente a quanto fatto prima, un operatore di scambio di isospin per 
cui P T ip(fi, Ti, r*2, o"2, T2) = ip(ri, cr±, T2, r*2, o"2, t±). Ovviamente anche per questo operatore 
valgono le relazioni 



Pi 



I 



p T = ±i 



Dunque l'operatore composto da tutti gli scambi equivale ad P\i = P r P a P T = —I : questa e 
la forma del principio di esclusione di Pauli. Abbiamo quindi costruito 6 tipi di forze 

Tipi centrali tensoriali 

Forze Ordinarie Vw{r) Vrw( r ')<S , i2 

Forze di scambio V M (r)P r , VB(r)Pg, Vg(r)P r(T VrM^S^-Pr 

Tutte queste forze compongono il potenziale Vg(r) che costituisce la correzione al potenziale 
nucleare di inter azione nucleone-nucleone. Vogliamo ora scrivere questo potenziale esplicitando 
la dipendenza dagli operatori a. Si utilizzano relazioni del tipcQ 



Prr 



1 



(1 + <7l • <7 2 ) 



Poiche gli operatori hanno forme simili si puo utilizzare anche una relazione del tipo 

1 



Pt 



-(l + f 1 -f 2 )^P r = (-l) 



r+l 



(1.7) 



(1.8) 



Utilizzando queste relazioni l'equazione di consistenza che definisce il principio di esclusione 
di Pauli si puo riscrivere come 



P r P a P T = -1 (-lY+S+l+r+l = 

La derivazione di queste relazioni e semplice : basta osservare che 

'+1 se S = 1 



S+7 



-1 



((Tl ■ CT 2 ) = 



Pa = (-1)' 



-3 se S = 

Piii formalmente si dovrebbe applicare l'operatore P„ agli stati di singoletto e tripletto. 



(1.9) 



(1.6) 
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ovvero t + S + r deve essere una combinazione dispari. Visto che l'operatore di Heisenberg 
equivale all'operatore di scambio di isospin, a meno di un segno , possiamo scrivere 

Pro = + n • f 2 ) (1.10) 

Utilizzando inoltre la relazione P r = —P a P T si ricava infine 

P r = + ^ ■ a 2 ){l + n ■ T 2 ) (1.11) 

Sostituendo questi risultati nel potenziale si ottiene finalmente 

V 6 =V c (r) + V a {r)ai ■ a 2 + V T {r)r x ■ f 2 + V T(T {r){a x ■ a 2 )(n ■ r 2 ) 
+ V T (r)S 12 + V TT {r)S 12 {fi ■ f 2 ) 

1.2 Teoria di Yukawa 



(1.12) 



Secondo la teoria di Yukawa l'interazione tra due cariche mediata da particelle, dette pioni. 
Graficamente , per l'interazione tra due nucleoni iVi ed N 2 , si utilizza la seguente notazione 




x 




Alcuni esempi di interazioni mediate da pioni sono riportate in figura. 

P n n V 

+ 







n p P n 



II primo quesito che Yukawa si pose fu quello di determinare la massa di questa particelle 
mediatrice : in questo caso la massa e evidentemente diversa da 0, a differenza di quanto 
succedeva con il fotone nell'interazione elettromagnetica. Questa particella non viene rivelata 
direttamente nel processo quindi deve essere ricavata con considerazioni generali. Quando il 
nucleone 1 emette la particella di mediazione x si sta violando la conservazione dell'energia 
, dunque questa particella puo esistere solo per un tempo di coerenza che si puo stimare 
utilizzando l'indeterminazione AEAt ~ H. Infatti , poiche AE = m x c 2 => At ~ h/m x c 2 . La 
distanza che questa particella puo percorrere si puo stimare cone 

. He H 

cAt ~ k = = Hint 

m x c z m x c 

Questo valore rappresenta l'ordine di grandezza del raggio di interazione ed e pari alia lunghez- 
za d'onda Compton. Supponendo i?i n t ~ 1,5 fm (~ raggio dell'interazione nucleare) si puo 
stimare 

2 he 200MeVfm 

m x c 2 = — = — — ~ 130MeV 

Hint 1, 5tm 
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1.3 Campo Pionico 

Vediamo come e possibile ricavare l'equazione che definisce implicitamente il campo pionico. 
Conviene prima stabilire un'analogia con il caso elettromagnetico. II campo prodotto da una 
carica puntiforme si puo scrivere risolvendo l'equazione di Poisson per una distribuzione di 
carica p(r) = e5 3 {r). L'equazione da risolvere e 



V 2 



- e 8 3 (f) <j>{r) 



1 



47T€o r 



Ovvero 



V(r) = &j> 



1 



(1.13) 



(1.14) 



4-7reo r 

Per trovare l'equazione analoga per il campo pionico consideriamo l'energia relativistica E 2 = 
p 2 c 2 + m 2 c 4 e sostituiamo E e p con i relativi operatori quantisticQ Si ottiene l'equazione di 
Klein-Gordon 



m 2 c 2 



c 2 dt 2 ' 



[1.16) 



dove 4> stavolta identifica il campo pionico. Si noti che ponendo m = si trova l'equazione di 
D'Alembert per le onde. Se invece consideriamo il caso stazionario dt4> = e l'equazione si 
riduce a V 2 </> — k 2 (p = , dove k = m x c/h : tutto questo vale nel caso di assenza di sorgenti. 
La sorgente del campo pionico e costituita dal nucleone : se supponiamo di posizionarlo nel 
punto r*o = l'equazione diventa 



k 2 4> 



-gd 3 (r). 



(1.17) 



La costante g e l'analogo della carica nel caso pionico. La soluzione dell'equazione e data dal 
potenziale di Yukawa 



9 e 



- kr 



4-7T r 



Vnn = g4> 



Un'ulteriore generalizzazione della (1.18) e data da 



fair) 



g e 



47r If — r\\ 



(1.18) 



;i.i9) 



dove A c = H/m^c e la lunghezza d'onda Compton e la sorgente del campo e data da g5 3 (f — r{). 
Si vede che il campo scritto prima e invariante per parita nonostante il pione abbia parita 
intrinseca negativa. Bisogna quindi modificare il termine di sorgente con g(B\ • Vi)<5 3 (f — fx). 



Sostituendo nell'equazione (1.17) si ottiene 



V 2_ 

h 2 



4>irif) 



-5(0*1 • Vi)5 3 (f-fi) 



;i-20) 



Questa equazione ammette la soluzione 

07r(f) = 



[g-fi 



f(5,-V,)f, 
47T r 



:i.2i) 
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Visto che {a\ ■ Vi) e uno pseudoscalare^J allora 0^ cambia segno sotto parita , come volevasi 
ottenere. II potenziale nucleare si puo quindi scrivere come 



„ P2-H.I 
g -> - e A c 

y NN = Vb PE (r) = -;-(<72 • V 2 )((fi • Vi)-^ 

47T |r 2 — r\\ 



(1.22) 



Da questa si ottiene immediatamenta^] 



4vr3 



( 1 + 3— + ] .V,, 



-r/A c 



;i.23) 



ove e l'operatore tensoriale definito dalla (1.2) e la sigla OPE sta per One Pion Exchange. 



Cambia segno per inversione dei 3 assi 



4 Passando al centro di massa si definisce una nuova variabile r = r\ — r*2 per la quale si definisce V r = V. 
Si ottiene quindi Vi = V e V2 = — V, ovvero (02 ■ V2)(o! ■ Vi) = — (<r 2 ■ V)(cfl ■ V) 
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Capitolo 2 

Fisica dei Pioni 



2.1 Proprieta generali 

I pioni possono essere di diversi tipi : ir + , tt° , ir~ . La massa di queste particelle e data 
rispettivamente da m w ± = 139, 57 MeV e m w o = 134, 97 MeV . 

Decadimenti Questi non sono stabili ma decadono usualmente secondo i seguenti decadi- 
menti 

• tt° — ► 7 + 7 , (98, 8%) 

• tt° ^7 + e+ + e-, (~ 1,2%) 

II tempo medio di decadimento e dato da r = (8, 4± 0, 6) • 10~ 17 s. Gli altri tipi di pioni 
decadono secondo 

• vr" — > pT + — > e~ + + v e , (~ 100%) , r ~ 2, 8 • 10~ 8 s 

• 7r~ — > e~ + v e , molto raro 

Massa La massa del pione si puo ricavare studiando questo decadimento del centro di massa. 
Vale infatti il bilancio 

m^c 2 = + p u c (2.1) 

Spin Si voglia misurare lo spin del pione. Dal primo decadimento illustrate , 7r° — > 7 + 7 , 
si deduce che lo spin deve essere intero. Si puo inoltre utilizzare il processo di scattering 
seguente 

p + p — > p + n + TT + = d + ir + 

ove d indica il deutone. Studiando anche il processo inverso e confrontando le due sezioni 
d'urto possiamo avere informazione sullo spin del pione. Si noti che la conseguenza 
dell'invarianza temporale per processi generici a -\-A — > b + B e b + B — > a + A impone 
che le sezioni d'urto siano le stesse a meno di cambiamenti di spin. Tale proprieta e 
detta del bil anci o dettagliato e nel seguito si dara una breve spiegazione del risultato 
ottenuto in (2.5) ,alla quale si rimanda immediatamente nel caso si voglia seguire una 
trattazione meno rigorosa. 

Si consideri un processo generico a + A — > b + B : il numero di transizioni per unita di tempo e 
dato dalla regola d'oro di Fermi 

2-7T 

w = —\(^ Bh \H\^ All )\ 2 p(E) (2.2) 
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dove p(E) = dn }^ e la densita di stati finali del sistema i unita di energia. Per una particella 
libera in un volume V si ha 



dn = (2s B + l)(2s 6 + l) 



AttVpI dp b 
h 3 



Visto che dE = Vb dpb si ha 



dn 
dE 



(2s B + l)(2s b + l)^^ 



Sostituendo nella (2.2) si ottiene 

V V\ 
w = — — — 

La sezione d'urto e definita come il rapporto 
n.transizioni 



(2s B + l)(2s b + l)\(iJjBh\H\ip Aa )\ 2 



N a v a a(a + A^b + B) 



n.nucleiA ■ tempo ' n.nucleiA 

Supponendo che ci sia un solo nucleo di A nel volume V si ottiene 

Vw _ V 2 p\ 



a(a + A -> b + B) 



-(2 SB + l)(2 S , + l)|(^Bb|ff|^aA)| Z 



(2.3) 



v a irh 4 v a v b 

Nel calcolare il rapporto tra le sezioni d'urto delle reazioni diretta e inversa conviene ricordarc 



che H deve essere Hermitiano, quindi K^Bbl-HI^Aa)! 2 = KV'Aal-HlV'Bb)! 2 - Si ricava quindi 

a(a + A -j. b + B) _ k 2 (2s B + l)(2s 6 + 1) 
a(b + B^a + A) ~ k 2 (2s A + l)(2s a + 1) 

Nel caso particolare del decadimento pp si ha 

a(p + p -> vr+ + d) kl{2s n + l){2s d + l) k% (2 S?r + 1)3 



(2.4) 



(2-5) 



a(ir+ + d^p + p) kj (2s p + l) 2 /2 kj 2 

II 2 al denominatore tiene conto del doppio conteggio effettuato , inoltre l'ultimo risultato 
e stato ottenuto sostituendo Sd = 1 e s p = 1/2. Sperimentalmente si osserva quindi che 
l'identita e vera quando s n = 0. 



PION LABORATORY ENERGY . MeV, { iwd — p * p ) 
SO 100 500 1000 2000 3000 



g 



• from p-p — - u**d 
from Tl'*c — . p * p 




300 500 TOOC ?00C 500C 1000C 

PROTON ^ABORATOR* ENERGY MeV (p.p— r**fl) 



2.2. RISONANZE 
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Parita La parita intrinseca del pione si ricava invece considerando il decadimento 

tt~ + d — > n + n 

Visto che l'interazione e forte e a corto range si puo considerare solo il contribute con 
£ = : per ricavare la parita bast a quindi applicare la conservazione IL = II/. In 
particolare visto che 

Ui = IL n U d (-l) ti = II* 

e che 

U f = U n U n (-l) e f = (-l) £ f 

si ottiene che IT,,- = (— 1)^/. Basta quindi ricavare £f utilizzando la conservazione del 
momento angolare. Poiche Ji = Jf si ottiene 

Ji = s v + s d + £i = 1 = Jf = s n + s n + £ f = S + £ f (2.6) 

La funzione d'onda totale dello stato finale ip = (px deve essere antisimmetrica per la 
statistica di Fermi-Dirac : gli unici stati possibili permessi da questa statistica sono 

1. 5 = 0,^ = 0,2,4... 

2. S = 1 , £ = 1,3,5... 

L'unica combinazione compatibile con Jf e data da 5 = 1 , £f = 1 , quindi H/ = — 1 = 

IV 

Produzione I pioni possono essere generati da scattering di nucleoni 

1 pione • p + p — > p + p + 7T° 

• p + p — > p + n + 7r + 

• p + n — > p + n + 7T° 

• p + n — > p + p + 7r~ 

2 pioni • p + p — > p + p + 7T° + 7T° 

• p + p >■ p + p + 7T + + 7T~ 

• p + P > p + 71 + 7T + + 7T° 

• p + p — >■ n + n + 7r + + 7r + 

Aumentando l'energia dei fasci si possono produrrre anche piu pioni. Vale la conser- 
vazione del numero barionico , per la quale il numero dei nucleoni si conserva. 



2.2 Risonanze 
2.2.1 Risonanze A 

Lo studio sperimentale dei processi di scattering pione-nucleone riserva alcune peculiarita 
fondamentali. La prima particolarita sussiste nell'andamento della sezione d'urto in funzione 
dell'energia del processo elastico (calcolata nel sistema del laboratorio) tt + + p — > 7r + + p : 
si nota infatti che vi sono dei picchi di risonanza che possono essere interpretati come stati 
eccitati metastabili. 
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II grafico della sezione d'urto in funzione dell'energia pud essere infatti fittato con una Loren- 
ziana , 

dove r rappresenta la larghezza a meta altezza. Come gia detto l'aumento di a intorno a 
~ 200 MeV puo essere interpretato come la formazione di uno stato intermedio che decadra 
poco dopo, ovvcro 

7T + + p > A + + > 7T + + p 

Dunque il pione ed il protone formano, interagendo , uno stato legato non stazionario molto 
instabile. Visto che T ha le dimensioni di una lunghezza si puo stimare il tempo medio di vita 
utilizzano il principio di indeterminazione 

t, hc 200MeV • fm n „ , n 2 o ,„ „, 

IY A ~ ft r A ~ -„„ , ~ 0, 6 ■ 10~ 23 s (2.8) 

Tc 115MeV-3- 10 23 fm/s ' v ; 

Cambiando le particelle dello scattering si ottengono gli altri tipi di particelle A + ,A°,A~. 
Ad esempio con il processo di scattering tt~ +p — > p + tt~~ l'andamento della sezione d'urto e 
lo stesso ma stavolta la particella associata alio stato metastabile e data da A . Le particelle 
A ++ ,A + , A , A - rappresentano i 4 tipi di risonanza possibili , tutti alia stessa energia. Da 
questo si ricava che la particella A ha isospin T = 3/2 e di conseguenza puo trovarsi nei 
diversi stati di proiezione degli isospin T3 = 3/2, 1/2, —1/2, —3/2 relativi rispettivamente alle 
particelle A ++ ,A+,A°,A" . Per queste particelle J n = (3/2)+. Le risonanze di tipo N* 
invece avvengono ad energie molto piu alte , intorno ai 1440 MeV, 1520 MeV : per queste 
J n = (1/2)+ e T = 1/2. Queste risonanze vengono indicate con l'asterisco perche fisicamente 
si possono interpretare come l'assorbimento di un pione da parte del nucleone che diventa 
quindi eccitato. Introduciamo ora il concetto di massa invariante , per la risonanza A. Si 
consideri quindi il processo ir + + p — > A ++ . II quadrimpulso relativo alia particelle A ha 
componenti 

J 1 



Pa = [ — .pa 1 
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La massa a riposo si calcola come 



P . ^ — r2 j n 2 _ n 2 _ PA C _|_ ( m AC ) _ ^2 



Consideriamo ora p^=Pn+Pp, ovvero = E n + E p , p& = p v + p n . In questo caso la 
massa a riposo e data da 

(.Ep + -EV) 2 , ^ , 2 22 

fL ~^ (Pp + Pit) = m A C 

ovvero tua = M{n + p) e la massa invariante. Questa espressione vale in ogni sistema di 
riferimento , dunque in particolare anche in quello del laboratorio. In quest'ultimo sistema si 
ha E p = m p c 2 e p p = dunque 



(m p c 2 + E n ) 2 - p\<? = m\c A 
L'energia E n e calcolabile con un semplice bilancio energetico 

El = {p^cf + (m^c 2 ) 2 ==>- E n =T n + m^c 2 
Visto che T w ~ 200 MeV si ottiene m A c 2 ~ 1232 MeV. 



(2.9) 



(2.10) 



2.2.2 Fotoproduzione dei pioni 

La reazione di fotoproduzione dei pioni e data da 

7+p^p + 7T° (2.11) 
La sezione d'urto del processo nel sistema del laboratorio ha l'andamento dato dalla figura 

1 




Si noti che la risonanza e presente intorno ai 340 MeV. Facendo questo tipo di scattering si 
possono ottenere risonanze mesoniche che verranno descritte nel seguito. 
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2.2.3 Risonanze p 

La risonanza di tipo p si incontra in processi del tipo ir + + p — > ir + + 7r° + p : in questo caso 
lo stato metastabile e dato dalla presenza della particella p + quindi la reazione puo essere 
riscritta come tt + + p — > p + + p — > ir + + tt° + p. Facendo un grafico della sezione d'urto in 
funzione della massa invariante M(7r + 7r°) si nota che e presente un picco, dovuto alia risonanza 
, intorno a 770 MeV. 




200 400 600 800 1000 1200 I'OO 



Analogamente a quanto visto per la A esistono 3 diversi tipi di particelle p , dati da p + , p°, p~ . 
Fittando con una Lorenziana si ottengono i valori T = 150 MeV e la massa di risonanza 
m p = 770 MeV. A titolo di esempio si riporta anche la reazione che porta alia formazione 
della particella p° 

~f + p->p°+p->ir + +ir~+p (2.12) 

Visto che questa particella esiste in 3 diversi stati di carica , si puo pensare come tripletto di 
isospin, quindi T = 1 e T3 = —1,0, 1. Inoltre si nota che J n = 1~. Si noti che , rispetto al 
pione, la particella p ha J / e massa diversa. 



2.2.4 Risonanze r] e u 

Si ottengono in processi del tipo 

7T+ +p -)• 7T+ +p + 7T + +7T~ + 7T° (2.13) 

In questo caso l'energia e sufficientemente alta per la formazione di 3 pioni. Se si considera la 
massa invariante per il sistema M(7r + 7r _ 7r°) e si grafica la sezione d'urto in funzione di questa 
si evidenziano due diversi picchi : il primo, corrispondente alle risonanze 7] , tra 500 e 600 
MeV ; il secondo , corrispondente alle risonanze co , in prossimita di 800 MeV. 
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In particolare dal fit si ottengono i valori di m n = 549 MeV e T v = 0,83 KeV ; J n = (T. 
Si osserva che T v e molto piccola rispetto al caso della particella A, quindi la vita media 
in quest'ultimo caso e molto piu alta : ~ 10~ 18 s. La risonanza decade quindi in modo 
elettromagnetico. Per le risonanze uj il fit restituisce m w = 783 MeV e r w = 9, 9 MeV mentre 
J n = 1~ . Le risonanze rj si possono ottenere an che in processi del tipo 



p + p -)• 7T° + 2-7T+ + 2tt~ 



(2.14) 



In questo caso il grafico della massa invariante del sistema, M(tt + tt -k ) contiene due picchi 



analogamente al caso iniziale (2.13). Viene spontaneo chiedersi se si possa scegliere qualsiasi 
combinazione di pioni per ottenere le stesse risonanze : in realta Tunica combinazione che pre- 
senta risonanza e quella analizzata con M(7r + 7r _ 7r°) ; altre combinazioni , come M(ir + ir + TT ) 
non danno luogo ad alcuna risonanza. Se ne conclude che ,per questo tipo di risonanze, esiste 
solo lo stato con carica nulla , dato dalla combinazione 7r + 7r~7r°, dunque lo stato sara sempre 
un singoletto di isospin, ovvero T = 0. Le proprieta viste si possono ricapitolare in una tabella. 



T 



m 



IT 


Or 


l 


vr+,vr u 


7T 


140 


V 


Or 





V 




550 


p 


1~ 


l 




P~ 


770 


LO 


1- 





LO 




783 



2.2.5 Interazioni mediante scambio di mesoni 

Per introdurre l'argomento successivo vogliamo ora discutere l'effetto di una risonanza nell'in- 
terazione nucleone-nucleone (vnn)- 

OPE=One Pion Exchange E gia stato analizzato il caso dell'interazione OPE. Un inter- 
azione di questo tipo si puo schematizzare con la figura seguente. 
inserire figura 



TPE=Two Pion Exchange Consideriamo quindi un processo di interazione tra due nucle- 
oni N. Se il primo nucleone e eccitato allora al suo interno si ha un processo che libera 
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un pione : il pione puo , a sua volta, eccitare il secondo nucleone che da luogo ad una 
particella A. Generalizzando si ha che i due nucleoni possono essere entrambi particelle 
A : in questo caso l'interazione e di tipo TPE , Two Pion Exchange. 

inserire figura 

p,uj meson exchange inserire figura In realta i tipi di risonanze, p, oj e i diversi meccanismi 
TPE , OPE , non sono arbitrari ma dipendono dalla distanza relativa tra i due nucleoni 
: graficando infatti il potenziale vnn in funzione della distanza relativa r si ha che , per 
distanze piccole si hanno piu frequentemente risonanze di tipo p, uj mentre , al crescere 
di r , si incontrano rispettivamente TPE e OPE. 

Interazioni tra 3 nucleoni Si puo ampliare il ragionamento esposto prima anche ad un 
sistema di 3 nucleoni legati da i potenziali v\2 , t>23 : in alcuni casi il problema puo 
essere ricondotto a due sottoproblemi di interazione a due corpi , mentre in altri casi 
la schematizzazione non e sufficiente. Per questi ultimi , quando si ha risonanza A nel 
nucleo , l'Hamiltoniana deve essere corretta , aggiungendo un termine di interazione che 
tiene conto dei 3 nucleoni. In questo caso si ha quindi 

v£j ijk 

Correzioni di ordine ulteriore migliorano i risultati sperimentalj^] Le due situazioni 
elencate precendentemente sono raccolte in figura inserire figura 



2.3 Mesoni e Barioni Strani 

Lo studio degli urti pione-nucleone porto alia creazione di particelle con proprieta particolari 
che vennero considerate inizialmente strane , visto che non rispettavano le leggi scoperte fino 
a quel momento. 

2.3.1 Kaoni 

II primo esempio di mesone strano e dato dal Kaone. Per questo esistono due stati di carica 
K + ,K° , con le masse = 493,7 MeV e m^o = 497,7 MeV. La somiglianza delle due 

masse porta ad intepretare i due stati di carica come stati di isospin della stessa particella , 
quindi T = l/2eT3 = ±1/2 ; J n = 0~. I mesoni K si possono produrre con il processo 

tT + p^n + K + + K~ (2.16) 

dove K~ = K + . Misurando la sezione d'urto di questo processo si ottiene a ~ 10~ 3 barn 
: questo valore e caratteristico delle interazioni forti. Nonostante questo si osserva speri- 
mentalmente che queste particelle riescono a sopravivvere per un tempo abbastanza lungo , 
tk = 1,2 • 10~ 8 s , quindi il processo responsabile del decadimento e un' interazione debole : 
ecco la stranezza. Si riportano i processi di decadimento. 

K + p+ + u„, (63%) ; K + -> tt+ + vr°, (21%) (2.17) 
1 Per ottenere un buon risultato si utilizzano correzioni fino al 18 mo ordine 
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La risonanza di energia piu bassa del Kaone e data dal K*(892) , T = 1/2 , J n = 1". Tale 
risonanza si forma con una reazione del tipo 

K + K* + K + + tt° + p (2.18) 

Questa particella decade forte ed ha una vita media molto piccola. 



2.3.2 Barione A 

Un tipo di barione strancj^je dato dalla particella A. Per questa m\ = 1116 MeV ; J n = l/2 + 
: possiede un solo stato di carica , T = 0. La particella si puo formare attraverso il processo 



7T 



+ n^A + K + (2.19) 



Osservando la sezione d'urto di questo processo si osserva che l'interazione e forte ma la vita 
media e di r ~ 2, 6 • 10~ 10 s quindi decade attraverso un'interazione debole. I modi principali 
di decadimento sono dati da 

A^p + 7T~ A^ra + vr (2.20) 

Decadimenti di tipo elettromagnetico per questo tipo di particella non sono mai stati osservati; 
similmente non e possibile ottenere un decadimento del tipo A — > p + K~ visto che Q < per 
questa reazione. Per risolvere i problemi legati a queste particelle venne introdotto un nuovo 
numero quantico , chiamato stranezza , che si indica con S : si assunse che le particelle come 
n,p, A, N* , 7r, p,oj,rj avessero 5 = 0. Per i mesoni e bosoni strani si assunse invece 5^0. 
In particolare si assegno S = 1 al K + . In seguito si postulo che la stranezza totale di un 
sistema si conserva nei processi in cui intervengono interazioni forte o elettromagnetiche , non 
si conserva invece per le interazioni deboli. Visto che il processo di produzione del K + e del 
K~ e dato da un interazione forte si ha S(K~) = —1 mentr^S(K°) = +1. Si noti che questa 
particella equivale ad un neutrone di stranezza diversa da 0. Ci si puo chiedere in che forma 
si possa scrivere l'interazione di questa particella con i nucleoni. L'interazione e forte , e a 
corto range, quindi lo scattering si puo descrivere in funzione dei soli parametri a (scattering 
length) e ro (raggio efficace). I parametri caratteristici sono 

[a= (-1,8 ±0,2)fm 

< v ' (2.21) 

|r = (3,16±0,52)fm 

Si osservi che l'interazione A — pe simile a quella N-N anche se piu debole e senza stati legati. 
Data la difhcolta di produrre fasci di particelle A polarizzati, la lunghezza di scattering del 
processo tiene conto di una somma di pesi statistici del singoletto e del tripletto. La particella 
mediatrice di questa interazione non puo essere un pione altrimenti la stranezza , nel processo 
A — t- p + ir~ non si conserverebbe. Si deve supporre che la particella mediatrice sia un Kaone 
, in modo che la stranezza si conservi. Si puo schematizzare la situazione in figura. 



2 I barioni strani vengono detti iperoni 

3 Questa ultima relazione e facilmente dimostrabile se considera il processo 7r~ + p — > p + K° + K~ , di 
interazione forte. II bilancio delle stranezze si puo scrivere come + = + 1 — 1 , dunque la stranezza di K° 
deve essere necessariamente 1. 
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P A 




A P 

Si noti che un processo di questo tipo non potrebbe avvenire spontaneamente vito che Q < 
, ma in questo caso il K~ e una particella virtuale che puo quindi violare la conservazione 
dell'energia in una frazione di tempo At compatibile con l'indeterminazione energia-tempo. 
Un altro modo di studiare l'interazione NA consiste nell'eccitare un nucleo in modo che al suo 
interno si formi una particella A : un nucleo di questo genere e detto ipernuclec^ 

Si noti che nel secondo caso l'interazione e mediata da 2 pioni, quindi si tratta di TPE ( 
Two Pion Exchange) . 



2.3.3 Barione S 

Di questo barione esistono 3 diversi tipi : 

• E+, m s + = 1189,4 MeV 

• E°, m s o = 1192,6 MeV 

• XT, m s _ = 1177,5 MeV 

Tali diversi tipi si possono interpretare come tripletto di Isospin , quindi T = 1 e J n = l/2 + . 
Queste particelle possono essere prodotte attraverso il processo di scattering 

ir- + p E~ + K + (2.22) 

Un altro processo possibile e 

K~+P^Z + + tt- (2.23) 

I decadimenti possibili sono E~ — > n + it con rg- = 1,5- 10~ 10 s e S + — > p + 7r°(— > n + 7r + ) 
con T£+ = 0, 8 • 10~ 10 s. Per quanto riguarda invece S° il processo di produzione e dato 
da K~ + p — > S° + 7T° : decade con il processo S° — > A + 7 ed ha una vita media di 
r s o =0,7- 10~ 19 s. Visto che in questo processo si conserva la stranezza e visto il tempo di 
decadimento l'interazione di decadimento sara molto probabilmente eletromagnetica. 



2.3.4 Barione E 

Le sue caratteristiche sono m- = 1318 MeV e T = 1/2 ; J n = 1/2+ . Esistono due stati di 
carica H° con m s o = 1315 MeV e H~ con massa mg- = 1321 MeV. Questo particolare barione 
viene prodotto tramite il processo 

K~ + p -»• 3- + K + (2.24) 

Se si calcola il bilancio delle stranezze si ottiene che deve essere necessariamente <S(3~) = —2 : 
in effetti la sezione d'urto dimostra che e un processo forte dunque questo tipo di ragionamento 
e applicabile. La particella 3 decade debolmente attraverso i decadimenti 

4 Una descrizione accurata degli ipernuclei e data nell'apposito capitolo. 

5 Si noti che al secondo termine non e presente p visto che si deve conservare il numero di barioni. 
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• ~° -»• A + vr° con r H o = 2, 9 • l(T 10 s 

• E~ ->• A + 7T- con r s - =1,6- 10" 10 s 

Si noti che non avviene il decadimento H° — > S° + 7r° poichc Q < 0. 



2.3.5 Barione 

E il piu massivo di tutti , infatti tuq- = 1673 MeV. A differenza degli altri J n = 3/2 + ; T = 0. 
La stranezza si ricava facendo il bilancio sulla reazione di produzione 



k~ + p -> n~ + K + + K° 



(2.25) 



Si ricava quindi S(Q ) = —3. La vita media e di tq- =0,8-10 10 s mentre i decadimenti 
possibili sono dati da 

• Q- -»• A + K~ 

• Or -»• H° + 7T- 

• 0- -> H" + 7T° 

Visto che i processi sono deboli in questi la stranezza non si conserva ma cambia di un'unita 
, come si vede esplicitamente. Si noti che non avviene il processo VT — > E° + K~ poiche , in 
questo caso, Q < 0. 

2.3.6 Ottetti e decupletti Barionici 

Vediamo di riassumere i barioni che hanno la coppia spin-parita in comune. Per J n = l/2 + 
si ha la seguente tabella 



J 



n 



1/2+ 



5 = p,n (T = i) 

S = -l A (T = 0), (T = l) 

S = -2 H-,H° (T= 1 ) 



Tale configurazione e detta Ottetto Barionico. Queste particelle possono essere riprodotte 
anche in un grafico. 

3 
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Raggruppando i barioni di spin-parita J n = 3/2 + si ottiene la tabellsj^] 

J n = 3/2+ 



s = 





A++,A+,A°,A- 


(T = 3/2) 


s = 


-1 




(T=l) 


s = 


-2 




= 1/2) 


s = 


-3 


Q- (T = 


0) 



Anche per il decupletto si puo fare un grafico simile. 




Lo stesso metodo grafico si puo applicare anche ai mesoni, ottenendo 





7T+ Tq 



2.4 Ipernuclei 

Nella sezione relativa alia particella A sono gia state descritte le caratteristiche fisiche. Si 
vuole ora descrivere 1'interazione che essa stabilisce con un nucleone, eccitando un nucleo : 



3 Le risonanze di tipo E*,H* sono analoghe al caso visto in (2.181 
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tale sistema e detto ipernucleo. La reazione di eccitazione e data da 

K~ + A X^iX + ir- (2.26) 

dove A e il numero totale di barioni , ovvero A = N + Z + N\. In laboratorio si riescono a 
produrre ipernuclei con AeE. All'interno del nucleo avviene la reazione 

K-+n^A + ir- (2.27) 

Si noti che questo nucleo non e stabile ma decade con la reazione 

A p + tt- (n + ir°) (2.28) 

Supponiamo ora di aver prodotto un ipernucleo A dell'ossigeno con A = 17 : in questo 
caso si puo utilizzare il modello a shell osservando che nelPultimo livello energetico (W5/2) e 
presente una particella A. In questo caso non e phi accettabile la combinazione J n = 5/2 + 
per la presenza della particella A : e giusta la combinazione J A = l/2 + . Si consideri invece 
l'ipernucleo A dell'eho-4 : visto che J = j n + j\ si ha che J = 0, 1. Anche questo nucleo non 
e stabile ma decade in un atomo di elio-4 ed un pione ir°. La reazione e riportata nella figura 
seguente. 

1.09MeV 1 + 



O.O^V \ 4 ge+?r o 



0+ 



E anche possibile misurare l'energia di legame , B , della particella A nel sistema dell'ipernucleo 
utilizzando questo decadimento ottenendo B\{\He) = 2,33 MeV. Questa energia deve essere 
confrontata con quella del normale nucleo Bi^He) = 28 MeV ricordandosi pero di dividere 
per 4 visto che B\ rappresenta l'energia , in media, per ogni nucleone. Si ottiene quindi 

^E^L = 7MeV > B A (\He) (2.29) 

quindi la particella nell'iperone e meno legata. Misurando indirettamente la massa del nucleo 
si puo avere una stima dell'energia di legame della particella A utilizzando il bilancio 

B A = [m(A) - Mi^X) - M(iX)} c 2 (2.30) 



I risultati sperimentali per B\ sono raccolti nel seguente grafico. 
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Capitolo 3 

Classificazione delle particelle 

La divisione bosoni-fermioni si basa sullo spin e sulla parita della funzione d'onda. Esistono 
altri modi per classificare le particelle : si possono infatti dividere le particelle a secondo 
delPinterazione che sviluppano con altre. 

Adroni Risentono delPinterazione forte e delle altre. 

Barioni Adroni a spin semiintero 
Mesoni Adroni a spin intero 

Leptoni Non risentono delPinterazione forte 

Le masse dei Leptoni variano di molto e non seguono un andamento precise 

Particella ^ v e t~ u t 

M [MeV] 0,511 105,7 <2 < 0, 27 1784 < 35 ' 

Si elencano ora le varie proprieta caratteristiche di queste particelle. 

3.1 Proprieta generali 

3.1.1 Parita 

La parita si conserva in processi con interazione forte , ed elettromagnetici ma non con la 
debole 

3.1.2 Numero barionico 

Si assegna a tutti i barioni il numero barionico B = +1 mentre agli altri B = — 1. In 
particolare ai mesoni si assegna B = , cosf come ai Leptoni. In base a tutti i processi 
osservati sperimentalmente si nota che B e sempre conservato. Questa conservazione si nota 
in processi del tipo 

• 7T + + p — > 7T + + 7T° + p 

• p + p^-p + p+p + p 

• e + + e~ — > p + p 

Infine la legge di conservazione del numero barionico permette di affermare che non possono 
avvenire processi del tipo p — > e + + 7. 
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3.1.3 Numero Leptonico 

Analogamente a quanto visto per il numero barionico esiste una legge di conservazione anche 
per il numero leptonico. Questo numero viene definito , a differenza di quello barionico, per una 
famiglia di leptoni. Esiste quindi il numero Leptonico elettronico L e , il numero Leptonico 
muonico L M ed il numero Leptonico tauonico L T . Ovviamente all'elettrone corrisponde il 
tripletto (L e , L^, L T ) = (1, 0, 0). I risultati di questa classificazione sono riportati nella tabella 
seguente. 



L e L^ L T 



e 


Ve 


1 








e+ 


Ve 


-1 










Vp 





1 












-1 





T 


V T 








1 


T+ 


V T 








-1 



Questa tabella rappresenta inoltre un valido strumento per bilanciare le reazioni con neutrini 
o antineutrini. Ad esempio nella reazione n — > p + e~ + 9 e e giusto utilizzare l'antineutrino. 
Analogamente si nota che il processo e~ + e~ — > ir~ + ir + non puo avvenire perche non 
rispetterebbe la conservazione del numero leptonico. Ha senso invece un processo del tipo 
e~ + e~ — > ir~ + 7r + . Si noti che la conservazione del numero leptonico va fatta separatamente 
per ogni tipo di numero leptonico: a titolo di esempio si consideri la reazione fi~ — > e~ + u e + u^. 
La conservazione del numero Leptonico elettronico daO = +1 — 1 + mentre quella del 
numero Leptonico muonico da +1 = + + 1. Vediamo ora come sia possibile distinguere 
sperimentalmente le particelle v e e v e . A questo scopo si consideri il decadimento f3 dato da 

2x N ^2 +1 Y + e- + v e (3.1) 

Si supponga poi di prendere le particelle v e e di utilizzarle per una nuova reazione v e + p — > 
n + e + . La stessa cosa si puo fare per v e provenienti dalla reazione e~ +p — >■ n+v e '■ stavolta il 
processo e dato da v e +p — > n + e + . Questo processo e poco probabile quindi si nota che le due 
particelle v e e v e forniscono risultati sperimentali diversi tra loro, dunque sono distinguibili. 



3.1.4 Isospin 

Si supponga di aver fissato T. Se si indica con Q la carica elettrica e con Q il suo valor medio 
allora l'isospin si puo scrivere come 

Mi"!) -;(«-«) < 32 » 

Ad esempio per T = 3/2 si hanno le particelle A++, A+, A , A". Si calcola Q tot /e = 2 mentre 
Qtot/e = 2/4 = 1/2 e si verifica che 

T 3 = (--- J =2-1/2 = 3/2 (3.3) 
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3.2 Modello a quark degli adroni 

Questo modello fu proposto all'inizio degli anni '60 : si suppose che gli adroni fossero com- 
posti dai quark e che questi ultimi avessero S = 1/2. Tale ipotesi e necessaria visto che si 
deve spiegare l'esistenza di adroni che hanno anche spin semiintero, dunque il costituente el- 
ementare deve avere spin semiintero. I quark non esistono liberi , ma solo in stati legati che 
formano i mesoni , q±q2 , ed i barioni , qiq2q3- Se si indicano poi i quark con la notazione 
u,d,s (up, down, strange) il numero di accoppiamenti possibili si puo calcolare con il numero 
di combinazioni senza ripetizioni : 



a 



n,k 



(n + k-l)\ 
fc!(n- 1)! 



C 3 ,3 = 10 



(3.4) 



Risulta quindi conveniente cercare di assegnare ai quark un numero barionico. Visto che un 
barione e composto da 3 quark segue che B q = 1/3. Per i mesoni invece possiamo supporre 
che essi siano composti da q\,q~2 , dove q\ rappresenta un antiquark. Supponiamo ora di voler 
determinare la carica elettrica dei quark : a tale scopo si consideri il decupletto barionico 
che ha carica totale nulla. Si puo calcolare la carica dei quark ragionando in questo modo 
: supponiamo , ad esempio , che la A ++ sia formata dai 3 quark , A ++ = uuu ; da questa 
si ricava Q u /e = 2/3. Analogamente si ricava A~ = ddd =>■ Qd/e = — 1/3 , fi~ = sss =>- 
Q s /e = — 1/3. Si noti che in questo modello si stanno introducendo particelle che hanno come 
carica complessiva un multiplo (non intero!) della carica delPelettrone. Riassumiamo quindi 
le principali caratteristiche dei quark nella tabella seguente. 



Flavour 


Spin 


B 


Q/e 


T 


T 3 


S 


n 


m(MeV) 


u 


1/2 


1/3 


2/3 


1/2 


1/2 





+i 


5 


d 


1/2 


1/3 


-1/3 


1/2 


-1/2 





+i 


10 


s 


1/2 


1/3 


-1/3 








-1 


+i 


150 



I quark possono essere riportati in un grafico (T3S) , come gia fatto per i mesoni e barioni , 
ottenendo il seguente. 




Si noti che , nonostante il quark non esista libero , e stato comunque possibile dare una 
stima della massa. Questo e stato possibile definendo due masse per il quark : la massa 
costituente e la massa di corrente ( misurata sperimentalmente in maniera indiretta ). Fino 
ad adesso non sono stati osservati quark isolati : questa particolaria viene attribuita alia natura 
delbinterazione quark-quark. Questa e , schematicamente, approssimabile come crescente con 
la distanza quindi per allontanare due quark dovrei fornire una quantita sempre piu crescente 
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di energia. Ad un certo punto avro fornito talmente tanta energia da formare una coppia qq 
dei quali, uno andra a rimpiazzare quello si sta tentando di estrarre mentre il rimanente andra 
a formare un mesone (\qq)) con quello che si sta estraendo , non permettendo l'isolamento. 
Tutte le particelle principali viste fino ad ora si possono scrivere in funzione dei quarkQ : 

• p = uud 

• n = udd 

• A = uds 

• 7T + = ud 

• 7r~ = ud = ud 

• 7T° = uu, dd 

Si noti in particolare che 7r° = uu, dd = ir°. Per i Kaoni si ha che 

• K + = us 

• K° = ds 

• K + = us = K~ 

• K0 = ds^K° 

Si considerino ora le particelle del tripletto di isospin T = 1. 

f7r+,7r ,vr- J n = 0" ,m 7r = 140MeV . . 

T = 1-S = 0=*>{ J ' n ' 3.5 

\p + ,p ,p- J u = l~ ,m p = 770MeV 

Si vuole capire come si spiega la differenza dei due stati formati dai quark ir + = ud ; p + = ud 
: in effetti queste due particelle sono composte dalla stessa combinazione di quark ma hanno 



masse completamente differenti , come visto in (3.5). II pione e Padrone piu leggero quindi si 
puo supporre che la coppia ud si trovi nello stato con £ piu piccolo possibile, ovvero I = e lo 
stato ha parita negativa. Con J = si puo comporre solo S = 0. Analogamente supponiamo 
per il mesone p + che sia £ = e dunque la parita sia nuovamente negativa. Per avere J = 1 
stavolta deve essere S = 1. Dunque i due stati sono fondamentalmente diversi , visto che il 
primo (it) rappresenta un'onda 1 5o mentre il secondo (p + ) un'onda 3 S\ : questi due livellj^] 
sono giustamente spaziati energeticamente. Si noti che anche in questa interazione si nota 
una dipendenza dallo spin. 

3.3 Massa dei quark 

Si riportano le masse dei quark gia visti nella sezione precedent^ 
• u = up =^ m u = 1,5-^3 MeV 



1 Si noti che ty° 7^ ss in virtu della conservazione dell'isospin , infatti T(ir°) = 1 e T(ss) = 0. 
2 Sperimentalmente si osservano anche gli stati con £=1,2; S = 0, 1 : la risonanza 61 a 1323 MeV e data 
dallo stato £ = 1, S = 0, J n = 1+. 

3 Si noti che il simbolo ■£■ non indica la divisione ma deve essere letto come tra. . . e. . . . 
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• d = down ==> m<f = 3-^7 MeV 

• s = strange ==>• m s = 95 ± 25 MeV 

• c = charm =^> m c = 1, 27 ± 0, 1 GeV; Q = \e 



b = bottom => m b = 4, 20 ± 0, 1 GeV; Q = -A 



e 



i = top =>• m t = 175 ± 7 GeV; Q 



In particolare le masse appena illustrate rappresentano la massa corrente dei quark. In effetti 
la massa dei quark si puo interpretare in due modi 

Massa Corrente : massa di un quark libero. Visto che i quark non esistono isolati, non ha 
molto senso parlare di massa corrente. Si pensa che, comprimendo piii adroni, si possa 
arrivare ad una situazione, detta di libertd asintotica , nella quale i quark si comportano 
come liberi : tale materia e detta quark- gluon-plasma. Questo stato e raggiungibile solo 
con altissime energie ( non possibili in laboratorio ) o ad altissime densita , ad esempio 
nelle stelle di neutroni. 

Massa Costituente : massa efficace che presentano i quark quando sono legati e formano 
gli adroni. Ad esempio si puo considerare la massa dei quark u e d come 1/3 della massa 
di un nucleone. Infatti , visto che m p ~ m n , m u ~ nid ~ 171^/3. 



3.4 Colore dei quark 

Si consideri la particella A ++ = uuu; J n = 3/2 + . Questa rappresenta Padrone piii leggero 
tra quelli con J = 3/2 quindi e plausibile assumere che t = e che lo stato sia combinazione 
degli spin | ftt)- Con queste ipotesi la funzione d'onda radiale , cosi come quella di spin, 
e simmetrica rispetto agli scambi di particella : questo non e in accordo con l'ipotesi che la 
particella A sia un fermion^] Inoltre non e rispettato il principio di Pauli visto che i quark si 
trovano tutti nello stesso stato quantico . Per ovviare a questo problema si introduce un nuovo 
numero quantico , detto colore , che puo variare tra G (green), R(red) e -B(6/we|^| Ovviamente 
esistono anche gli opposti dei colori , R,G,B : il colore va pensato come una estensione della 
carica con la differenza che , nel caso elettromagnetico , esiste solo la carica e l'anticarica 
mentre , nel caso dei colori, esistono 3 tipi di colori e i relativi opposti. Si noti che gli adroni 
non mostrano caratteristiche di colore : in questi i colori si mescolano quando i quark vanno 
a formare la particella in modo da avere un colore totale risultante nullo=BIANCO. Valgono 
le proprieta 

• B+R+G=0 

• B+B=R+R=G+G=0 
La funzione d'onda complessiva di un adrone si scrive quindi come 

V' = V'spin' i /'spazio' i /'flavourV ; colore (3-6) 



4 Questi hanno, per definizione, funzione d'onda totale antisimmetrica. 

5 Per distinguere i quark di diversi colori verra apposto un apice : ur,ug,ub- 
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La parte relativa al colore cambia nei barioni e nei mesoni. 

f -4= \RGB + GBR + - RBG - GRB - BGR] barioni 

^colore = < V r - (3.7) 

\ ^ + GG + 55] mesoni V 7 

L'esistenza della A ++ a 1232 MeV con J n = 3/2 + puo essere spiegata solo con l'aggiunta di 
questo numero quantico. Si noti che il colore di un quark , analogamente ad altre grandezze 
fisiche, puo variare nei tempo. 

3.4.1 Interazioni tra Quark 

La particella mediatrice delle interazioni tra quark e il gluone , particella senza massa e con 
spin 1. II campo che lega i quark e il campo di colore : si vede quindi che il colore non e 
soltanto un numero quantico correttivo , ma rappresenta la fonte dell'interazione forte tra i 
quark ( analogo della carica per l'interazione E.M.). II gluone deve portare due tipi di carica 
colorata : 1 colore e 1 anticolore. FIGURE Si noti che 

1. Le interazioni forti NON possono cambiare il flavour di un quark. Possono solo avvenire 
riarrangiamenti di quark nelle particelle. 

2. Coppie quark-antiquark possono essere create e annichilate. 

3. Le interazioni deboli possono cambiare il flavour dei quark , mediante i bosoni vettori 
W ± e Z°. Si ha che mw — 80 GeV/c 2 , mz — 90 GeV/c 2 . I processi in cui i quark 
cambiano flavour si possono scrivere come 

u^d + W + ; u^d + W~ (3.8) 
s^u + W~ ; s^u + W + (3.9) 

Ad esempio nei processo n — > p + e~ + uno dei quark d che compongono il neutrone 
cambia flavour con la reazione d — > u + W~ , ottenuta dalla prima delle precedenti 
equazioni mediante lo scambio di W + . 




3.4.2 Prove dell'esistenza dei quark e del loro colore 

Prendiamo uno scattering e + e~ (elettrone-positrone) ad alte energie. NelPimagine seguente si 
osservano sciami adroni. Questo si puo spiegare osservando il processo dal punto di vista dei 
quark. La reazione e del tipo e + + e~ — > 7 — > q + q : la coppia quark-antiquark , allontandosi 
in direzione opposta , aumenta l'energia creando coppie quark-antiquark che ricombinandosi 
creano adroni. 
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Figure 18.1T (a) A two-jet event, in which q and q each produce a shower of 
hadrons. (o) A relatively less-common three-|et event, in which one ol the quarks 
radiates a gtuon which also creates a shower of hadrons. The visible tracks show 
the trajectories ot the hadrons, bent into arcs by a magentie field Experimental 
results are from DESY. 



L'esistenza del colore venne invece provata confrontando le sezioni d'urto dei due processi 
e + + e~ —7- Adroni e e + + e~ — > u + + II risultato sperimentale di (3.10)) puo infatti 
essere spiegato solo con l'aggiunta del colore come nuovo grado di liberta . Teoricamente si 
otterrebbe infatti un risultato tre volte piu piccolo : l'errore sta nel considerare i soli 3 gradi 
di liberta relativi ai flavour dei quark mentre aggiungendo i colori si arriva a 9. 



R 



cr(e + + e — > Adroni) 
cr(e+ + e~ — > u + u~) 



~ 2 



(3.10) 
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Capitolo 4 

Decadimenti a 



L'emissione spontanea di una particella a e rappresentata dal seguente processed 

zXN^zZiYN-2 + <* (4.1) 

II decadimento a e un effetto della repulsione Coulombiana : in effetti quest 'ultima diventa 
significativa per nuclei pesanti in quanto , mentre la forza Coulombiana cresce con Z 2 , la 
forza di legame dei nuclei cresce approssimativamente come A. L'emissione della particella 
a si spiega osservando che V\He e un nucleo doppiamente magico, quindi molto legato, ed 
il bilanciamento energetico e per questo molto favorevole.Si vogliono cercare i nuclei avvan- 
taggiati per questo tipo di decadimento studiando il segno di Q : se Q > la reazione sara 
avvantaggiata, non avverra altrimenti. La conservazione dell'energia impone 

M x c 2 = M Y <? + T Y + m a c 2 + T a (4.2) 

E quindi facile ricavare il Q di reazion^] 

Q = M x - M Y - M a = T a + T Y (4.3) 

II decadimento avviene spontaneamente solo se Q > 0. Dalla conservazione dell'impulso segue 
che il nucleo Y e la particella a si muovono in direzioni opposte con impulso uguale in modulo 
'■Pa = Py- Sperimentalmente si osserva che per Y ed a vale T <C mc 2 . Si pud quindi usare la 
cinematica non relativistica ponendo T = p 2 /2m per queste particelle ricavando 

+ Q -T. + 7>-4-(>±±*)+T.- TT * (4.4) 




2m a \ M Y J 1 + M a /M Y 

II valore di T a e misurabile ( ad esempio con un spettrometro magnetico) quindi il valore di Q 
e calcolabile utilizzando la formula precedente. I Q calcolati con questo metodo sono riportati 
in tabella : si noti che il primo valore di Q > corrisponde alia particella a. 

Particellaemessa Q (MeV) 

n -7,26 
p -6,12 



A He 5,41 
5 JIe -2,59 

1 La particella a rappresenta un nucleo di \He. 
2 Si ponga c = 1. 
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Procedendo con nuclei ancora phi massivi si trova che nuclei come 9 Be e 12 C hanno Q > 0. 
Studiando i nuclei massivi che decadevano a venne messa in evidenza un'importante proprieta 
: la vita media variava di molti ordini di grandezza a seconda del Q . Si consideri ad esempio 
232 Th con Q = 4,08 MeV ed un altro isotopo 218 Th con Q = 9,85 MeV. Nonostante i due 
Q differiscano di un fattore 2 si osserva che le due vite medie variano da r = 1,4 • 10 10 yr 
per il primo a r = 1,07 • 10~ 7 s per il secondo. Questo comportamento si puo spiegare con 
Pintroduzione della regola di Geiger-Nuttal. 

Q ~ 7^- (4-5) 

L'andamento sperimentale del tempo di dimezzamento e ben visibile nella seguente figura. 




I I 1 i I I 

4 5 6 7 8 9 10 

Q (MeV) 

In generale si osserva che , per variazioni di Q = (4 4- 10) MeV il tempo di vita varia tra 
10 -9 s < r < 10 17 yr. Per stabilire il segno di Q si puo utilizzare la formula semi-empirica di 
massa visto nel primo modulo. Si tratta quindi di calcolare 

Q = M X -M Y - M a 

conoscendo 



M x = Zm p + Nm n - B(Z, A) 
M Y = (Z - 2)m p + (N - 2)m n - B(Z - 2, A - 4) 

M a = 2m p + 2m n — B a 

L'espressione per Q si semplifica 

Q = B a + B(Z -2, A -A)- B{Z, A) 
dove B a = 28, 3 MeV. Utilizzando la formula semiempirica di massa 



B(Z,A) = a v A-a s A 2/3 



ac 



Z 2 

ZV3 



'■sym 



(A - 2Z f 
~A 



+ 6 



per il primo ed il secondo nucleo e sostituendo nella(4.9) si ottiene 
Q = B a + a s \a 2 ^ 3 -(A- 4) 2 / 3 l 



+ a c 



z 2 A -i/3 _ ( Z -2) 2 (A-4)- 1 



-1/3 



(4.6) 
(4.7) 
(4.8) 

(4.9) 
(4.10) 
(4.11) 
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Visto che A>4( stiamo considerando nuclei pesanti) si puo sviluppare il termine di superficie 
( e tutti gli altri termini ) 



1-1 + 



3A 



= -asA-V* 



(4.12) 



as ^2/3 _ {A _ 4) 2/3 ^^2/3 

In generale si utilizza lo sviluppo 

(i4 _ 4) ._^( 1 _*)' =! v( 1 -^) 

Utilizzando lo stesso procedimento per gli altri termini si ottiene infine (stiamo considerando 
nuclei pari-pari) 



(4.13) 



Q(Z, A) -28, 3 - 4a v + jja^" 1 / 3 + 4a c ZA-^ (l - 

(2Z\ 2 
+3a p A~ 7 / 4 



(4.14) 



Questa espressione per Q(A, Z) definisce, a Q fissato , delle curve nel piano Z — N : la curva 
divide le due regioni dello spazio dove Q > o Q < , ovvero indica quali sono gli elementi 
stabili e quelli che tendono a decadere a. 



4.1 Modello di Gamow-Gurney-Condon 

Venne introdotto per spiegare l'andamento della vita media in funzione di Q. Nel modello si 
suppone che la particella a sia gia preformata nel nucleo : la probability di uscire dal nucleo e 
pari alia probability di penetrazione di una barriera di potenziale ed e lagata alia vita media 
del decadimento. II potenziale e di tipo Coulombiano 




(4.15) 



A r 
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Per dare una stima di a si puo sommare il raggio nucleare a quelle- della particella a , a ~ 
Ry + R a = (6 -T- 8) fm . La profondita della buca e data invece da Vq = (35 40) MeV. Visto 
che Q = (4 -=-8) MeV per questi tipi di reazione, manca da stimare solo l'altezza della barriera 



B = V(r = a) 



Z Zy ( 



4vT£o 



\4lT£ohc ) 



Z a Zy 



ahc 



2Z Y _ 200 Z a Z Y 
137 a 



(4.16) 



Scegliamo a = 7 fm e Zy = 90 . Per questi valori si ottiene B ~ 37, 5 MeV. In generale i valori 
critici di B si trovano nel range B = (30 40) MeV, ovvero Q < B per ogni reazione, quindi 
classicamente la particella non potrebbe mai passare : il passaggio e dovuto all'effetto tunnel. 
Indichiamo con A la probability di decadimento, A = 1/r. In questo modello si suppone 
che A sia proporzionale alia probability di penetrazione della barriera, moltiplicata per il 
numero di tentativi che la particella fa per passare, ovvero la frequenza con cui la particella 
tenta di penetrare la barriera. Questa si puo stimare con / = v a /a dove v a si ottiene dalla 
conservazione dell'energia 



-mva = Q + Vq 



2(Q + V ) 



f 



2(Q + Vq) c 
a 



(4.17) 



2 V mc" V mc* 

Con una stima si ricava (Q = 5 MeV, Vq = 35 MeV, a = (6 8) fm) v a ~ 0, 15 • c =>- 
/ = 10 21 Hz. Se definiamo poi con b il punto di inversione classica del moto si puo stimare la 
larghezza della barriera : 

Si vuole calcolare or a la probability di penetrazione. Per fare questo si considera prima una 
barriera phi semplice , rappresentata in hgura. 



V(r = b) = Q=>b = ahc- 



b ~ 44 fm (b - a) ~ (35 H- 40) fm 



(4.18) 



V 








E 




Ar 






■4 K 





La probability di penetrazione per questa barriera e nota ed e data da 

1 , J2m(V - E) 

P = - 2 con k b = ^ y - '- 4.19 

Conviene poi sviluppare 

e x -e~ x . , 2 e 2x + e~ 2x -2 

smhx = smh z x = (4.20) 

2 4 v ' 

Sappiamo che m a ~ 4000 MeV , V ~ 35 MeV, Q ~ 5 MeV , quindi sostituendo si trova 
fcft ~ 2,4 fm -1 . Del resto Ar ~ (b — a) ~ 30 fm , quindi kbAr ~ 52 » 1, quindi il termine 
dominante al denominatore e solo Pesponenziale crescente. Si ottiene quindi 



(4.21) 
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a meno di un fattore costante che non si discosta molto da 1. Per dare una stima piu re- 
alistica bisognerebbe risolvere l'equazione di Schroedinger per il potenziale realistico definite- 



dalla (4.15). II calcolo si puo fare in maniera piuttosto esatta utilizzando l'approssimazione 
semiclassica di WKE0 . II risultato del calcolo e dato da 



p _ V(fl-Q)(Q-Vb) _ 2G 
B-V 



(4.22) 



dove 



G 



2m 



[V(r)-Q] 1/2 dr 



(4.23) 



e il fattore di Gamow. Si deve quindi svolgere Pintegrale J a ° [V(r) — Q] 1 ^ 2 dr . Sviluppiamo 



V(r) -Q = ahc 2 ^- - Q = Q ( b - 1 
r V r 



(4.24) 



Si definisce ora la variabile adimensionale x = r/b dr = bdx quindi l'espressione si 
trasforma in 

1/2 

dx, con xq = a/b = Q/B (4.25) 



G= V^Q b 



Si sostituisce ora x = cos 2 ■ 
L'integrale diventa quindi 





"1 








- 1 




X 




\[x 




COS 



dx = 2 cos (if — sin( 



, (p = arccosy/xo 



G 



^/2m~Qb 

H, 

s/2rnQb 

h 

yfMQb 

H, 





cos 2 <p 


J<t>o \ cos2 <P 


j"^ sine? 


-2 cos q 




Jo cos (, 


9 


r<t>o 




/ sin 2 


^d0 


Jo 




[00 - sin 


4>o cos <; 


F(x ) 





2 cos ip(— sin (/?) d(/> 



(4.26) 



II risultato dell'integrale si scrive quindi come 

1/2 



F(x ) 





"i-l" 


f 




>XQ 





dx 



arccos x. 



f-xV 2 (l-xo)V 2 



(4.27) 



3 In meccanica quantistica l'approssimazione WBK (Wentzel-Kramers-Brillouin) conosciuta anche come ap- 
prossimazione WKBJ (Wentzel-Kramers-Brillouin-Jeffreys) e un'approssimazione semiclassica nella quale si 
impone che la funzione d'onda sia una funzione esponenziale che varia lentamente e quindi viene sviluppata in 
serie di potenze della costante di Planck. 
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Per x < 1 si hf0F(x o ) ~ vr/2 - 2x l J 2 = vr/2 - 2 y / Q/B. Visto che A = fP = r~ l , e che 
P oc e~ 2G , si ricava 



In r = In 



a / mr 



c V 2(Q + V&) 



2G 



Basta esplicitare G ricordandosi della (4.27) e della seguente. 



G 



-bF(x ) = — b 



7T Q 

- -2\/ — 
2 VB 



Combinando queste due equazioni si ricava quindi 



a 



Z a Z Y V2 



vr 1 1 



In r = In 



a / mc^ 



c V 2(Q + Vb) 



+ 2aZ a Z y \/2 



vr 1 /l 

27^" Vs 



(4.30) 



(4.31) 



(4.32) 



Facendo un grafico di questa funzione si puo studiare l'andamento della vita media relativa- 
mente ai parametri tipici del sistema che si sta considerando. 



4.2 Processi favoriti 

Si consideri il nucleo 220 Th : di questo sono stati osservati due decadimenti : 

• 220 Th^ 216 Ra+a con Q a = 8, 95 MeV e r a = 9, 7 us 

• IfTh^f Po+ 12 C con Q c = 32, 14 MeV e r c = 3, 3 • 10 6 s. 

Per capire quale dei due decadimenti sia favorito si possono elencare le differenze sostanziali 
tra i due: 

1. La barrier a Coulombiana cambia nei due casi : nel primc^] si ha B a ~ 2-88 = 176 
MeV mentre per il secondo Be ~ 6 • 84 = 504. Quindi e avvantaggiata l'emissione della 
particella a. 

2. Si nota che Qc > Qa quindi sembrerebbe avvantaggiata l'emissione del Carbonio. 

3. E stata fatta l'ipotesi che la particella a fosse preformata nel nucleo : tale ipotesi e 
sempre valida nella seconda reazione, a patto per di intendere che la probabilita di tale 
preformazione e sicuramente minore di 1. Infatti e meno probabile che nel nucleo si 
preformi il Carbonio-12 piuttosto che una partiella a. Matematicamente questa osser- 
vazione si concretizza nell'aggiunta di un termine w nell'espressione della vita media 
: A = fP • wq- Ovviamente wq <C w a quindi, nuovamente , e favorito il processo di 
emissione della particella a. 

4 Si sono utlizzati gli sviluppi 

7T 1 3 

arccosa; ~ x x + . . . (4.28) 

2 6 y J 

(l-z) 1/2 ^l-f + ... (4.29) 

arrestandosi al primo ordine. 
5 Si ricordi che B oc Z a Zy- 
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4.3 Struttura fine degli spettri a 

II decadimento a e un processo a 2 corpi quindi l'energia della particelle a e ben determinata 
dalla conservazione dell'energia. Ci si aspetterebbe quindi che E a si mantenga costante al 
variare di N a . In realta l'energia cambia , come si vede nella figura seguente nella quale e 
rappresentato lo spettro a relativo al decadimento del 25l Fm 




-.-of : < *-»'ev 



Questo comportamento si puo spiegare considerando che il nucleo figlio Y si trova in uno stato 
eccitato : il passaggio alio stato con E a minore rappresenta la diseccitazione del nucleo che 
passa ad un altro stato energetico. II processo e quindi del tipo 

X^Y* + a^Y + -f + a (4.33) 

La particella a , essendo formata da 2 neutroni e 2 protoni tutti nello stato Is e con gli spin 
accoppiati a 0, avr momento angolare totale dovuto solo al momento angolare orbitale. 

4.3.1 Momento angolare e parita 

Si consideri ora la reazione standard dei decadimenti a,X->F + «esi supponga che 
il nucleo X si trovi in uno stato J n = + . II nucleo Y puo decadere in vari stati en- 
ergetici + , 2 + , 4 + , 2~, 3 + : per valutare quale decadimento sia favorito si puo utilizzare la 
conservazione del momento angolare e della parita. Applicando queste condizioni si ha 

Conservazione momento angolare : ij = If + £ a , visto che I a = 0. Allora £ a = {If + 
ii), ... ,\If — Ii\ . In particolare se I{ = =>■ l a = If. 

Conservazione parita : Hj = IL/(— l) ia . Si distinguono due casi 

• ITj = Ily ==> £ a pari 

• Hi = —Hf ==> £ a dispari 
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Si vede quindi che sono permesse le transizioni + — > + , + — > 2 + , + — > 4 + con momento 
angolare rispettivamente £ a = 0,2,4 mentre non sono concesse le transizioni + — > 2~ , + — > 
3 + . Se si osservano le intensita dei diversi stati si nota che quelli ad energia piu grande ( gli 
stati eccitati ) hanno intensita sempre minori. Questo si puo spiegare con due effetti 

1. Se Qo e il Q della reazione quando il nucleo figlio decade nello stato fondamentale ed E* 
e l'energia dello stato eccitato allora Q* = Qo — E* e dunque Q* < Qo , quindi lo stato 
eccitato ha una vita media piu lunga e dunque ha meno probability. 

2. Nel considerare la probability di attraversamento della barriera abbiamo supposto £ = 
, quindi il termine di potenziale centrifuge 6 - ^mr^ non era stato considerate Questo 
termine aggiunge al pontenziale un termine efficace che , per I ^ innalza la barriera 
di potenziale, diminuendo cosf la probabilita di tunnelling. 

Infine se si studia il decadimento del 238 Pu si nota che le energie dei livelli 8 + , 6 + , 4 + , 2 + , + 
stanno nei seguenti rapporti tra di loro 

Ea Er Er Es , . 

= 3, 3; = 2, 07; = 6, 82; = 1, 68 (4.35) 

2 ft4 ft2 ft6 

Questa separazione e dovuta alia presenza di uno spettro rotazionale. Le bande di questo 
spettro (bande rotazionali) hanno energia 

*=^t!> (4 .36, 

dove I rappresenta il momento di inerzia. II rapporto tra i vari livelli energetici e quindi del 
tipo 

E n n(n + 1) 



E m m(m + 1) 
questa formula e in buon accordo con i dati sperimentali. 



(4.37) 



6 Si ricordi che le parte radiale dell'equazione si Schrodinger per un campo centrale e data da 

R"+*R'+ f k 2 -^^)R = (4.34) 



Capitolo 5 

Fusione Nucleare 



£ il meccanismo attraverso cui le stelle producono energie, inoltre e alia base di creazione di 
molti degli elementi presenti sulla terra. 



He CO NeMgSi Fe 

I — i — i — i I i — I rj i 1 1 1 r 




MASS NUMBERS A 



Osservando la figura precedente si nota che , per estrarre energia dal nucleo , si possono far 
fondere tra loro nuclei leggeri , con A piccolo : questo significa che vengono utilizzati nuclei 
nella parte sinistra della curva relativa alia binding energy. II processo di fusione si schematizza 
con la reazione 

a + X^Y + b (5.1) 

II Q di reazione si scrive come Q = m a + mx — fny — mb '■ ovviamente sono possibili le reazioni 
con Q > 0. La reazione di fusione piu semplice p+p — )- 2 He non avviene a causa dell'instabilita 
dell' 2 He. Le reazioni piu probabili si dividono in vari tipi 

Reazioni p-p E un esempio la reazione 

p + p -> 2 H + e + + v e con Q = 1,44 MeV (5.2) 
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CAPITOLO 5. FUSIONE NUCLEARE 



La sezione d'urto di questo processo e molto piu piccola rispetto al normale scattering 
pp e l'interazione e debole. La barriera Coulombiana per questi tipi di reazione si puo 
calcolare come 



Br, 



1 



~ 900KeV 



pp 47re 2Rp 

Reazioni D-D Vengono fatti fondere due atomi di deuterio. 

2 H + 2 H — > 4 He + 7 con Q = 23,8 MeV 



(5.3) 



(5.4) 



Perche si conservino l'energia e l'impulso c'e bisogno di un gap energetico aggiuntivo , 
rispetto ai processi p-p , fornito dal fotone. Risutano piu probabili processi del tipo 

• 2 H+ 2 H^ 3 He+n con Q = 3,3 MeV 

• 2 H+ 2 H^ 3 H+p con Q = 4 MeV 

Inoltre si calcola .Bdd ~ 360 KeV. 
Reazioni D-T Vengono fatti fondere un atomo di deuterio ed uno di trizio. 

2 H + 3 H — » 4 He + n ; Q = 17, 6 MeV (5.5) 

Nel grafico seguente vengono riportate le sezioni d'urto , al variare dell'energia per i tipi di 
processi appena elencati. 
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5.0.2 Temperatura 

Volendo avere reazioni termonucleari di fusione, ovvero reazioni in cui l'energie dei nuclei di 
partenza e solo di agitazione termica, quabe la temperature necessaria per attivare il processo 
di fusione ? 

fe~ k B T 

\k B ^8,6- l(T 5 eV/K 



5.1. SEZIONE D'URTO 
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Se voglio, ad esempio, ksT ~ 1 KeV serve una temperatura di T ~ 10 7 K. Nel sole si ha 
Tq = 1,57 • 10 7 K , che corrisponde ad un'energia di 1,35 KeV. Si noti che non serve avere 
un'energia termica maggiore della barriera Coulombiana visto che questa puo essere passata 
per effetto tunnel. 



5.1 Sezione d'urto 



La sezione d'urto per la fusione nucleare si puo definire a partire dalla sezione d'urto definita 
per le interazioni nucleari e la probability di penetrazione della barriera Coulombiana. Si 
ottiene quindi 

CTfus = O-nucl • P (5.7) 

Nel caso di scattering NN a basse energie e gia stato visto che si puo considerare il solo 
contributo di t = ottenendo 



o(l = 0) 



1 



sin 2 6 (E) 



(5.8) 



quindi <j ^ \ J hp 1 ^ l/^.Poniciiiio quindi cxnuci 

= S{E)/E , dove E e detto fattore astrofisico. 

Si ha quindi 



<7fus = E e 



(5.9) 



dove G e il fattore di Gamow definito nella (4.23) : conviene approssimare questa formula per 
il particolare caso della fusione. 



„ V2mE, 
G ~ — - — b 
n 



7T 

2 



(5.10) 



Innanzitutto si consideri il rapporto E/B : visto che le grandezze caratteristiche 
900 KeV , B DD ~ 350 KeV t^E ~ k B T Q ~ 1,35 KeV il rapporto \jEjB e trascurabile 
rispetto a ir/2 quindi 



„ V2mE7r u 
G ~ — b 

h 2 

Utilizzando la definizione di (4.16) G si puo riscriver^] come 

/ 2mc 2 vr Z a Z x 
G ~ aZ a Z x \ — — — = vtq — - — 
V E 2 u/c 

Si e usata anche la relazione E ~ mv 2 /2. 



(5.11) 



(5.13) 



5.2 Reaction Rate 

Per stimare la potenza prodotta da un reattore a fusione nucleare si deve conoscere il tasso 
di reazione che , a differenza della sezione d'urto, dipendera anche dalla concentrazione delle 
specie che si considerano. Si definisce quindi 

,. reason. 
Volume • lempo 



1 Nel seguito verra adottato il simbolo per indicare una quantita relativa al sole. 
2 infatti 



Z « Z * = akc^L (5.12) 



4tt£o E v/c 
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Per calcolare questo rapporto ci si riferisce sempre a processi del tipo a + X — > b + Y . Si 
supponga quindi di inviare un fascio di particelle a , a velocita v a , verso un bersaglio X. Si 
indica con rtj la densita delle due specie di particelle , ovvero 

n x = -~ ; n a = ^ ( 5 - 15 ) 

Si puo quindi scrivere 

R = nx<y^a = nx — -• — = n a nx(JVa (5.16) 
s-t v a 

dove a e la sezione d'urto di fusione del processo. D'altronde la sezione d'urto e anche definita 
come (vd. Appunti del primo modulo) 

Rb (5.17) 



I a N X 



dove I a rappresenta il numero di particelle a inviate (per unit di tempo) , Rb il numero di 
particelle b prodotte (sempre nellunita di tempo) e con Nx il rapporto tra il numero di nuclei 
bersaglio X e la superficie F dove si fa la misura. Ovviamente tutte queste quantita hanno 
senso sperimentalmente solo se si considerano le rispettive medie : per la velocita occorrera 
la distribuzione di Maxwell-Boltzmann data da 



f m \ 



3/2 



hm{v) = 47r (^ 27r fc B r J v e 2fcsT ' dovev = v a -v x (5.18) 
Quindi 

/•oo 

(R) = n a n x (crv) = n a n x I afMB{v)&v (5.19) 



per a ^ X. Nel caso particolare in cui le particelle che interagiscono sono della stessa specie 
si ottiene 

(R) = \n 2 a (av) (5.20) 
Si devono quindi studiare le due funzioni integrande nel termine (av) . La sezione d'urto e 



data da ( per la ( |5.9| ) e la (5.13)) 
S(E) 



a 



E 



/ 27TaZ a Z x \ S(E) f ,f^- 27 7 1 1 ,, , 

exp < ; > = — — — exp < — Trav2mc L Z a Zx—F= } (5.21 

{ v/c j E { Je) 



Se si definisce poi /3 = ira\/2mc 2 Z a Z x =± 31, 29(KeV) 1 / 2 V2 Z a Z x , dove m u = 931, 494 



MeV e l'unita di massa atomica , la (5.21) diventa semplicement^] 



a= S[El e -P/VE (5 _ 22) 
E 



La distribuzione di Maxwell-Boltzmann puo essere scritta in funzione dell'energia relativa con 
un semplice cambiamento di variabili : /mb oc Ee~ E l kBT . Se si studia Pandamento delle due 
funzioni ed il loro prodotto si ottiene il grafico seguente 



3 Esiste anche la notazione equivalente e 27rr >^ E ^ con rj(E) 



v/c 



5.2. REACTION RATE 



47 




Si vede quindi che la funzione integranda della (5.19) e trascurabile al di fuori di una piccola 
zona dove il prodotto delle due funzioni raggiunge un massimo ( detto Picco di Gamow ) . Alia 
reazione contribuiscono quindi solo le particelle appartenenti alia coda della distribuzione di 
Boltzmann . La zona del picco di Gamow corrisponde ad un range di energie molto basse, 
dimcilmente raggiungibili in laboratorio. Per ottenere l'andamento a queste energie si cerca 
quindi di ricavare l'andamento sperimentale ad alte energie in modo da trovare la curva che 
interpola bene i dati e ricavare l'andamento anche per basse energie , sotto i 100 KeV, del 
fattore S(E). 




Nell'esperimento L.U.N. A. del 1998 sono state raggiunte energie molto basse ed e stato 
osservato che il fattore astrofisico tende nuovamente ad aumentare 
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questo effetto e dovuto alia presenza della nube elettronica attorno ai nuclei ( screening). 
Raggiungendo energie sempre piu basse si osserva che S(E) si mantiene pressoche costante : 
utilizzeremo questa proprieta nel capitolo successive 



5.2.1 Proprieta del picco di Gamow 

Si vuole ora stimare l'energia del picco di Gamow in funzione della carica e della massa ridotta 
dei nuclei. Per fare questo si consideri la forma del rate medio 

1/2 / 1 \ 3/2 

MP, 



(<jv) 



(- 



i \ 3 / 2 roo ( 

5r) /„ S(E) ^{ 



E 



k B T 



dE 



(5.23) 



Per trovare il massimo si pud massimizzare la funzione integranda ( supponendo S(E) per 
quanto detto prima ) 



_d_ 
dE 



exp 



E 
~k~B~f 







E=E 



Con semplici calcoli si ottiene 

'Ph B T\ 2/3 



E n 



1.22-[^^ 

TYlqt, 



o 



1/3 



KeV 



(5.24) 



(5.25) 



dove si e utilizzata la notazione Tq = T/10 6 K. Dopo aver ricavato il massimo si pud ottenere 
anche la larghezza utilizzando il metodo di Laplace 



exp 



E 

kef 



C exp 



E — Eq 
A/2 



Le costanti C e A si calcolano facilmente fittando la curva con una Gaussiana 
quindi 



A = -(E k B T) 



1/2 



E = (10 -r 30) KeV; A = (10 20) KeV 



(5.26) 
si ottiene 
(5.27) 



da confrontare con le energie tipiche k B T = (1 -j- 3) KeV. 



4 Si supponga di dover calcolare 7 = /_ °° g(x) e F ^ dx , dove g(x) dipende debolmente da x . Nel pun- 
to di massimo xq si puo sviluppare F(x) = F(xo) + F'(xo)(x — xo) + (l/2)F" (xo)(x — xq) 2 — F(xo) + 
(1/2) + F"(xo)(x — xq) 2 . Data la debole dipendenza da x di g Pintegrale si puo quindi riscrivere come 
I ~ g{xo) e~ F ( x °^ exp ^—F"(xo) ^ c ~% ^ j da;. Quindi si puo approssimare e~ F< - x ^ come una Gaussiana. 



5.3. RISONANZE BELLA FUSIONE 



49 



5.3 Risonanze della fusione 

Si prenda il processo p+ 12 C — ^ 13 N + 7, che avviene nel sole. Osservando l'andamento della 
sezione d'urto di fusione in funzione dell'energia si nota una deviazione dal comportamento 
standard , che consiste in un picco a ~ 0, 5 MeV. Si interpreta questo picco come una risonanza, 
analogamente a quanto fatto per le A : questa particolare risonanza rappresenta uno stato 
eccitato metastabile del 13 iV. 



5.4 Nascita degli elementi 

Se si osserva un grafico dell'abbondanza relativa agli elementi costituenti il sistema solare si 
osserva come l'idrogeno sia presente per il 90 % e l'elio per circa il 9%. 




MASS NUMBER A 

Questa abbondanza relativa pud essere spiegata con la teoria del Big Bang. Si teorizza che 
nello stato primordiale la materia fosse altamente densa e calda : queste caratteristiche fisiche 
la costrinsero a reagire per formare nuovi elementi. AM'inizio i quark erano separati ( non 
esistevano solo in stati legati ) e formavano il cosidetto stato di quark- gluon-plasma ( QGP 
) ; in seguito i quark si legarono per formare adroni e diedero luogo ad elettroni, protoni e 
neutroni. In questo stadio la temperatura era di circa T ~ 7, 5 • 10 9 K. La densita era ancora 
molto alta quindi il rate di reazione era elevato e sufficiente a superare la barriera Coulombiana 
: tramite fusione si crearono i nuovi elemementi 4 He, 2 H, 3 He, 7 Li. Man mano che l'universo si 
andava espandendo la temperatura e la densita diminuivano quindi la barriera Coulombiana 
si alzava e non c'erano piu le condizioni per creare nuovi elementi tramite fusione. Tutti gli 
altri elementi sono stati prodotti dalle stelle , dove la concentrazione era abbastanza elevata 
da poter dar luogo a nuove reazioni di fusione , e portati nell'universo grazie alia loro presenza 
nella nube primordiale. 
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5.5 Catena pp - Processi di fusione protone/protone 

Questi sono i principali processi che avvengono nel sole, del quale ricordiamo le caratteristiche 
fisiche. 

• M Q = 1,989 • 10 33 g 

• Rq = 6, 96 • 10 5 km 

• T| ore = 1,57- 10 7 K 

• Tq P ~ 6000 K 

• p™ re ~ 130 g/cm 3 

• L Q = 3, 846 • 10 33 erg/s 
II processo 

1 H+ 1 H-K 2 H + e+ + i/ e (5.28) 

e gia stato analizzato. Dato che il positrone tende a reagire con un altro elettrone il Qtot 
della reazione va calcolato sommando quello di entrambi i processi : per la prima reazione si 
ha Q = 0,42 MeV mentre per la seconda Q e + e - = 1,02 MeV dunque in totale Q to t = 1,44 
MeV. La sezione d'urto di questo processo e molto piccola visto che l'interazione responsabile 
e quella debole : prendendo infatti energie tipiche di 1 KeV si ha er(~ 1 IKeV) — 10 -33 
Barn. Tuttavia il sole tramite questi processi e in grado di produrre molta energia visto che 
ha molti protoni a disposizione. Per stimare quindi il rate delle reazioni nel sole prendiamo 
innanzitutto il rate per protone ( per la 5.20| ) : 

rat6 = = l nJav) pp ~ 5 • 10- 18 reaz/sec (5.29) 

n. protoni n p 2 

Basta quindi moltiplicare questo risultato per il numero di protoni per ottenere il rate totale. 
Per ricavare il numero di protoni si puo fare una stima grossolana dividendo la massa del sole 
per la massa dei protoni : l'ordine di grandezza del risultato e corretto. 



M & 2 • 10 33 r 7 

K- — = -\ a 1 n 24 - 10 P r0t0ni 5 - 30 

m p 1, o • 10 



Ovviamente non tutti questi protoni sono disponibili per la fusione : si puo supporre che solo 
1/10 di questi abbia le condizioni termiche opportune per reagire. Dunque 



A/; burn ~ ^M p ~ 10 56 => AA rcaz = 5 • 10~ 18 • 10 56 ~ 10 38 reaz/sec (5.31) 



Quindi nonostante la sezione d'urto sia molto bassa il rate di reazione e molto grande. Si puo 
stimare anche la massa dell'idrogeno consumata al secondo 

~ 10 14 g/sec ~ 10 8 ton/sec (5.32) 



Riprendiamo ora in esame il processo di produzione della (5.28). 



1. Per quanto riguarda 2 H sono probabili i seguenti processi 



5.5. CATENA PP - PROCESSI DI FUSIONE PROTONE/PROTONE 
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(a) 2 H+ 1 H-> 3 H+e + + v e : processo debole, a piccola. 

(b) 2 H+ 2 H— » 4 He+7 : abbondanza del deuterio scarsa , quindi probability molto bassa 
che la reazione avvenga. 

(c) 2 H+ 1 H— » 3 He+7 : e la phi avvantaggiata. 

(d) 2 H+ 2 H— > 3 B.+p : abbondanza del deuterio scarsa , quindi probability molto bassa 
che la reazione avvenga. 

Si sceglie quindi la reazione 



2. Per quanto riguarda 3 He sono possibili i seguenti processi 

(a) 3 He+ 1 H— ^ 4 Li+7 non avviene perche il litio non e in uno stato legato. 

(b) 3 He+ 1 H^ 4 Li^ 2 H+2 1 H non e possibile poiche Q < 0. 

(c) 3 He+ 1 H-» 4 He+e + + v e viene detto processo hep ed ha un Q = 19, 28 MeV. E rile- 
vante in quanto, nonostante la sezione d'urto molto bassa ( dovuta ad un processo 
debole ) , genera i neutrini piu energetici e quindi maggiormente rilevabili. 

(d) 3 He+ 2 H— > 5 Li non e possibile perche il litio non si trova in uno stato legato 

(e) 3 He+ 3 He^ 4 He+2 1 H con Q = 12, 82 MeV. 

II processo piu avvantaggiato e quindi 



2 H + 1 H -»• 3 He + 7 



Q = 5, 49MeV 



(5.33) 



3 He + 3 He -»• He + 2 X H 



34) 



Ricapitolando dunque la catena di processi piu probabile nel sole e data da 



i 



L H + 1 H ^ 2 H + e+ + v e 
2 H + 1 H ^ 3 He + 7 
3 He + 3 He ^ 4 He + 2 X H 



(5 
(5 
(5 



35) 
36) 
37) 



II bilancio totale della reazione e dato da 



Particelle 4 l H He + 2e+ + 2u e 



Energia Q tot = 26, 72 MeV 



Questo ramo di reazioni , il piu probabile, viene chiamato ppl. Esistono anche i rami ppll, 
pplll : per avere un'idea della divisione di questi processi si faccia riferimento alPimmagine 
seguente. 
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Jh + 1h 




^He + ^He -» *He + 2 [H 
(PPI) 



Be + e~ 



^He + *He 



99.7% 




B + y 



(pp id 



^B *Be + e + + i/ e 



(PPJI1) 



Si osservi die il bilancio delle particelle , cosf come (piu o meno) quello energetico , e lo stesso 
nei tre cicli. Visto che nelle reazioni sono coinvolti neutrini occorre tener conto dell'energia 
irradiata da questi : si definisce quindi un Q e ff che tenga conto di questa perdita di energia. 
Risulta quindi naturale porre Q e ff = Qtot — ^neutrini- Nel caso dei rami analizzati prima si 
ottiene 

Processo Q c g [MeV] 



ppl 


26,2 


ppll 


25,66 


pplll 


19,17 



L'avanzamento della reazione dipende , oltre che dalla concentrazione dei reagenti, dalla 
temperatura. Questa dipendenza e ben visibile se si cerca di ricavare la variazione della 
concentrazione di deuterio utilizzando 



dt 



n d = R, 



VP 



R 



1 



■pd 



n2 p { av )pv ~ n p n d {av) T 



(5.38) 



Le medie (erf) infatti sono funzioni della temperatura T. L'equazione (5.38) e risolvibile 
ricavando prima la soluzione in condizioni stazionarie : all'equilibrio la variazione sara nulla. 
Dopo aver illustrato tutte le caratteristiche microscopiche dei processi che avvengono nel sole 
e nelle stelle si possono fare alcune interessanti stime. Ad esempio si puo calcolare la binding 
energy del sole e vedere se e trascurabile rispetto all'energia totale. L'energia di legame dei 
protoni e data da 

B & = N PQ m p - M G& (5.39) 
In particolare, per il sole si ha che Mq & (? ~ 2 • 10 33 g • c 2 dove lg • c 2 = 9 • 10 20 erg. Si ottiene 



quindi M Q c z = 1,8 • 10 erg. La binding energy puo essere scritta come 



D 



3 GMq 
5~R^ 



1 3 2GM Q M Q c 2 3 .Rsch 



25 



R P 



10 Rr. 



M c 2 







(5.40) 



5.6. CICLO CNO 
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ove i?sch rappresenta il raggio di Schwarzschilcj^JPer il sole B ~ 10 34 erg , da confrontare con 
l'energia totale ~ 10 erg. Nelle stelle di neutroni il rapporto tra queste due energie e di 1 a 10 
: la stella nella compattazione libera una quantita di energia molto piu alta di quella del sole. 
Utilizzando i dati a disposizione si puo fare una stima del tempo di vita solare supponendo 
che l'energia venga prodotta con il meccanismo illustrate E gia stato visto che nel ciclo ppl si 
ha Q ~ 26, 7 MeV. Inoltre e noto il numero di protoni disponibili per la fusione, ovvero 1/10 



della (5.30). Per calcolare l'energia liberata dal sole basta moltiplicare il numero di protoni 



disponibili per il Q di reazione ottenendo 

^"burn 

Etot = - 1051 er g ( 5 - 41 ) 

si noti che il 4 al denominatore e stato aggiunto visto che 4 protoni si uniscono per formare 
una particella a. Infine per stimare la vita del sole si puo calcolare 

*0 = = o o 1 ^ 6 / - 2, 6 • 10 17 s * 8 • 10 9 yr (5.42) 
L Q 3, 8 • W 66 erg/s 



5.6 Ciclo CNO 



Un altro tipo di ciclo presente nelle stelle e quello di tipo CNO : in questo caso la stella deve 
avere a disposizione carbonio quindi deve trattarsi di una stella di seconda generazione. II 
ciclo complessivo e mostrato in figursj^] 



15 N+ p + ^12 C+ 4 2He * 





-* 15 N+p +^160+y 



12 C+p + ^ 13 N+y 








13 N 


-^i3C+e + +v e 



160 + p+^17 F+y 


17p 






17 F 


:^i7o+ e + +v e 



13 C+p + ^ 14 N+7 - 



14h 



N+p + ^ 15 0+y - 17 0+p + ^ 14 N+ 4 2 He 



II. 



Si noti che il 1 2C agisce da catalizzatore per la reazione perche viene sempre riutilizzato : se 
questo non e presenta la reazione non puo iniziare. II bilancio totale del processo da 



4 : H 



He + 2e + + 2v e 



(5.43) 



5 II raggio di Schwarzschild rappresenta il raggio critico per il quale la velocita di fuga relativa ad un corpo 
diventa c : oggetti con r < r sc h sono definiti buchi neri visto che neanche la luce riesce a sfuggire airimmenso 
campo gravitazionale prodotto dalla loro supermassivita. 

6 Con la notazione p + in figura si intende, ovviamente , il protone, ovvero equivalentemente 1 H come indicato 
fin'ora. 
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II ciclo CNO sembrerebbe favorito rispetto al ppl visto che la prima reazione di quest'ultimo e 
un processo debole. Nonostante questo il ciclo CNO fornisce solo una piccola parte dell'energia 
proveniente dal sole : per risolvere questo dubbio basti osservare che la barriera Coulombiana 
nel primo processo e molto piu alta dunque la reazione e meno favorita. II ciclo CNO e 
favorito solo a temperature basse ; dopo una certa temperatura il ciclo ppl diventa dominante 
, come si osserva bene nella figura seguente. 




5 10 15 20 25 30 

(771 6 K) 



Queste predizioni teoriche si possono veriflcare sperimentalmente studiando i neutrini prove- 
nienti dal sole. Infatti questi non interagiscono con la materia e conservano le informazioni 
riguardanti il processo che li ha generati. Proprio per questo motivo tali esperimenti devono 
essere condotti sotto terra , alia profondita maggiore possibile, in modo da separare gli eventi 
buoni da quelli non rilevanti sperimentalmente. 



5.7 Modello Solare Standard 

II modello solare standard rappresenta il miglior modello ad oggi disponibile per descrivere 
il comportamento fluidodinamico , chimico e fisico del sole. Analizziamo schematicamente le 
varie equazioni che compongono questo modello. 

1. Inizialmente bisogna scrivere l'equazione di equilibrio idrostatico per una distribuzione 
di materia a simmetria sferica non in rotazion^J 



dP _ n m(r)p(r) _ 

" (5.45) 



— —G r2 — > gradiente di pressione 
= 4:irp(r)r 2 — > distribuzione della massa 



Per risolvere questo sistema serve un'equazione di stato per il materiale stellare , ovvero 
bisogna conoscere la variazione di pressione in funzione della densita del materiale e delle 
coordinate spaziali. 



7 La prima equazione si ottiene con il bilancio 



dS[P(r) - P(r + dr)] - F g {r) = conF 9 (r) = G ™^^ — (5.44) 
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2. Visto che il corpo e composto da vari materiali bisognera conoscere la concentrazione 
di ognuno di essi. La seconda equazione per il modello e quindi data schematicamente 
dalla seguente 

P = P(p,X,Y,Z) (5.46) 

dove sono state definite X: abbondanza di X H ; Y: abbondanza di 4 He ; Z: abbondanza 
dei metalli. La pressione avra due contributi : uno dovuto alia materia stessa e uno che 
dipende dalla radiazione (oc T 4 per Stefan-Boltzmann). 

P = Pgas + ^rad = KpT + ^aT 4 (5.47) 

3. Bisogna inoltre conoscere la luminosita 

L q = e nud {r)4irr 2 p(r)dr (5.48) 
Jo 

e la sua legge di variazione 

= e n ucip(r)47rr 2 (5.49) 

dove e rappresenta la quantita di energia prodotta nel punto r per unita di massa e 
unita di tempo. Si osservi che la prima equazione rappresenta la condizione al contorno 
per la risoluzione della seconda. Ovviamente questa dipendera da 

£nucl = Enucl({<Tv) ij,p, X, Y, Z, T) (5.50) 

Si dovra scrivere un'equazione anche per il trasporto di energia che potra essere sia radiativo 
che convettivo. Risolvendo tutte queste equazioni si possono determinare gli spettri dei neu- 
trini emessi dal sole , verificando le ipotesi teoriche. Si nota che il calcolo di questi spettri e 
molto sensibile ad una piccola variazione di temperatura. Ad esempio nel caso della reazione 
8 B ^ 8 Be + e + + v con N v ~ T 13 . Quindi i neutrini possono essere utilizzati come termometro 
per studiare la temperatura del sole. 



5.8 Neutrini 

I neutrini sono di fondamentale importanza per lo studio dei processi che avvengono nel sole, 
come gia illustrato nel paragrafo precedente. Per avere un'idea del numero di neutrini prodotti 
dal sole che giungono sulla terra su una superficie di 1 cm 2 . II numero di neutrini prodotti al 
secondo dal sole e dato da 

M v = 2^ ~ 2 • 10 38 neutrini/sec (5.51) 

Si deve quindi calcolare il flusso di neutrini, definite come 

<& u = ~ ~~ttt ——— ~ 7 • 10 10 neutrini/cm 2 • sec (5.52) 
And 2 sec 

dove d e la distanza terra-sole , d = 1 U.A.= 1,496 • 10 13 cm. La misura di questo flusso 
non e semplice : alia difficolta di trovare un sito adatto per il rivelatore si aggiungono i tempi 
di attesa per verificare una reazione indotta dai neutrini. In figura e riportato lo spettro di 
neutrini e il tipo di rivelatore che e stato utilizzato per ricavare l'andamento sperimentale. 
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Neutrino Energy (MeV) 



5.8.1 Rivelazione dei neutrini 



Homestake Uno dei primi esperimenti fu quelle- di Homestake (Davis) (1967). II rivelatore 
utilizzato era al Cloro : C2CI4. II processo che permetteva di visualizzare i neutrini era 
il seguente v e + f £ CI — > f | Ar + e~ . Questo processo ha bisogno di neutrini abbastanza 
energetici con E u > 0,814 MeV. L'Argon prodotto decade (/3) successivamente con 



-1/2 



35 days. I risultati vennero riproposti con una nuova unita di misura , definita 



da 



lSNU(Solar Neutrino Unit 



lneutrino catturato 
10 36 nuclei bersaglio • secondo 



(5.53) 



I risultati dell'esperimento sono riportati in figura. 



2.0 



1.5 



a 



r< 



1.0- 



0.5- 



u 

a 
o 
a: 



* 1 - 



t h ^ f 

1122 :n » 11V .1? 



SNU 



t m 



4 TUI 



- 



1970 i 1971 i 1972 j 1973 1 1974 j 1975 1 1976 1 1977 1 1978 j i979 



10 



mentre nelle figure seguenti vengono riportate alcune foto dell'attrezzatura utilizzata per 
le rivelazioni 
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II valore previsto era di 8 SNU mentre nel grafico vennero rivelati valori molto piu bassi. 
In generale si ebbe che (/)^ p ^ — \4>^ e °^ ■ Inizialmente la discrepanza tra i valori speri- 
mentali e quelli teorici suggerf che il modello solare potesse risultare da un'elaborazione 
errata dei processi solari. In realta il modello venne confermato anche da esperimenti 
successivi, quindi la spiegazione doveva essere un'altra. 

KamiokaNDE/SuperKamiokaNDE La dicitura sta ad indicare che l'esperimento viene 
condotto nella miniera di Kamioka, mentre NDE rappresenta il tipo di esperimento : 
Nuclear Decay Experiment. In questo si riesce a rivelare anche la direzione dei neutrini 
e il rivelatore e ad acqua. In effetti viene sfruttato lo scattering elettronico 

v e + e~ — > e~ + v e (5.54) 

La vasca utilizzata, insieme al rivelatore, e ben visibile nella seguente foto. 
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In questo particolare scattering l'elettrone si muove ad una velocita maggiore della luce 
nel mezzcJ^J Viene quindi prodotta radiazione , detta radiazione Cerenkov. L'angolo 
caratteristico della radiazione definisce il cono nel quale viene emessa la radiazione ed e 
dato da 

(C / TLX 
(5.55) 
v J 

L'apparato e sensibile anche ai neutrini appartenenti alle altre famiglie v il -,v T , Anche 
in questo esperimento venne rivelata una discrepanza tra dati teorici e sperimentali : 

4 exp) ^|4 theo) . 

Gallex/GNO Vengono utilizzati rivelatori al Gallio , che hanno bisogno di neutrini con 
E v 3> 233 KeV. La reazione utilizzata nel rivelatore e la seguente 

v e + 71 Ga -> 71 Ge + e~ (5.56) 

SNO(Sudbury Neutrino Observatory) Viene utilizzato un rivelatore ad acqua pesante 
(D 2 0). 

5.9 Fasi finali dell'evoluzione stellare (cenni) 

Una volta esaurito l'idrogeno nelle zone centrali della stella, dove la temperatura e abbastanza 
alta da permetterne la fusione in elio si giunge alia situazione seguente : un nucleo ricco di 
4 He e rivestito da uno strato di idrogeno a temperatura minore. In questa configurazione non 
sono possibili nemmeno le reazioni di fusione dell'elio , quindi la stella , non avendo phi fonte 
di energia, si raffredda e si contrae. L'energia associata al nucleo della stella e data da 

2 

Tfl 

E G ~ —G — (5.57) 
r 



Quando questo si contrae la variazione di energia e data da 



Tfl 

AE G = G^Ar < (5.58) 



a questa variazione e associato un aumento delPenergia interna. Infatti, per il teorema del 
viriale 

(E kin ) = -^(E G )^AE kin = -^AE G (5.59) 

dunque l'energia cinetica aumenta , quindi aumenta anche l'energia interna del sistema. In 
questa fase la stella si riscalda , quindi si ricreano le condizioni per la fusione dell'idrogeno 
in una corona sferica immediatamente adiacente alia sezione dove e presente l'elio. Se la 
temperatura aumenta ancora inizia a fondere anche l'elio con conseguente produzione del 12 C. 
Dopo questa fase l'elio si esaurira , il nucleo si raffreddera nuovamente e ricomincera il ciclo 
di fusione con i prodotti di quello precedente. Ad ogni nuovo ciclo della stella i materiali 
utilizzati sono sempre phi pesanti; il limite e costituito dal ferro che rappresenta l'elemento 
phi legato in natura per il quale la fusione produce energia : gli elementi al limitaj sono 2§Ni 



8 In un mezzo con indice di rifrazione n la luce ha velocita c' = — < c. 

9 Viene considerato solo l'isotopo 56 Fe in virtu della sua abbondanza: in realta l'isotopo 58 Fe sarebbe piii 
legato del precedente ma la sua scarsa abbondanza ( ~ 0, 29%) fa si che questo non venga considerato nella 
trattazione della fusione stellare. 
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e 56 Fe. Man mano che il processo va avanti gli elementi piu pesanti si distribuiscono al centro 
della stella cosicche questa assume una configurazione a strati : andando verso l'esterno si 
trovano elementi sempre piu leggeri. Sorge spontanea la questione sul destino della stella 
quando questa non puo piu sintetizzare elementi. In questo caso la stella non ha piu energia 
dunque il nucleo continua a contrarsi fino a che questa implode. In particolare solo le stelle 
con M > 8Mq riescono ad arrivare fino alia sintesi del ferro, altre si arrestano prima. 




Irm,< _ 



Si supponga ora che una stella riesca a sintetizzare tutti gli elementi fino al Carbonio. L'arresto 
del processo di sintesi e dovuto alia pressione di degenerazione degli elettroni che e in grado 
di sostenere la stella. II limite di massa per cui la stella viene sostenuta dalla sola pressione 
di degenerazione e detta massa di Chandrasekhar : se la massa della stella e minore di questo 
valore allora e stabile , altrimenti collassa. 

5.9.1 Fusione dell'elio 

II processo e dato da 4 He+ 4 He-> 8 Be con r = 7 • 10 17 s. Prima che il Berillio decada 
spontaneamente c'e una probability non nulla che avvenga la reazione 



4 He + 8 Be ^ 12 C + 7 con Q = 7, 366MeV 



Successivamente 



4 He + 12 C O + 7 con Q = 7, 161MeV 
Si puo continuare aggiungendo altre particelle a 

4 He + 16 0^ 20 Ne + 7 

Quando la temperatura non e piu sufficiente per la fusione la reazione si arrestera. 



5.9.2 Fusione del Carbonio e Ossigeno 



T ~ 7 • 10 8 K 



12 



C+ 12 C 



T ~ 10 9 K 



24 Mg + 7 

12 C+ 12 C ^ 23 Na+ Ijj 
12 C +12 C ^ 20 Ne + 4 He 

j-16 0+ 16 ^ 32 g + 7 
16q +16 O _». 31 p + 1 H 

16 +16 o -». 28 Si + 4 He 



(5.60) 
(5.61) 

(5.62) 



(5.63) 



(5.64) 
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5.10 Gas di Fermi 

L'argomento e stato sviluppato nel primo modulo quindi la descrizione sara limitata all'ap- 
profondimento dei temi gia trattati. Per il gas di fermi perfetto valgono le leggi 

pV = n M RT R = k B N A N = n M Af A (5.65) 

Si definisce con uiq la massa di una particella , quindi 

N M 

M = ttiqN n = — p = — = rriQ-n (5.66) 

ove n indica la densita numerica e p quella di massa. La pressione si scrive come 

p = k B nT = kpT (5.67) 

mentre l'energia 

E=*n M RT^p=*e (5.68) 

dove e e la densita di energia per volume. Per il gas di Fermi non relativistico si ottengono le 
relazioni 



v ;3 3 h 2 /6tt 2 Y /0 5/3 2 h 2 /6tt 2 Y /0 5/3 , cfin . 



In questo caso la densita di massa e pari a p = m§n. Nel caso di gas perfetto relativistico 
bisogna tener conto della massa a riposo , scrivendo l'energia della particella come 

e(k) = ^(hck) 2 + (m c 2 ) 2 (5.70) 

Si ricava quindi l'espressione dell'energia per il caso relativistico 

E = (2^3 / d 3 ke(k)T(k) (5.71) 

dove 

m = \ 1 ' 0<k<kF (5-72) 
I altrove 

Si ottiene quindi 



E = ^ Jdk- k 2 ^{hck) 2 + (m c 2 ) 2 ^ e = ^ dk- k 2 m c 2 ^ (J^J +1 (5-73) 

Si definisce ora hk/moc = X c k = y e x = \ c kp dove A c e la lunghezza d'onda Compton. 
Sostituendo si arriva all'espressione 



" m f[ d,VVPTT-^^j[" <Hf-»VF+i (5.74) 



e = 



2tt 2 h? 

Svolgendo l'integrale , con semplici passaggi si arriva alia forma 



i 4 5 

8 2^^ X(X) ^ X(X) (5 - 75) 
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ove 



X 



(x) = x(l + x 2 ) 1/2 (l + 2x 2 ) -In x + (1 + x 2 ) 



2\l/2 



i 4 5 

1 rrigC 
82^2 ~p~ 



(5.76) 



Ricordandosi poi che x oc kp oc n 1//3 si vede che c'e una relazione esponenziale tra la densita 
di energia e la densita di particelle. La pressione si ottiene calcolando 



de 

p = n- e 

on 



P 



dE 
dV 



Utilizzando la (5.67) e la (5.73) si ottiene quindi 

1 v 



P 



3 (2tt) s 



d d k ■ hkvF{k) 



(5.77) 



(5.78) 



Conoscendo l'impulso e l'energia delle particelle la velocita v e univocamente determinata 
dalle relazioni 

{p e = mojv p P c hck 

e = m^c 2 



v 



e(k) 



(5.79) 



Utilizzando la (5.70) e introducendo la variabile y = X c k = hk/moc si ottienf 10 

7](j)(x) 



1 v mie 3 f x 



3 2vr h 3 



dy 



con 



x(l + x 



2U/2 



1 +ln 



x + (1 + ^)1/2 



(5.80) 



(5.81) 



Le equazioni (5.75) e (5.80) forniscono l'equazione di stato del sistema per il gas di fermi 
relativistico. Si puo calcolare anche la densita del gas : a differenza del caso non relativistico 
qualsiasi forma di massa si manifesta in energia, quindi 



1 

c z 



n 



x{x) 



(5.82) 



Risulta utile considerare i limiti di questo sistema. 

Limite ultrarelativistico Nel caso ultra relativistico rriQ — > o , equivalentemente, i>1 



hckp 3> rriQC . In questo limite la ( 5.70 ) si riduce a hck. Dunque e 



k F . Utilizzando la relazione per la densita 



n 



f^pk F si ottiene 



(4vr): 



La pressione si puo ricavare con le usuali formule termodinamiche 



Jd?k-hckT{k) 



(5.83) 



de 



1 /67r J 



2\ 1/3 



^ 71 dn £ 4 I v 



hen 4 / 3 



(5.84) 



10 



II fattore rj e equivalente a quello introdotto precedentemente, a meno di una costante moltiplicativa, 
ininfluente ai fini della trattazione. 
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Limite non relativistico II limite opposto e costituito dal caso non relativistico. Si ha 
quindi x <C 1 hckp <C m^c 2 . Per la densita di energia , cosf come per la pressione , 
si puo sviluppare in serie di Taylor rispetto ad x ottenendo 

8/3 3 \ 

x (x)~ -ix 3 + —x 5 - —x 7 + ...) (5.85) 

Arrestando lo sviluppo alia quinta potenza si ottiene 

8 5 2 h 2 /6vr 2 \ 2/3 5/3 
p = r}— x 5 = n 5/3 (5.87) 



15 5 2mo \ v J 

8 ( 3 3 
( x H 

3 V 10 



K' + TTT^ ) = £ + £i (5-88) 



dove si sono defmiti 



e = ^kpmoc 2 = moc 2 n ; e x ee r/^^x 5 = ^7^7- ( — ) n 5/3 (5.89) 



8_3_ 5 _ 3__7^_ /6^\ 2/3 
(i- 2 r " > '310^ ~52m V " / 

Si noti che £0 rappresenta il contributo della massa a riposo alia densita di energia mentre 
s\ il contributo dell'energia interna. 

5.11 Classificazione delle stelle 

1. M < 8M G : Nana Bianca 

2. 8M Q <M< 25M : Fe-core=>Supernova=>Stella di Neutroni 

3. M» 25M : Fe-core =>• Buco Nero 

Per il sole M = 1, 989 • 10 33 g , R Q = 6, 96 • 10 5 km p ~ 1, 4 g/cm 3 . 



Stella 


Massa 


R 


P 


R Scwh /R 


Sole 


M 


Rq 


~ 1, 4 g/cm 3 


4 


• 10" e 


Nana Bianca 


< 1,4M 


10- 2 R Q 


10 6 g/cm 3 


4 


• 10- 4 


Stella di Neutroni 


1^-2M 


~ 10km 


5 • 10 14 g/cm 3 




0,3 



5.12 Equazione di stato per core stellare 

Supponiamo innanzitutto che il materiale sia perfettamente ionizzato e che la temperatura di 



questo sia sufficientemente piccola da esser considerata nulk 1 ^ La pressione si puo scrivere 
come somma di due termini : P = Pi + P e ~ P e visto che Pj -C P e , dove / rappresenta gli ioni 
e e gli elettroni. Prima di ricavare l'equazione di stato occorre sviluppare alcuni argomenti. 



Ovviamente questa condizione deve essere rapportata con l'energia del livello di fermi : dovra essere 
kT < e F . 
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5.12.1 Limite per gas non relativistico sulla densita 

La densita totale del sistema si trova considerando la densita di tutti i nuclei e degli elettroni. 
In realta gli elettroni hanno una componente oc x 5 che contribuisce alia densita ma , in questo 
caso, questo termine e trascurabile. La densita si puo quindi calcolare come 



12 



p = (mA)nj ~ (A)mHni (5.90) 

Ovviamente n e = {Z)nj , quindi 

(A) 

P = J^ m H n e (5.91) 

Si puo definire Y e = cosicche p = 1 ^ L n e . La densita numerica degli elettroni e nota dal 
modello del gas fermionico(z^ = 2) : 

1,3 _ m H 1 /m e c\ 3 3 



Conoscendo la forma esplicita di p e possibile studiare i limiti i > 1 ex <C 1 con facilita. 
Conviene partire dal caso in cui x = 1 : si ponga per semplicita Y e = 0, 5. Si ha quindj 1 " 3 ] 

p * = p (x = l) = t™^) 3 - 1, 102 • l(T 6 MeV/c 2 • fm 3 ~ 2 • 10 6 g/cm 3 (5.93) 



3vr 2 0, 5 V h 

Quindi per p < p* gi elettroni non sono relativistici , viceversa sono ultra relativistici. Nel 
caso del sole p® = 1,3- 10 2 g/cm 3 quindi gli elettroni non sono relativistici. 

5.12.2 Condizione di gas ideale 

La condizione da porre perche il gas possa considerarsi ideale e che il termine di interazione 
nelPenergia sia piccolo rispetto all'energia di Fermi : |Vc ui| ep. Per ricavare questo limite 
esplicitamente si consideri innazitutto il volume medio degli ioni, calcolabile come 

V 1 4 o 
"< = N, = n, = 3" r » (5 - 94 > 

ove ro e il raggio medio del nucleo. La condizione di gas ideale si puo quindi scrivere come 

< - — k 2 F (5.95) 



47r£o ro 2m 

1 7 4 -\ ( 7 \ 1 ^ 3 

Si ha che — = — = g-vrrQ =4> ro ~ ( — J Quindi sostituendo nella diseguaglianza si ottiene 

Ze 2 i/q h 2 9 /o 9 / m e c 2 \ 3 Z 2 , 



La relazione finale per la densita si ricava calcolando 



p> -^-AZ ~ 11,3 -ZA g/cm 3 (5.97) 

Ad esempio per il Carbonio-12 il valore critico e di circa 8 • 10 2 g/cm 3 mentre per il Ferro-56 
la densita critica e di circa 1, 65 • 10 4 g/cm 3 



12 In questa approssimazione e stata trascurata l'energia di legame dei nucleoni presenti nel nucleo ( 6 -r- 7 
MeV). 

Per l'ultima parte dell'uguaglianza si utilizza il fattore di conversione 1 MeV/c 2 = 1,78266 ■ 10" 27 g ; 
lMeV/c 2 /fm 3 = l,78266-10 12 g/cm 3 . 
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5.12.3 Correzioni elettrostatiche all'equazione di stato per un gas di fermi 
perfetto 

Per il calcolo della correzione dovuta all'interazione elettrostatica si immagini di dividere il 
sistema in molte cellette (gli atomi del gas) : ognuna di queste risulta globalmente neutra. 
Questo tipo di schematizzazione e dovuta a Wigner-Seitz . La densita di carica e costante ed 
e data da 

Pe = -^ (5.98) 

Quindi l'energia di interazione elettrone-elettrone e ione-elettrone si calcola come 

r° q(r) , 3 Z 2 e 2 „ f r ° Ze , 3 Z 2 e 2 

E ee = ^dq=- ; E Ie = —dq = -- 5.99 

Jo r 5 r Jo r 5 r 

L'energia totale si ottiene sommando 

£coui = E ee + E Ie = = -^Z 2 e 2 y /3 Z-V3 n i/3 (5.100) 

Quindi l'eneregia media e data da 

^=-^V(f) V3 ^/3 n l/3 (5.101) 

Utilizzando le usuali relazioni termodinamiche e possibile trovare la pressione dovuta al ter- 
mine di interazione Coulombiano. 

2 d(E Coul /Z) 3 2 2 /4vr\ 1/3 2/3 1/3 1M > 
PCoul = n e = -— Z e — Z ' nj (5.102) 



dn e 10 V 3 / 

Sviluppiamo questo risultato nei limiti gia considerati di elettroni ultrarelativistici e nonrela- 
tivistici. 

Elettroni ultrarelativistici 

po = ^(3ir) 2 kcnj/ 3 (5.103) 



Quindi la pressione totale e data da 



P = P0 + PCoul = Po H = Po 



2 5/3 /3\ 1/3 0/3 
X -5 UJ 



(5.104) 



Ad esempio per C , Z = 6 , pcoui/f?o — 0, 015 ~ 1, 51% mentre per Z = 26 pcoui/.Po — 
4%. 

Elettroni nonrelativistici 

"^^---f 1 -^) (5 ' 105> 
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5.12.4 Equazioni Politropiche e di Lane-Emden 

Dopo queste considerazioni e possibile scrivere il sistema die definisce implicitamente l'e- 
quazione di stato per il core stellare 

f dp _ Q m(r)p{r) 

If „ r2 (5.106) 

\^=^r 2 p{r) 

Questo sistema si risolve fissando un valore della densita centrale p(r = 0) = p c e integrando 
fine- al raggio esterno. Le condizioni al contorno sono quindi date da 

= 

r=0 (5.107) 

p(R) = per un certo r = R 

5/3 

Elettroni nonrelativistici In questo caso p oc n e , quindi 

Y e 1 , 1 . 

n e = p = p dove p e = — (5.108) 

m H Pem H Z 

Con queste sostituzioni la pressione diventa 

P = l^ 2 ) 2,3 ^jhr 3 P" 3 - W 3 (5-109) 
° " l e p e rrijj 

Elettroni ultrarelativistici La pressione si scrive come 

?=i<^ 2 ) 1/s ^V 4/3 < 5 - 110 > 

* Pe 

In generale un'equazione di stato del tipo p = Kp 1 e deta equazione politropica , dove 7 = ^±1 
ed n e detto indice politropico. Si puo dimostrare le equazioni del sistema (5.106) possono 
essere combinate in una sola equazione. Infatti 

dp ^m(r)p(r) r 2 dp _ . . 

= —G J p — - — = -Gmr 
dr H p(r) dr ii 

_ G dm(r) _ d^ / ■r 2 _dp\ = _ G 4vrr , 2p(r) 
dr dr \p(r) dr J 



-AttGp (5.112) 



Si ottiene quindi 

1 d (r 2 dp\ 
r 2 dr \ p dr J 

f .w/ 1 — n 

Passando alle variabili adimensionali r = a£ , p = Pc$ n {€) , dove a 2 = p c n , si ottiene 14 

l'equazione di Lane-Emden 

Id / 2 d0\ 

Le soluzioni di questa equazione sono le funzioni di Lane-Emden $ n (£) : tali funzioni sono di 
vitale importanza in quanto il raggio stellare puo essere espresso in funzione di queste. 

14 Si osservi che, con questa sostituzione, p = p c i9 n+1 



i£ 2 ^) =-no (5.113) 
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5.13 Raggio stellare 



Rappresenta il valore di r per il quale la densita della stella si annulla : p(R) = = = 
Pc$ n (£l) = 0. II primo zero delle funzioni $ n rappresenta il punto in cui si annulla la densita 
, quindi definisce implicitamente la superficie del raggio stellare. Questo sara quindi dato da 
R = a£i dove £i e il primo zero delle t} n . Visto che a e determinato dai parametri del problema 
si puo ricavareUn metodo di risoluzione utilizzato per il risolvere numericamente Pequazione 
e quello di Runge-Kutta . 



R = a£i 



(n + l)k 
AttG 



1/2 l-n 1— n 

£ n 2n — r> „ 2n 



(5.114) 



Si noti che R oc p c . 



5.14 Massa stellare 



del sistema (5.106). 



La massa della stella si puo ottenere dall'equazione di stato , integrando la seconda equazione 

10 



f R 

M = m(R) = m(a£i) = Air dr ■ r 2 p(r) 

Jo 



(5.115) 



Cambiando variabile si ottiene 



M = Ana 3 p c / ?1 e$ n (0 de = ^a 3 p c ( -f^f) (5.116) 
Jo V d ? 



dove l'ultimo passaggio segue dalla (5.113). I valori delle costanti che compogono le equazioni 
di Lane-Emden sono facilmente calcolabili con metodi numerici: viene riportata una tabella 
con i valori delle costanti per alcuni ordinj 15 [ 



TABLE 4 



The Constants of the Lane-Emden Functions* 


■ 






»«/* 




A'„ 


H'„ 







2- 4494 

2.7528 

3- 14159 

3- 65375 

4- 35 28 7 

5 - 35528 
6.89685 

8 01894 

9 5358i 
14-97155 
31 .83646 

169,47 

OO 


4 . 8q88 
3-7871 
3 141 59 
2 . 7 1406 
2 .41105 
2 . 18720 
2 .01,824 
1 . 94980 
1 . 89056 
1 -79723 
1 73780 
1-7355 
1 . 73205 


I .OOOO 
I. 8361 

3 28987 
5 W071 
11.40254 
23.40646 
S4-182S 
88.153 
152 884 
622 408 
6189.47 
934800 

CD 


0.33333 
0.O2156 




. 1 19366 
0. 26227 
0.392699 
0.770140 
1. 63818 
3 . 90906 
1 1 . 05066 

20365 
40.9098 
247 558 
4922 . 125 

3.693X10 6 

03 


o-S 

0.53847 

o-S 

0.53849 
0.60180 
. 69956 
0.85432 
0.96769 
1 .13087 
1 . 66606 
3-33ioo 
16.550 

CO 


0-5 

I . o 


2 . 270 

O.63662 

O.42422 

O.36475 

0351SO 

0.36394 

O.37898 

40104 

O.47720 
O.65798 
I -340 

co 


«-5 

2.0 

2- 5 

3- o 

3 25 

3 5 

4° 

4 5 

4- 9 

5- 


132.3843 
10.4950 

3 . 82662 
2.01824 

1 S47I& 
1 . 20426 
0. 729202 
0.394356 
0.14239 




* The values for n — 0.5 and 4.0 are computed from Emden's inteerations of 0«; for n — 3.35 an unpublished integration by Chandrasekhar has been used, 
corresponds to the Schuster- Emden intejml. For the other values of n the British A nocialion Tables, Vol. II, has been used. 



Sostituendo il valore di a , gia ricavato in precedenza, nell'equazione (5.116) e utilizzando la 

2 d'd n 



(5.114) si ottiene infine 
M = 4vr 



n+l)K 
4vrG 



3/2 3^n 
2n 



Pc 



15 Una trattazione accurata viene data nella appendice [c| 



3 — n 
A n 2n 



A„. „3=Hl 
Rl-n 



D 



3 — 71 

l-n 



(5.117) 
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Si considerino ora alcuni casi particolari. Sen = 3/2 =^7 = 5/3: in questo caso la massa 
assume la dipendenza M 3 / 2 oc R~ 3 . Nel caso, invece , dira = 3=^7 = 4/3 1a massa diventa 
indipendente dal valore della densita centrale ( quindi dal raggio ): M^(R) =costante. Dato 
che si puo passare da un caso all'altro dobbiamo poter combinare i due andamenti : visto 
che p cresce al diminuire del raggio , andando a regimi ultrarelativistici, e ovvio pensare che 
Mz(R) sia valida per fl->0e viceversa M 3 / 2 (i?) per R grandi. Quindi la massa non crescera 
sempre ma avra un valore limite definito da 



M Ch = A 3 = 4tt 



4K 
4^G 



3/2 



-d^ 



G 



1 1 

2 9 



5,8 -Ma 



Me 



(5.118) 



Per il C il valore \i e = 2 , quindi Mq^ = 1,45 • M . La massa limite per la nana bianca 
appena calcolata rappresenta una schematizzazione molto semplice che non tiene conto delle 



interazioni elettrostatiche calcolate nel capitolo 5.12.3 Visto che pcoul e negativa le correzioni 
elettrostatiche diminuiscono il valore della massa limite : maggiore e Z e minore sara la massa 
limite a cui il rapporto M/Mq tende al variare della densita centrale. L'andamento descritto 
e ben visibile nel seguente grafico. 

modello di Chandrasekhar 




9.5 10.0 10.5 11.0 
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Capitolo 6 

Decadimento (5 



Negli anni successivi alia scoperta della radioattivita fu osservato un diverso comportamento 
delle particelle emesse dalle sostanze radioattive durante il decadimento. In molti casi gli 
strumenti di rivelazione mostravano la presenza di tracce simili a scie: quando era applicato 
un campo magnetico, le tracce provenienti da alcune sostanze radioattive avevano la curvatura 
rivolta verso un lato, altre avevano una curvatura rivolta nel verso opposto alle prime. Ai raggi 
associati alle tracce deviate da una parte fu convenzionalmente dato il nome di raggi alfa , a 
quelli con tracce deviate dalla parte opposta quello di raggi beta e ai rimanenti, il nome di raggi 
gamma .La scia emessa nel caso dei raggi beta e dovuta al fatto che viene emesso un elettrone 
nel corso del decadimento beta. II motivo per cui i tre tipi di raggi sono deviati in modo diverso 
dipende dalla diversa carica elettrica che hanno le particelle emesse: positive (se si tratta di 
particelle alfa) nel caso decadimento alfa, negative (elettroni) nel caso del decadimento beta 
e neutro (se si tratta di fotoni) nel caso del decadimento gamma. Inizialmente i decadimenti 
/3 risultarono decadimenti anomali visto che non veniva rispettata ne la conservazione del 
momento angolare ne la conservazione dell'energia : non era stato ancora scoperto il neutrino. 
Con l'introduzione di questa nuova particella da parte di Pauli si riuscf a rispondere a tutti 
gli interrogativi posti : il neutrino ha quindi q u = 0, m u = 0, s v = 1/2. I tipi di decadimento 
(3 sono dai da 

Decadimento /3~ : z^-N — > z+iX' iV _ 1 + e~ + V e 
Decadimento j3 + : z^-N — > z-iX' N+1 + e + + v e 

Cattura Elettronica(e) : e~+ z A '^N — > z-i^-'n+i + Infatti dopo la cattura elettronica di 
un elettrone il nucleo X' si trovera in uno stato eccitato e quindi decadra , emettendo 
un'energia pari alia binding energy dell'elettrone. 



Per scrivere il Q-rate di questa reazione si analizzi innanzitutto il decadimento del neutrone 
libero 



Ovviamente il Q si puo scrivere anche come somma delle energie cinetiche dei prodotti : 
Q = T p + T e + T u . Stimando questa quantita si osserva che T p ~ 0, 3 KeV quindi e trascurabile 
nella somma, pertanto Q ~ T e + T v = T emax = Qq Sperimentalmente si misura Q[j xp = 



6.1 Q-rate 



n — > p + e + v e con Q = (m n — m p — m e — m v )c? 
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(0, 782 ± 0, 013) MeV ; utilizzando questo valore si puo ricavare un limite superiore alia massa 
a riposo del neutrino 

Q = c 2 (m n — m p — m e ) — m u c 2 = Q[j xp — m u c 2 

Si assuma per il momento m u = visto che risulta molto piccola. I Q di reazione si possono 
facilmente calcolare nei tre casi 

Decadimento /3~ Si ha 

Qp- = M(Z,A) -M(Z + 1,A) -m e (6.2) 

dove 

(m(Z,A)=M(Z,A)+Zm e -B e (Z,A) 

| m(Z + l,A) = M(Z + 1, A) + (Z + l)m e - B e {Z + 1,A) 



Facendo la differenza di queste ultime due equazioni e sostituendo in (6.2) si ottiene 

Qp- =m(Z,A) -m(Z + 1,A) (6.4) 

infatti la differenza delle due binding energy puo essere trascuratsQ visto che e nell'ordine 
dell'energia di legame dell'elettrone esterno (~ 10 eV). 

Decadimento /3 + 

Qp+=m(Z,A)-m(Z-l,A)-2m e (6.5) 

Si noti la differenza tra i due decadimenti /3~ e /3 + : il secondo puo avvenire SO- 
LO quando la differenza delle due masse e superiore a due volte la massa a riposo 
dell'elettrone. 

Cattura Elettronica 

Q £ = m(Z,A)-m{Z-l,A)-b e (6.6) 

6.2 Teoria di Fermi del decadimento (5 

Si supponga che il sistema nel caso del decadimento /3 sia descrivibile tramite due Hamilto- 
niane : l'Hamiltoniana iniziale Hq e l'Hamiltoniana finale, dopo il decadimento , H' . Uti- 
lizzando la teoria delle perturbazioni si puo calcolare la probabilita di passare dallo stato 
deH'Hamiltoniana iniziale a quello dell'Hamiltoniana finale. La probabilita e data da 

9-7T 

W=^ f \H'\i; i )\ 2 p(E f ) (6.7) 

questa analisi va sotto il nome di Regola d'oro di Fermi. La funzione ipi e associata al nucleo 
iniziale, e autofunzione, e sara identificata dalla terna \rpi) = \ Jj, mi, Lo stato finale ipf , 
senza considerare l'interazione tra le particelle, si puo scrivere come 

= |*/)|0e)|^> (6-8) 

dove le funzioni d'onda rappresentano , rispettivamente, lo stato del nucleo , dell'elettrone e 
del neutrino. La p(Ef) rappresenta la densita di stati , ovvero il numero degli stati compresi 
tra l'energia Ef e Ef + dEf, ovvero 

i i dn 

1 I Q dei decadimenti /3 sono , generalmente, nell'ordine del MeV. 



6.2. TEORIA DI FERMI DEL DECADIMENTO (3 



71 



. Per la funzione d'onda del neutrino si pud prendere un'onda piana, normalizzando su un 
volume V. 

< r -|^)= * e 1 ^ (6.9) 



Per l'elettrone si supponga che questo sia, in prima approssimazione, libero in modo che la 
funzione d'onda si riduce nuovamente a 

<^ e ) = JLe 1 ^ (6.10) 
v V 

Gli elementi di matrice si calcolano come 

^f\H'\^) = 1 J e- ^ e" ifcV ^}(n, . . . , f A )H'A(fi, r A ) d 3 r e d 3 r, d 3 n . . . d 3 r A 

(6-11) 

Fermi ipotizzo che l'interazione fosse a raggio zero e che il neutrone si trasformasse, dentro al 
nucleo , in p + e~ +v e . Questa inter azione si puo scrivere come 



H' = gO'J 3 {r e -r i )5 3 {f 



r; 











. . . , r A ) d 3 r e d 3 rv d 3 ri . . . d 3 ?"^ 




ikv 




.,f A )0'j 3 (f e 


- n)5 3 (r u - fi)tpi(fi, ■ ■ ■ , r A ) d 3 r e d 3 ?v d 3 r x 




ik u 




. ,r A )0' x if)i{fi, 


. . . , r A ) d 3 ri . . . d 3 r A 



con g una costante che tiene conto delle dimensioni.L'operatore 0' x e l'operatore che trasfor- 
ma neutrone in protone, quindi sara proporzionale all'operatore di salita dell'isospin. Le 5 
invece tengono conto che l'interazione e a cortissimo range. Con queste ipotesi l'integrale e 
ulteriormente semplificabile : 

V 
_ 9_ 
V 

(6.12) 

Per semplificare ulteriormente si cerca di stimare i valori tipici di impulso per le particelle. Si 
supponga che T e = 1 MeV; si puo stimare il k e dalle relazioni 

E e = sjipec) 2 + (m e c 2 ) 2 ^E e = T e +m e c 2 =>- Pe c -1,4 MeV k e = jr = ~ 0, 7-10" 2 fin" 1 

(6-13) 

Si noti che , nell'argomento dei due esponenziali della funzione integranda, ri — > oo nell'in- 
tegrazione : in realta pero la identicamente nulla fuori dal nucleo, ovvero per valori di 
rj > R ( R e il raggio nucleare), quindi l'integrale e non nullo solo per i valori 

\k e - n\ < k e n < k e R ~ 7 • 10~ 3 fm • R ~ 0, 01 4- 0, 03 < 1 

Quindi i due esponenziali si possono approssimare come 

e-^slse-^ (6.14) 

L'approssimazione (6.14) e detta delle Transizioni permesse^J Utilizzando la (6.14) nell'inte- 
grale (6.12) si ottiene infine 

" V 



{^ f \H'\^ i ) = ^ I y}(f 1 ,...,r A )0' x iP l (f 1 ,...,r A )d 3 r 1 ...d 3 r A (6.15) 



2 Si osservi che tale approssimazione e giustificabile anche in virtu dello sviluppo in serie. Infatti 

oo 
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Questa quantita rappresenta semplicemente gli elementi di matrice dell'operatore 0' x quindi 
si puo scrivere come 

WH'm = ^(* f \0'M) = f M Pi (6.16) 
Conviene ora concentrarsi sulla densita degli stati p(Ef) , avendo gia semplificato il termine 



matriciale della (6.7) . 



6.2.1 Spettro energetico degli elettroni 

Si indichi con dn e la densita numerica numerica di elettroni : e gia stata ricavata la formula 



V 2 V o . 

dn e = ^ 2 ^^ 3 47rfc e dk e = j^^Pe d Pe ( 6 -17) 



Analogamente per i neutrini 



V 2 

dn u = -j-^47rp u dp u (6.18) 



Quindi il numero di elettroni e di neutrini si puo ottenere come 

V 2 



dn e dn v = ^(infpl dp v p 2 e dp e (6.19) 



La probability di transizione tra i due stati si puo quindi scrivere come 



dW= 2 ^ f \H'm 2 ^^ (6.20) 
Utilizzando l'approssimazione delle transizioni permesse 

L'energia finale e data da Ef = E e + p u c , avendo considerato m u = 0. Se si immaginia di 
fissare un valore di E e , si ricava 

dp u _ 1 
dE f ~ c 



quindi si puo sostituire nella (6.21), ottenendo 



dW = -^ r g 2 M Yl p 2 e dp e pl = Cp 2 e dp e p 2 u (6.22) 

Dunque la funzione di distribuzione per l'impulso si puo scrivere utilizzando la forma di W e 
le relazioni 

\ N / P f P '; C f: dP < *»c = Q-T.* r . = \lQ-Tj (6.23) 
[Q = T e + T v = T e +p u c c 

Quindi 

N e {p e )= C ^p 2 e {Q-T e ) 2 (6.24) 



che , per £ = si riduce a 



i h -r ■ f i \ Sill hV 

J"-= 3o (kr) = — ->1 
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L'energia cinetica per elettrone e della forma T e = \J (p e c) 2 + (m e c 2 ) 2 — m e c 2 , che si puo 



sostituire nella (6.24) per avere la forma completa della distribuzione. 



N e (p e ) = -^pl [Q - y 7 {Pec) 2 + (m e c 2 ) 2 - m e c 2 J 



(6.25) 



La figura seguente mostra l'andamento dello spettro per la (6.25) e per lo spettro dell'energia 



cinetica, ottenuto dalla (6.25) con un cambio di variabili 



N e (T e ) = ^ (T 2 + 2m e c 2 T e ) 1/2 (Q - T e f{T e + m e c 2 ) 



(6.26) 




z.c 

7, (MeV 



Sperimentalmente si osservano alcuni discostamenti da quanto visto teoricamente. Ad esempio 
per il nucleo 29CU35 esistono i due decadimenti 

29Cu 35 ^ ||Zn + e -+X7 e (6.27) 
29CU35 -> 1 Ni + e+ + v e (6.28) 



Gli spettri dei due decadimenti sono shiftati, come si vede dalPimmagine seguente. 
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■ L 2 C 

7, IMtV 



Questo effetto e dovuto alle interazioni elettrostatiche che modificano l'Hamiltoniana : in 
particolare la funzione d'onda non e piu approssimabile con un'onda piana. Per risolvere 
questo problema Fermi introdusse un fattore adimensionale che tenesse conto delle interazioni 
elettrostatiche : F(Z',p e ) e detta funzione di Fermi e si puo esprimere in forma approssimata 
come 

„ x 2TraZ' 

F = — dove x = ± — (6.29) 



dove il segno + e per il decadimento /3 e viceversa. Aggiungendo questa correzione la densita 



della (6.24) diventa 



N e {p e 



C 



p 2 e (Q-T e ) 2 F(Z', Pe ) 



(6.30) 



Per avere un riscontro con i dati sperimentali conviene prima trasformare questa formula , 
linearizzando l'andamento ( plot di Kurie ). 



NM - VC (Q-T e ) (6.31) 



plF(Z',p e ) C 



Sperimentalmente si verifica che l'andamento previsto e rispettato , a meno di scattering a 
basse energie. 
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6.2.2 Vita media 

Si vuole ora trovare una relazione per la vita media. 



dW = l2 J^g 2 \M Fi \ 2 F(Z', Pe )p 2 e (Q - T e f d Pe (6.32) 
Integrando su tutti gli impulsi permessi si ha 

W = j dW = ^$^r F (Z',Pe)pl(Q-T e ) 2 dp e (6.33) 
Per risolvere l'integrale conviene prima defmire una quantita adimensionale data dg[^] 

- {m J M -> ST f{z '^ {q - T * )2dp < <6 - 34> 

Definendo Eq = S™ ax si ottiene l'espressione finale per la probability di transizione 

W = ^ M ^ Z '^ = \ = ^ 2 < 6 - 35 > 

Si definisce una vita media comparativa data da 

2ir 3 h 7 1 

T(Z', Eo)t 1/2 = In 2 5 4 (6.36) 

Dopo aver dato questa definizione si possono dividere i decadimenti in alcune classi : in parti- 
colare le transizioni con log J-t = 3 -j- 4 s sono dette super-permesse. Dai risultati sperimentali 
e possibile ricavare g ~ 0, 88 • 10~ 4 MeV- fm 3 e quindi G = ^rg — 1CP 5 nel caso di inter- 
azione nucleone-nucleon , molto piccolo se confronto , ad esempio, con g 2 dell'interazione forte 
pione-nucleone (~ 1) o anche con l'interazione elettromagnetica tra due protoni (~ 1CP 2 ). 

3 La quantita ( m e )3( m c '2y2 & stata aggiunta per mantenere Padimensionalita . 
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6.2.3 Limite per massa nulla del neutrino 

Se si toglie l'ipotesi m u = si puo fare una stima dell'upper limit per la massa del neutrino. E 
gia stato visto che 

dW ~ pi dpePv ^' 



dE f 

Per buona parte dello spettro , se siamo nella condizione m v c 2 <C Q, si ha che l'energia del 
neutrino e molto maggiore di m u c 2 :in queste zone quindi e una buona approssimazione quella 
ultra relativistica, per la quale E u ~ p u c. Nella zona vicina all' end point , cioe dove £ p ->0e 
quindi E u ~ m v c 2 , l'energia cinetica del neutrino diventa piccola.In questo limite (p e — > p™ ax ) 
il neutrino si puo quindi considerare non relativistico : l'energia si puo quindi approssimare 
come 



E v 



Vv . 2 

h m v c 



L'energia totale e data da Ef = E v + E e ^> ^ = ^ => dW ~ pIpI^ &Pe- Si ottiene poi 



Pv 



Jhn~ v {Q - T e fl 2 



(6.37) 



La differenza sostanziale tra il caso di massa non trascurabile e il suo opposto si nota se si 
computano i limiti 



m u = 



lim 



dN e { Pe 
dpe 



m u / 



dN e (p e ) 
urn — ; = — oo 



dp e 



Sfruttando questa proprieta si puo ricavare sperimentalmente un upper limit per la massa del 
neutrino. La situazione e ben riportata graficamente nella figura seguente 




2966 32 keV/r 



T967 32 WV/< 



2500 «eV 
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6.3 Regole di selezione 

E gia stata analizzata la regola di selezione per le transizioni permesse : si richiedeva che la 
funzione d'onda fosse approssimabile con 

e (r) = JLe^~l (1 = 0) 



Transizioni Permesse Avendo £ = , dalla conservazione del momento angolare totale , si 
ha che Ii = If + S , dove S = s e + s v . A seconda del valore di S le transizioni permesse 
vengono divise in varie classi. 

• S = 0. Transizioni di Fermi : Per queste Ii = If , Ui = Uf. Sono transizioni 
del tipo 0+^0+. 

• S = 1. Transizioni di Gamow- Teller : Per queste L = If + 1, . . . , \If — 1| =4> 

AI = \If — Ii\ = 0, 1 e ITj = Uf. In questo caso la transizione + — > + non puo es- 
serci perche si dovrebbe accoppiare un momento angolare con spin 1. Appartiene 
a questo tipo la reazione He— > 6 Li+e~ + v e che e di tipo + — > 1 + 

Transizioni Proibite Questi decadimenti non sono proibiti , ma sono meno probabili delle 
transizioni permesse. Si incontrano frequentemente quando lo stato iniziale e finale 
hanno parita opposte , e quindi quando e violata la regola per le transizioni permesse. 
Per soddisfare il cambio di parita , Pelettrone ed il neutrino dovranno essere emessi con 
un valore di £ (momento angolare rispetto al nucleo) dispari. Ad esempio si consideri il 
momento angolare che assume un elettrone nello scattering con un nucleone. Si supponga 
che il nucleone abbia raggio R; l'elettrone incide con impulso p e = 2 MeV/c quindi riesce 
ad assumere un momento angolare pari a £ = p e R = 10 MeV-fm/c , ovvero , in unita di h, 
£/h = 10/200 = 0, 05 <C 1 = k e R. Dunque le transizioni proibite diventano importanti 
solo quando il nucleo non puo fare una transizione permessa, altrimenti il loro contributo 
e insignificante. 

• £ = 1. Transizioni Prime Proibite : Per queste ovviamente Ii = If + S + £ 
quindi AI = 0, l,2ellj = — Hi. Si distinguono in base al valore di S. 

— S = . Transizioni prime proibite di Fermi : per queste AI = 0, 1 e 
U f = -Ui. 

— S = 1 . Transizioni prime proibite di Gamow- Teller : per quest^] 
AI = 0, 1,2 e U f = -Ui. 

Alcuni esempi di queste transizioni sono 

17 N^ 17 (f ^f) (6-38) 
76 Br 76 Se (I" -> 0+) (6.39) 

• £ = 2. Transizioni Seconde Proibite : Per queste AI = 2,3 e ITj = ITj. Alcuni 
esempi di queste transizioni sono 

22 Na -> 22 Ne (3+ -»• 0+) (6.40) 
137 Cs^ 137 Ba (I + (6.41) 

4 Si noti che le transizioni + — > CP sono necessariamente di Gamow- Teller. Infatti se I = If = l'unico 
modo per avere I + s — con I = 1 e S = 1 . 
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Capitolo 7 

Modello vibrazionale del nucleo 



7.1 Introduzione 

II problema della collocazione dei nucleoni nel nucleo non e di facile risoluzione. Evidenze 
sperimentali mostrano che si puo assumere per i nucleoni una distribuzione secondo il modello 
a shell , come gia fatto per gli elettroni esterni alio stesso nucleo. Questo comportamento e 
giustificato teoricamente dalle seguenti osservazioni 

1. A distanze piccole l'interazione nucleone-nucleone e fortemente repulsiva, dunque i nu- 
cleoni non possono avvicinarsi molto. 

2. I nucleoni sono fermioni quindi devono riempire i livelli energetici nello stesso modo degli 
elettroni. 

La combinazione di queste due proprieta fa si che i nucleoni si distribuiscano secondo un 
modello a shell. Supponiamo invece che il nucleo venga eccitato : inizia un moto di vibrazione. 
Per nuclei non troppo pesanti (A < 150) il modello basato su vibrazioni attorno alia forma 
sferica di equilibrio funziona abbastanza bene; per nuclei pesanti la forma non e phi sferica ma 
ellissoidale. Per questi ultimi nuclei si avranno quindi altri modi di eccitazione, in particolare 
moti rotazionali attorno ad un asse NON di simmetria. Si osservi che i nucleoni non possono 
saltare facilmente da una shell all'altra dell'interazione forte. Se forniamo meno 

energia di quella necessaria al salto succede che i nucleoni, leggermente spostati dalla posizione 
di equilibrio , tentano di tornare nelle loro configurazioni fondamentali, dando luogo a moti 
collettivi, cioe modi vibrazionali attorno alia forma di equilibrio.il primo stato eccitato del 
nucleo sara quindi uno stato collettivo. Si puo cercare il primo atomo che ammette uno stato 
eccitato 

• A = 2 : L'unico sistema con stati legati e dato dal Deuterio. L'energia di legame e di 
circa -2,2 MeV. Questo sistema NON ha stati eccitati. 

• ^4 = 3: Esistono due elementi , 3 H e 3 He . Le energie di legame sono rispettivamente 
di -8,5 MeV e -7,7 MeV. Si noti che il trizio e radioattivo e decade naturalmente in 3 He. 
Entrambi questi sistemi non hanno , pero , stati eccitati. 

• A = 4 : E presente la particella a , data da 4 He. L'energia di legame e molto alta, 
circa 28,3 MeV : per rompere l'atomo e trasformarlo in un Deuterio phi un protone 
all'infinito si devono dare 20 MeV di energia. Questo suggerisce che il sistema non abbia 
stati eccitati : l'unico stato eccitato possibile si trova sopra alia soglia quindi si tratta 
di una risonanza. 
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• ^4 = 5: Non esiste nessun sistema di 5 nucleoni che sia legato. Infatti se si provasse 
ad avvicinare un nucleone alia particella a, cercando di formare 5 He la forza nucleare, a 
corto range, non sarebbe sufficiente per vincere il principio di Pauli che non permette il 
posizionamento di un solo nucleone in un nuovo livello della shell. 

• ^4 = 6: Esistono gli elementi 6 Li , 6 He. II primo e phi leggero del secondo , mentre 
quest'ultimo decade j3~ nel primo. Questi elementi NON hanno stati eccitatfj] 

• A = 7 : II 7 Li rappresenta il primo nucleo con uno stato eccitato. 

• ^4 = 8: Non esiste un nucleo con 8 nucleoni. 

7.2 Oscillazioni Nucleari 

Si consideri il nucleo come una goccia sferica di liquido nucleare. Nello stato eccitato questa 
goccia iniziera a vibrare. La perturbazione puo essere descritta considerando le oscillazioni 
del raggio intorno alia sua posizione di equilibrio; in formule si ha 



Con questa parametrizzazione e con diversi valori di [i e A si possono descrivere varie defor- 
mazioni della goccia. 



• A = : Si ha una deformazione lineare poiche R($,(p) = R(l + costante) , quindi la 
sfera si contrae o si gonna. 

• A = 1 : Si ha un moto di oscillazione del centro di massa , quindi la circonferenza si 
muove senza deformarsi. 

• A = 2 : Si ha un'oscillazione quadrupolare. 

• A = 3 : Si ha un'oscillazione ottupolare. 

Si noti che per A < 6 il tentativo di eccitare il nucleo porta ad una demolizione dello stesso. 






7.2. OSCILLAZIONI NUCLEARI 
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X = 1 
(Dipole) 





X = 2 
(Quadrupo(e) 



X = 3 
(Octupole) 



Visto che si vogliono descrivere stati di bassa energia si deve imporre che il volume si mantenga 
costante nella deformazione, quindi si deve imporre l'uguaglianza 

4-7T 



-Rn 



r 2 dr dQ 



(7.2) 



Si ottiene quindi 



r 2 dr dJ7 



r>3 

r'drdQ = I dtt— (<?,¥>) 



3 



dO 



XfJ 



(7.3) 



i + 3 ^2 a \ti Y \,t* + 3 1^2 



A/i 



A/t 



Infatti, visto che si tratta di piccole oscillazioni, il cubo puo essere sviluppato fino al secondo 
ordine. Svolgendo l'integrale si ha 



dft 



1 + 3 ^2 a A M y A, M + 3 1^2 a x f j, Y x, /j 



A/t 



A// 



3 
3 

47T 



4vr + 3^a A/ ,y dOr AjM (tf, ^Yq.o v 7 ^ + 3 ^ |a AM | 2 / dQ\Y x J 2 ($, tp) 

Xfj, " A/i 



'Air 



4-7T + 3 • — — a o + 3 ^2 \ a ^(' 



4tt 
3aoo 



Xfi 



(7.4) 



A/f 



Si noti che tra la seconda e la terza uguaglianza e stata usata la condizione J dQY£ „Yg m 
SxeSfim ottenendo 



Xfi Xfi 



Xfi 



(7.5) 



47T V47T 
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: l'unico termine che contribuisce alia somma , infatti , e il termine zeresimo. Non resta 
,quindi, che porre Puguaglianza 

^o 3 = f ( 4 - + 3 ^|+ 3 EM 2 ) (7-6) 

Visto che A = tiene conto delle compressioni ed espansioni del volume si puo eliminare il 
termine nella somma con aoo : e stato supposto che il volume venga mantenuto costante. 
Inoltre non si conta il termine con A = 1 in quando tiene conto di semplici traslazioni. Si 
ottiene quindi 



oo A 



A=2/u=-A 

avendo riconosciuto lo sviluppo (1 + x) 3 ~ (1 + 3x). Si ricava quindi 
R R ° 

7.3 Modello dinamico 



A=2,u=-A 



(7.7) 



R(#, ip) = R a ( 1 + E a ^A^ <P) ) ( 7 - 8 ) 

V A=2 ii J 



Si deve ora costruire un modello dinamico per la goccia di liquido nucleare. II problema 
consiste nel trovare l'Hamiltoniana del sistema in modo da calcolare gli energetici da evolvere 
nel tempo. Dato che si tratta di piccole oscillazioni l'energia cinetica e data da 



2 

A/i 

Per giustificare questa formula occorre costruire il sistema che descrive il fluido nucleare come 
incompressibile e irrotazionale. Vale quindi il seguente sistema 

+ (7.9) 

La velocita puo quindi essere scritta come gradiente di una certa funzione scalare (p : v = V(f> 
, dove (f> rappresenta il campo scalare delle velocita. Supponendo p costante (ipotesi di fluido 
incomprimibile) si ottiene 

fclpW -P(W) = ,V(V0)=pV^=-^ = O (7.10) 



Quindi si ottiene l'equazione di Laplace 



V 2 = 



La soluzione dell'equazione di Laplace sono date dalle potenze oc r A ; in particolare 

^ = E 6 V^a,m(^) (7-11) 
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Le soluzioni irregolari oc r~( A+1 ) sono state scartate a causa della loro irregolarita nell'origine 
: dato che il punto appartiene al fluido occorre, infatti, che il campo di velocita sia ivi ben 
definite Le condizioni al contorno da porre sono 



v r = g = R(#,<p) 
in r = Rq 



(7.12) 



Sviluppando la (7.12) si ottiene 



d<f> 
dr 



r=R o Xfi 



dR 

~dt 



= Ro E ®^ Y ^ 



Confrontando i coefficienti delle somme si ha 



^b\^Ro X 1 — Roa-xn =>• bx^ — a Xfl — 



-X+2 



(7.13) 



(7.14) 



Quindi la soluzione particolare e data da 



R 



-A+2 



X/x 



Conoscendo <j) si puo facilmente calcolare l'energia cinetica, infattQ 



T= 1 p [ v 2 {r)d 3 r = l pf |Vc/>| 2 d 3 r =-pf 4>*V(f>dS 
2 Jv 2 J 2/5 



Passiamo quindi al calcolo esplicitcP] 



(7.15) 



(7.16) 



R~ x ' +2 

E ~^T> a AV rA Y X'v> 



X'fi' 



Hi, 



£ 

X'fi' 



R -X'+2+X 



V+2+A' 

•* v* 

a X'fi' Y \'fi' 



w=E 



R 



-A+2 



-ax^\r x - l Y Xll 



Y^Ro X+2+X - l ax»Y M 



Xfj, 



Ro 



Xfi 



(7.17) 
(7.18) 



Integriamo sulla superficie infinitesima ( a raggio Ro fissato ) : dS = Rq dQh , dove h e il 
versore normale alia superficie. 

I ^V^dS = E E TT^'m'^m^ / *W ^ = E E T7 A AvMaV = E T 1 ^ 



Riprendendo la (7.16) si ha quindi 



2 r ^ A 

Xfi 



(7.19) 
(7.20) 



L'ultima uguaglianza segue dal teorema della divergenza. 
3 Utilizziamo due indici diversi A,/i e A',// perche nel prodotto non si potra subito confondere i due termini 
in un'unica somma. 
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Per l'energia potenziale bisogna considerare due termini : uno Coulombiano e uno di 



superficie : V = Es — E§' + Ec — E^ con 



E s = a § dS = a J J 1 + | VR\ 2 R 2 dQ 



II potenziale ha dunque la forma 



'4 0) = <?^Rl = aAnr 2 A 2 / 3 = a s A 2 / 3 



doveC^ = Cf-C c x 



A/Lf 



Sono stati deflniti inoltre per comodita le costanti(al variare di A) seguenti 

Cf = aRl{\-\){\ + 2) 



C? 



3 (eZ) 2 A - 1 



2vr R 2A + 1 
L'Hamiltoniana finale assume quindi la formaQ 

•Hcoll = n {B\\a\n\ 2 + Ca^A/J 2 ) 



(7.21) 
(7.22) 

(7.23) 
(7.24) 

(7.25) 



Questo rappresenta una somma finita di oscillatori armonici indipendenti. Calcolando il 
momento coniugato tt\^ alia variabile oa^ si ha 



dT 
da X u 



B x a^ =► H coll = \Y. (ir'^l 2 + C x \ax^ (7.26) 

Xfj, 



Con queste definizioni risulta immediata la regola di commutazione [a\^,ii\^\ = ih. Si noti 
che in H co \\ tutti i modi sono disaccoppiati ( modi normali ), in virtu della (7.26). Utilizzando 
il formalismo di seconda quantizzazione si possono introdurre gli operatori di creazione e 
distruzione PXfj,,Pxa > defmiti da 



L'Hamiltoniana si trasforma quindi in 

HcoU = E h^wx^xti + V 2 ) 
X/j, 

I livelli energetici si calcolano quindi banalmente con 



con lu\ 




1 7.27: 



(7.28) 



hu\ ( N + 



2A + 1 



(7.29) 



Infatti fj, rappresenta il numero magnetico, quindi assumera (2A+1) valori. Con questo modello 
si riesce quindi ad avere una stima accurata dei livelli energetici. Nello stato fondamentale 



II pedice coll sta per collettivo. 



7.4. SPETTRO VIBRA ZIONALE E SPETTRO ROTAZIONALE 
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qualsiasi nucleo descritto da questo modello si trova nello stato J n = + . Gli stati successivi 
si possono costruire applicando P^ At |0 + ) . A seconda dell'oscillazione (A = 0, 1, 2 . . .) variera 

quindi il primo stato eccitato : nel caso ottupolare, ad esempio, p\^\0 + ) = 3 _ . Dal modello 
collettivo si possono ricavare altre importanti proprieta. A tal proposito si supponga di voler 
costruire lo stato di piu fononi dallo stato fondamentale |0) con momento angolare definito : 
si dovra applicare p\^p\^\0) ■ Dato che i fononi sono bosoni ci sara simmetria per scambio di 
stato, ovvero 

£ (A^AH^Mbj^JO) = \J,M) = Y, ^2-Api\J,M)p{^p{^) (7.30) 

Utilizzando la proprieta dei coemcienti di Clebsch-Gordan si ha quindi 

(\Hi;Xfi 2 \J,M) = (-l) J (A/x 2 ;A/ii|J,M) (7.31) 

Quindi, per quanto visto prima, J deve essere pari. Dunque gli stati eccitati di due fononi 
soncj^ 4 " , 4 + . La situazione descritta e riportata nella figura seguente 

Hexadera|»le (2 A ) 
4 H 

<C = 10, 2, 3 

0, 2, li, 4, (y ^tZjt F 4 - °* OctujH^c (23) 

3-phonon 



2-phonon 

2+ 

1-phorran 



Pure Vibrational Hipher-order 

vibrational states in vibrntinnH 

re;il nuclei 

00 (b) (c) 



7.4 Spettro Vibrazionale e Spettro rotazionale 

Per verificare la bonta del modello elaborato si puo osservare se sopra al secondo stato eccitato 
dell'elemento e presente una nuova banda di stati. Altrimenti si puo misurare una transizione 
elettromagnetica : il valore di questa sara proporzionale al numero di nucleoni che fanno la 
transizione. In particolare un indice di validita del modello sara dato dal rapporto tra l'energia 
del secondo stato eccitato e quella del primo stato eccitato, al variare di A. II valore di questo 
rapporto e riportato nel grafico seguente. 

In realta intorno alio stato 4 + ci sono altri stato + e 2 + ad energie diverse anche se molto rawicinati : 
questa anomalia rappresenta una delle deviazioni dal modello. 
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Visto che il rapporto E{A + ) / E{2 + ) rappresenta un indice della validita del modello, da questo 
grafico si possono trarre importanti informazioni. Si osserva innanzitutto che , a fronte di 
un valore aspettato di ~ 2 , il rapporto varia tra 1 e circa 3,5 : il modello e applicabile solo 
in un certo range di nuclei. In particolare nella zona in cui A ~ 90 -f- 120 si ha un rapporto 
costante di 3,33 : questo significa che, dopo lo stato 2 + compare uno stato 4 + . Questo si 
verifica perch aggiungendo neutroni il nucleo non ha piu forma sferica ma tende a schiacciarsi 
diventando un'ellissoide. Questo tipo di forma geometrica per i nuclei comporta la comparsa di 
uno spettro rotazionale. Gli autovalori per questo spettro si trovano facilmente considerando 
l'equazione di Schroedinger. L'Hamiltoniana del sistema rotante e data da 

#rot = ~ (7.32) 
dove / e il momento di inerzia del sistema. Risolvendo si trova 

~V = EjiP R 2 Y JM = J {J + l)h 2 Y JM ^Ej= J i J -plh 2 con J = 0, 2, 4 . . . (7.33) 

Quindi E\ ~ A(A + 1) e dunque il rapporto tra due stati sara del tipo 

E n n(n + 1) 
E m m(m + 1) 



Ovviamente gli spettri rotazionali e vibrazionali possono combinarsi tra di loro dando luogo 
a spettri roto- vibrazionali. 



Capitolo 8 

Fissione Nucleare 



Nella fissione il nucleo di un elemento pesante - ad esempio uranio-235 o plutonio 239 - decade 
in frammenti di minori dimensioni, ovvero in nuclei di atomi a numero atomico inferiore, con 
emissione di una grande quantit di energia e radioattivit. A differenza della fusione nucleare 
questo processo e conveniente solo per elementi molto pesanti, come mostra il grafico della 
binding energy. 




" *0 « B 30 60 90 120 ISO 180 210 210 

MASS NUMBERS A 

Si possono identificare due tipi di fissione 

Fissione Spontanea Avviene spontaneamente in un nucleo senza eccitazione. 

a Xn ^ M y Ni + M w N2 +k-n (8.1) 

dove k e un intero. Ovviamente dovra valere che A\ + A2 + k = A e Z\ + Z2 = Z. I 
neutroni emessi spontaneamente vengono chiamati neutroni pronti e vengono emessi in 
un tempo di circa 1CT 16 s. 

Fissione indotta II nucleo puo essere eccitato con dei neutroni in modo che raggiunga lo 
stato X* , detto nucleo composto. 

^X N + n^ A+l X* -^Y + W + k-n (8.2) 



< 
LU 



87 



ss 



CAPITOLO 8. FISSIONE NUCLEARE 



Se si definisce con v = (k) si puo dare una stima dei neutroni liberati : ad esempio nell' 235 U 
vengono emessi v = 2, 423 + 0, 066 • E n , con E n energia dei neutroni ( E n < 1 Mev ). L'emis- 
sione di questo numero di neutroni comporterebbe una reazione incontrollata che porterebbe 
ad una sicura esplosione : per reazioni controllate occorre quindi progettare un sistema di 
assorbimento dei neutroni. Osservando il grafico della Binding Energy si nota che tra il nucleo 
limite per la fissione , 56 Fe e i nuclei piii pesanti con A = 240 c'e circa 1 Mev di differenza 
: questo significa che nella fissione viene liberato circa 1 Mev per nucleone, ovvero Q ~ 200 
MeV. Queste energia viene liberata sotto forma di energia cinetica dei frammenti. 



8.1 Reazioni favorite 

Utilizzando il modello a goccia e possibile avere una stima per prevedere quali sono i valori di 
A per i quali una reazione di fissione risulta favorita. Si consideri , per semplicita , il caso di 
una reazione simmetrica 

iX N ^2 A z f 2 Y N/2 ,Q = M(A,Z)-2M(^,j) (8.3) 

Visto che , in generale, M(A,Z) = Zm p + Nm n — B(A,Z) il Q di reazione e la differenza 
delle due binding energy : 

Q = 2B{~)-B{A,Z) (8.4) 
Sfruttando ora il modello a goccia si ha che la binding energy puo essere scritta com^] 

B(A, Z) = ay A - a s A 2 / 3 - a c -^ - a sym (A (8.5) 



Sostituendo la (8.5) nella (8.4) si cancellano tutti i termini tranne quello di superficie e quello 
Coulombiano. 

Q = a s (l - 2 1 / 3 )A 2 / 3 + a c (l - 2" 2/3 )J^ = -0, 26 • a s ■ A 2 / 3 + 0, 37 • a c ■ (8.6) 

Si noti che il termine di superficie inibisce la fissione mentre il termine Coulombiano la fa- 
vorisce : in effetti aumentando la superficie si creano elettroni meno legati ma si guadagna in 
repulsione Coulombiana. Per ricavare un limite utile per la fissione occore quindi porre Q > 

f > ^^ 0,7027.^18,2 (8.7) 
A 0, 37 ac ac 



avendo sostituito i valori di as = 18, 56 MeV e ac = 0, 717 MeV. Questa rappresenta una 

88c 
38 k 



condizione necessaria per la reazione ma non sufficiente. Ad esempio si prenda il nucleo IgSn : 



per questo Z 2 /A ~ 16,41 quindi la fisisone non e possibile. Vi sono poi alcuni elementi come 
^ 6 Pd con Z 2 /A ~ 19,96 e |] 8 Pb con Z 2 /A ~ 32,32 che rimangono stabili. II nucleo ^ 8 U ha 
Z 2 1 A ~ 35, 56 e stabile perche decade molto lentamente con una vita media di r ~ 10 16 yr . II 
Q di questa reazione e di circa 177 MeV. Gli ultimi casi sono rappresentati da nuclei artificiali 



1 I1 termine di pairing verra ora trascurato anche se dopo si osservera che questo e di fondamentale importanza 
nelle reazioni di fissione : da questo dipende , ad esempio, la fissionabilita dell' 235 U piuttosto che dell' 238 U. 
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ij^Fm con Z 2 /A ~ 39,06 e r ~ 2,63 hr e infine f^No con Z 2 /A ~ 40,32 e r ~ 1,2 ms. I 
tempi di vita media per i vari nuclei sono riportati graficamente nell'immagine seguente. 




(8.8) 

La deviazione dall'andamento lineare e dovuta ad effetti di shellj^] Dal confronto con i tempi 
di vita media si vede che la condizione Q > non e sufficiente a garantire la fissionabilita 
degli elementi. Infatti si deve pensare che i due nuclei figli siano formati per fluttuazioni nel 
nucleo iniziale e che, quindi, debbano oltrepassare una barriera di potenziale del tipo 



U(r) 



e 2 Z 2 



U(Ri + R 2 ) ~ 250 MeV 



(8.9) 



dove R\ ed R2 rappresentano i raggi nucleari dei nuclei figli {R\ + R2 = 2R nel caso simmetrico 
)• 



2 Swietecki propose una legge sperimentale che tenesse conto di questa correzione nella massa degli elementi 
M(A, Z) = Zm v + Nm n — Bldm — Bsheii (LDM sta per liquid drop model) . Per la vita media si estrapola 



una legge del tipo 



log t = log T + kSM 



dove k ~ 5 e M e la massa dell'elemento. 
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Separation distance 



La barriera puo essere attraversata per effetto tunnel ( effetto con probability rilevante solo 
nel caso di nuclei fissili spontaneamente ) oppure fornendo energia di eccitazione. 




-■?;B'a' " 



8.2 Deformazioni del nucleo 

Per una trattazione piu accurata della fissione occorre considerare le deformazioni che il nucleo 
subisce durante il processo. In effetti questo subisce deformazioni della forma , originalmente 
sferica, come mostrato in figura. 

Si indichi con E l'energia totale del nucleo. Ovviamente l'energia dipendera dai parametri di 
deformazione : E{e) = M{e) (c = 1) con e = (ei,£2, • • • , £n)- Se il nucleo e sferico quindi 




8.2. DEFORMAZIONI DEL NUCLEO 
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Si = , Vz , altrimenti viene deformato. II raggio del nucleo si puo scrivere , per la (7.1) 
come 

oo +A 

1 + E E a V*A, M (^) )> (8.10) 

A=l n=-\ 



R(0, <p) = A a R { 



Utilizzando questa relazione si puo esprimere il volume in termini dei parametri di defor- 
mazione 

V(e)= dip sintfdW r 2 dr = -ttR^ (8.11) 
Jo Jo Jo 3 

Nel caso di simmetria assiale 

R(#,(p) = A a R |l + f^anPn (cos 0) | 

L'energia totale del nucleo si puo quindi scrivere come 

E(e) = M(e) = Zm p + Nm n - B(e) (8.12) 
E(0) = M(0) = Zm p + Nm n - B(0) : nucleo sferico (8.13) 

Definendo U(e) = E(e) — E(0) = —B{e) + -B(O) si nota che il nucleo risulta stabile quando 
B(0) > B(e) , ovvero per U(e) > 0. Questo rappresenta il primo criterio per decidere se un 
nucleo puo fissionare spontaneamente : se U(e) e negativa il nucleo tendera piu facilmente a 
fissionare spontaneamente. Per scrivere esplicitamente i B{e) utilizziamo il modello a goccia 
liquida 

B(e) = a v V - E s (e) - E c {e) - E sym (e) (8.14) 
B(0) = a v V - E s {0) - E c (0) - E sym {0) (8.15) 

Le ipotesi sono che il volume si mantenga costante, che la carica elettrica sia uniformemente 
distribuita e che non ci siano addensamenti di neutroni o protoni (Z/N ~ costante) : in questo 
modo il termine di volume e di simmetria si cancellano quando si calcola la differenza tra B(0) 
e B{e). Restano da calcolare i termini 

E s {e) = (75(e) (8.16) 
E c{e) = ~^p\ [ d 3 r / d 3 r'^-^ (8.17) 



v Jv 



e analogamente per il caso sferico 



Eg(0) = aAnRl = 4iTarlA 2 / 3 = a s A 2 / 3 (8.18) 

Si ha che U(e) = Es(e) — Es(0) + E(j(e) — Ec(0). Si supponga quindi , per semplicita , che il 
nucleo venga deformato in un ellissoide di rotazione con un piccolo parametro di deformazione 
e. In questo caso il raggio si scrive come R($) = i?o(l + a2-P2(cos$)). Gli assi dell'ellisse sono 
dati da 

(a = Rn(l + e) 
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Si noti che , con questi dati, e soddisfatta la condizione V = ^irab 2 = ^ttRq. Si pud mostare che 
la superficie aumenta proporzionalmente alia deformazione della sfera , quindi S = 47ri?Q(l + 
|e 2 + . . .) = S(e). II termine Coulombiano invece , durante la deformazione , decresce come 
e 2 /5. AlPordine phi basso in s ( escluso lo zeresimo ovviamente ) si trovano quindi le condizioni 

E s (s) ~ E S (0)(1 + ^e 2 ) (8.21) 

E c (e) ~ E C (0)(1 - ^e 2 ) (8.22) 

Dunque si ricava 

U{e) = 2 -e 2 E s (0) - \e 2 E c {V) = \(2a s A 2 ^ - a c ^)e 2 (8.23) 

Se il primo termine e piu grande del secondo la differenza di energia e positiva : si acquista 
energia nella deformazione del nucleo e quindi maggiore e la deformazione, maggiore ener- 
gia viene acquistata. Questo tipo di nucleo e instabile sotto deformazione e tende quindi a 
fissionare SPONTANEAMENTE. La condizione per la fissione si puo quindi scrivere come 

\acZ 2 A- 1 ' ? > > la s A 2 / 3 > 1 (8.24) 

5 5 2as A 

Questo ci suggerisce di introdurre un parametro di nssionabilita dato da 

l_Ecm_ac L Z 2 _ 
X -2E s (0)-2a s A (8 " 25) 

II valore x = 1 definisce un regime intermedio (valore critico) per la nssionabilita 

La goccia rimane stabile , rispetto alia deformazione, quando x < 1 ; viceversa la goccia e 
instabile. 

Elemento x 



ig 6 Pd 0,385 

^ 8 Pb 0,624 

IfXJ 0,686 

f^Fm 0,754 

2 f 2 No 0,778 



Si noti infine che il solo parametro e non e sufficiente a descrivere la forma del nucleo quando 
questo si fissa : si dovrebbe infatti descrivere con ulteriori parametri la tipica forma strozzata 
assunta dal nucleo durante il processo. Se si aumentano i parametri U diventa quindi una 
ipersuperficie : i punti di sella ( in 2+1 dimensioni ) rappresentano i punti di instability 
scelti dal nucleo , quello ad energia piu bassa sara il favorite Deformando il nucleo si ha 
una ridistribuzione dei livelli del modello a shell : questo causa un diverso andamento della 
U(s) con conseguenti nuove configurazioni (meta)stabili. II potenziale , dunque, si modifica e 
assume il tipico andamento a doppia buca ben visibile in flgura. 



8. 3. DISTRIB UZIONE DEI FRAMMENTI DI FISSIONE 
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\ 



DnltMior ' \ 



Nella seconda buca di potenziale si trovano gli isomeri di forma : elementi con vita media 
molto piu breve (infatti la barriera da oltrepassare e molto piu sottile) . 




MX' 



8.3 Distribuzione dei frammenti di fissione 

In figura e graficata la distribuzione in massa dei frammenti ottenuti dalla fissione termica 
dell' 235 U. 




70 K 9C ioc no :k iso 140 i*: ISO 

Mass numoe- 
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Si nota die la distribuzione e simmetrica , come ci si aspettava ma che il caso di fissione 
simmetrica ha poca probability di avvenire. Questo si puo spiegare notando che , in questo 
caso, la barriera Coulombiana ha un massimo visto che l'interazione EM presenta un massimo 
funzionale in Z\ = Z2 : in questi casi la fissionabilita spontanea e svantaggiata. Anlizziamo 
ora un altro grafico. 




'0 80 90 100 110 120 130 1*0 150 160 
Mass HimDer 



Se si confrontano le distribuzioni di frammenti per alcuni isotopi come 
si nota che la curva e shiftata : in realta la traslazione per elementi piu pesanti non avviene 
in modo perfettamente proporzionale. Infatti la distribuzione dei frammenti leggeri si sposta 
leggermente in avanti mentre quella dei frammenti pesanti rimane praticamente identica : 
questo e dovuto ad effetti di shell, infatti nei frammenti pesanti c'e una sovrapposizione di 
numeri magici , quindi si ha un elemento prodotto doppiamente magico, ovvero molto legato e 
stabile che tende a formarsi piu facilmente : il picco dei frammenti pesanti presenta quindi un 
picco su questo elemento e tende a spostarsene il meno possibile.Gli effetti correttivi dovuti al 
modello a shell sono ben visibili, oltre che nell'immagine precedente, anche in quella seguente 
dove viene rappresentata la variazione dell'energia di attivazione per la fissione in funzione di 
A. Si noti che la previsione teorica del modello a shell (linea piu sottile) presenta dei picchi 
proprio sui nuclei che hanno numeri magici. 



8.4. NEUTRONI DELL A FISSIONE 
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8.4 Neutroni della fissione 

I neutroni emessi durante la fissione hanno uno spettro di energia del tipo 

fn(E) = k (^^) 1/2 e-°' 775 ' £ (8.26) 

Questa e simile ad una Maxwell-Boltzmann ed e un tipico spettro evaporativo: il sistema ha 
phi neutroni di quanti ne possano tenere i frammenti quindi questi evaporano via subito dopo 
la fissione. Questo andamento e visibile in figura 




Heutron Energy [MeV] 



In questa seconda figura e invece rappresentata la probability di emissione (P u ) di un certo 
numero di neutroni (y) in relazione al valor medio (V) , per varie fissioni. 
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Figure 13.7 Distribution ol lission neutrons Even though Ihe average number of neulrons t. changes with the 
lissionmg nucleus. Ihe distribution aboul the aveiage is independent ot Ihe original nucleus From J Terrell in Physics 
and Chemistry ot Fission. Vol 2 (Vienna IAEA ig6t>). p. 3 



Si nota che , nonostante il numero medio di neutroni cambi a seconda del nucleo fissionato, 
la distribuzione attorno al valor medio rimane pressoche identica e indipendente dal nucleo 
originale. Oltre ai neutroni pronti , che costituiscono lo spettro evaporativo appena descrit- 
to, vi sono anche dei neutroni ritardati : questi sono neutroni emessi in certi casi dopo un 
decadimento (3 del frammento della fissione. Se si immagina di proiettare il processo nel piano 
N — Z si puo descrivere l'emissione di neutroni ritardati nel modo seguente : dalla fissione 
si generano due frammenti piu neutroni pronti ; in seguito i frammenti decadranno f3 per 
raggiungere l'isobaro (A =cost.) stabile. In generale pero dopo il decadimento j3 accade che il 
nucleo risultante si trovi in uno stato eccitato che non ha ancora raggiunto la valle di stabilitd 
(sempre nel piano N — Z) : in questo caso puo essere emesso un neutrone dopo un ritardo~sec. 
I neutroni coinvolti nella fissione nucleare vengono divisi in : 

Neutroni Termici : E ~ 0, 025 0, 1 eV 

Neutroni EpiTermici : E ~ 1 4- 10 eV 

Neutroni Lenti : E ~ 1 4- 10 KeV 

Neutroni Veloci : E ~ 100KeV4-10 MeV 



8.5 Sezione d'urto per fissione indotta 

Per la scelta di un materiale fissile ad alta efficienza conviene considerare varie caratteris- 
tiche del materiale come Abbondanza relativa, tempo di dimezzamento e sezione d'urto. Per 
quanto riguarda le prime due caratteristiche si nota, dalla tabella seguente, che tra gli isotopi 
dell'uranio sono disponibili Puranio-235 e il 238. 
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A Stabile Abbondanza (%) 



233 
234 
235 
236 
237 
238 
239 



•1/2 



Decadimento 



/ 

/ 



0,770 
99,275 



0,1592 Myr 
0,245 Myr 

23,45 Myr 
6,75 day 

23,5 min 



Nel grafico seguente viene riportata la sezione d'urto dei due isotopi \J , dS U. 
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Questa curva contiene moltissime informazioni. Si nota inanzitutto che il comportamento della 
sezione d'urto per energie comprese tra 1'eV e il migliaio di eV e assimilabile alle risonanze 
dovute all'assorbimento del neutrone. Inoltre si osserva che utilizzando come materiale fissile 
235 U conviene utilizzare neutroni di bassa energia : la sezione d'urto per neutroni termici e 
circa 600 volte la sezione d'urto per neutroni veloci. Nel caso si vogliano utilizzare i neutroni 
ad energie del MeV prodotti dalla fissione per indurre nuove fissioni questi dovranno essere 
rallentati ad energie di neutroni termici in modo da aumentare la sezione d'urto. Purtroppo 
l'uranio-235 e poco abbondante in natura quindi bisogna utilizzare l'uranio-238 utilizzando 
neutroni ad energie molto alte. Un moderno reattore consiste quindi di due parti 



Materiale Fissile : Da origine alia fissione : Uranio. 



Moderatore : Rallenta i neutroni per portarli all'energia termica. Per moderare , trascu- 
rando processi anelastici tra neutroni e moderatore , servono moderatori con A tale 
che 

• E > 0, 1 MeV : A piccolo 

• E > 1 MeV : A grande 

Queste due parti possono essere accorpate in un unico apparecchio : in tal caso si parla di 
reattore omogeneo, altrimenti eterogeneo. 
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8.6 Energia della fissione 

Si consideri il processo 



n+ 235 U 



236 tt* 



U* — > fissione 



(8.27) 



La differenza tra i livelli energetici dei vari isotopi e visibile in figura 

236 U - = 23S U + n 



'Ji 



:3i u 

Si indichi con E cx la differenza tra la massa dei due isotopi (eccitato e no) delPuranio : 

E CX = M{ 236 U*)-M( 236 U) (c=l) 

Trascurando l'energia cinetica del neutrone ( neutroni termici ) l'energia dello stato eccitato si 
puo scrivere come M( 236 U*) ~ M( 235 U) + m n . L'energia E ex puo essere vista come somma 
di AE , il gap di energia tra i livelli di 236 U e 236 U* , e 5 , ovvero l'extra binding energy 
dovuta al termine aggiuntivo di pairing. Si ha quindi che 

E ex ~ M ( 235 U) + m n -M{ 236 U) = S n = 6, 54 MeV (8.28) 

dove S n e definita come separation energy del neutrone. Visto che _E a tt( 236 U) = 6,2 MeV< 
E ex allora 236 U puo essere fissionato senza energia spesa per i neutroni. Per la fissione dell' 
238 U invece si ha E cx = 4, 8 MeV mentre E att = 6,6 MeV. 

-39,j L£s , 




Quindi per fissionare quest'ultimo elemento servono neutroni di energia~ MeV. Si vede quindi 
che il termine di pairing 5 gioca il ruolo di inibitore o catalizzatore della fissione : l'energia 
di legame dell' 236 U e incrementata ( ovvero, l'energia dello stato fondamentale e diminuita ) 
di un valore 5 (~0,56 MeV) ; l'energia di eccitazione e quindi incrementata di 5. Nel caso 
delh 238 U l'energia dello stato fondamentale , prima della cattura del neutrone, e diminuita di 5 
e quindi l'energia di eccitazione e ridota anch'essa di 5. La differenza in energia di eccitazione 
tra questi due isotopi e quindi di circa 2<5 ~ 1, lMeV. 



8. 7. RALLENTAMENTO DEI NEUTRONI 
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8.7 Rallentamento dei neutroni 

Si vuole ora studiare la dinamica dell'urto tra i neutroni ed i nuclei del materiale moderatore. 
A tal scopo si indica con m la massa del neutrone e con v\{t) la sua velocita ; analogamente 
M e la massa del nucleo su cui il neutrone impatta e V2(i) la sua velocita (#2(0) = 0). 

Sistema del Laboratorio : Dopo l'urto il neutrone avra una velocita e la sua direzione 
sara deviata dai un angolo i?l rispetto alia sua direzione iniziale ; analogamente per il 
nucleo ti 2 , 0l- Trascurando la Binding Energy M ~ Am quindi 

TTl 1 

Vcm = -icjvi = ——rvi (8.29) 

m + M 1 + ^4 




Sistema del Centro di Massa : In questo sistema le due velocita sono u\ e ui prima del- 
l'urto, mentre dopo questo si hanno le velocita e it 2 con angoli rispettivamente di i? cm 
e /3 cm rispetto alia direzione iniziale. 




Con semplici conti di cinematica si possono esplicitare le relazioni tra i due sistemi: 

til = til- ti cm = j^tii = -^tii 

U2 = V 2 - V cm = -jtf&Vl = -ATT^l 

Inoltre 




(8.31) 



per l'elasticita dell'urto. L'energia persa dai neutroni durante l'urto e data dalla differnza 
tra l'energia iniziale Eq = ^mvf e l'energia finale E = ^m(u' 1 ) 2 . Con una semplice analisi 
vettoriale (vd.figura seguente) 
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si not a che 



1 — ^cm dunque 

(v'x) 2 = (u[ + V cm ) 2 



+ u'i + 2v cra u\ COS 



Sostituendo la (8.29) e la (8.30) nella (8.32) si ottiene , per il rapporto E/En 

1 + A 2 + 2 A COS tfcm 



E 



v? 



m 2 + M 2 + 2mM cos -& c 
(M + m) 2 



(A + l) 



Conviene introdurre una nuova variabile 



A - I 



< ft < 1 



A + l, 

in modo da riscrivere il rapporto E/Eq in termini di questa. 

E 1 



-En 



[(1 + a) + (1 - ft)costf c 



.32) 



3.33) 



(8.34) 



3.35) 



Da qui si vede subito che a < J| < 1 quindi Penergia varia tra -Emm = olEq < E < Eq = -E max 
, ovvero a = E^^/En. Si puo ricavare infine la relazione che lega l'angolo di scattering nei 
due sistemi di riferimento. Osservando che v[ cost?/, = v cm + u[ cos^cm si ricava , utilizzando 



la (8.29) , 



COS l?L 



+ u[ cos i? cm A cos $ cm + 1 

COS "i? cm + 1 



A+l 




E 



A-l 




(8.36) 



8.8 Distribuzione di probability per l'energia dei neutroni scat- 
terati 

Si supponga di essere nel caso di scattering a basse energie, per cui I = e la distribuzione 
angolare e isotropa nel riferimento del centro di massa. Se si indica con dW la distribuzione 
di probability che il neutrone venga trovato , dopo l'urto , con un angolo i? cm e data da 



dW = p($ cm )<M c 



Ldn 2^sintf cm dtf cm 1 



sin i5 cni di? 

C 



-p(E) dE 



3.37) 



47T 47T 2 

Infatti si deve calcolare il rapporto tra l'area sottesa dall'angolo solido individuato da i? cm 
e tutto l'angolo. Nell'ultima uguaglianza il segno e stato scelto negativo perche al crescere 
dell'angolo di scattering l'energia dei neutroni diminuisce. Si ottiene quindi , utilizzando la 



(8.35) 



dE 

dtfr.m 



En 



[(1 — ft) sin$ c 



dA: 



Quindi 



dE E (a - l)sini? c 
2 



= _ p(^,)M cm = 1 

■ '■ 2 En{a — 1) smv 



dE 2 E n (a- Dsintf™ (1 - a)E 

Quindi la probability e costante e non nulla solo per energie comprese tra a En e En. 



3.38) 



(8.39) 



8.9. MODERATORE IDE ALE 
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1. Teniamo fisso Eq e variamo A . In questo caso si vede che al variare di A la curva si sposta 
: ad esempio per l'idrogeno cxEq = 0. Se si graficano le varie situazioni si ottengono una 
serie di rettangoli che diminuiscono progressivamente la base e aumentano in altezza al 
crescere di A. 

2. Teniamo fisso A e variamo Eq : una rappresentazione di questo tipo e detta degli urti 
successivi. In questo caso Eq > > E^ > .... Facendo un grafico anche in questo 
caso si nota che la distribuzione tende a spostarsi verso lo 0. Piu realisticamente , al posto 
di una distribuzione di probability costante ( rettangolo) possiamo inserire una curva 
continua e notare che il picco della curva si sposta verso lo ad ogni urto successive 

Da quanto visto sembrerebbe conveniente definire una quantita , detta decremento logaritmico 
medio dell'energia dei neutroni , data da 

e=(log§>= r%(S)log§di?=l + ^loga = l+ ( 4^ 1)2 log^ (8.40) 
E J aEo E 1 - a 2A A + l 

Si noti che per l'idrogeno A = 1 £ = 1 , mentre invece per A—} oo, £ — > 0. Inline definendcj^] 
(log E) = log Ei si ha che 

£ = log E - (log E) = log E Q - log E x =>■ log Ei = log E - £ => Ei = E e" ? (8.41) 

II valor medio dell'energia del neutrone sara Ei = Eq dopo un urto, E2 = Ei = Eq e~ 2 ^ 
dopo due urti.. etc.. Questo ragionamento si pud iterare, ottenendo, per n urti , 

E n = E e"< n = {log |^ ^ (8.42) 

In particolare se si considera come energia finale quella dei neutroni termici e come iniziale 
Eq ~ 1 -T- 2 MeV la relazione appena ottenuta fornisce una stima di urti necessari per termaliz- 
zare i neutroni. Si puo quindi costruire una tabella dove si evidenzia l'efficienza del materiale 
moderatore. 



Elemento 


A 


a 




n 


H 


1 





1 


18 


D= 2 H 


2 


0,111 


0,725 


28 


12 C 


12 


0,716 


0,158 


114 


238jj 


238 


0,983 


0,00838 


2172 



8.9 Moderatore ideale 

II moderatore ideale dovrebbe essere poco costoso , facile da reperire, stabile dal punto di 
vista chimico, molto efficiente nell'assorbire i neutroni , compatto e con sezione d'urto di 
assorbimento molto alta. II moderatore piu comune e costruito con la grafite , altrimenti si 
utilizza l'acqua. Lo svantaggio dell'acqua e che la sezione d'urto della reazione n + p — > D + 7 
e tale che vengono persi molti neutroni. Per ovviare a questo problema si puo utilizzare acqua 
pesante o uranio arricchito. Lo svantaggio dell'acqua pesante e che risulta molto costosa. II 



3 L'esistenza di Ei e assicurata dal teorema della media. 
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prototipo di un reattore che sfrutta questo tipo di funzionamento e dato dal reattore di tipo 
CAND. Si introduca un fattore koo definite dalla seguente 

^ N(i + 1) n. neutroni termici alia (i + l)-esima generazione ^ 
x n. neutroni termici alia generazione z-esima 

La situazione ideale sarebbe quella in cui k^ = 1 : infatti se k^ > 1 il nucleo diventa 
instabile ed il reattore diventa supercritico. Analogamente , per k^ = 1 il reattore e detto 
critico, altrimenti subcritico. Si dovranno implementare meccanismi che non facciano uscire 
il reattore dalla zona k^ < 1. Si noti che questo fattore si puo scrivere come il prodotto di 
quattro fattori diversi. Inizialmente , al fine di semplificare la trattazione, conviene considerare 
che il reattore abbia dimensioni infinite. Procedendo per step progressivi si ha 

1. N(i) neutroni termici provocano la fissione. II numero di neutroni finali e dato da N(i)rj 
neutroni veloci. II parametro r] si puo calcolare tenendo conto di tutti i processi che 
partecipano alia fissione ma che fanno perdere neutroni. 

n. Neutroni che provocano fissione 
rj oc 



x 



n. Neutroni assorbiti 

Sezione d'urto macroscopica di fissione 



Sezione d'urto Macroscopica di fissione + Sezione d'urto di fissione per assorbimento 

= v 

+ S a 

(8.44) 

dove v e il numero medio di neutroni pesanti. Si ha quindi che 

V = ^ v S _f v (8-45) 

^F + ?->a 

dove la sezione d'urto macroscopica T,p e definita come T,p = ppap ove pp e la densita 
del materiale fissile e op e la sezione d'urto usuale. Questo rapporto e calcolabile : ad 
esempio per l'uranio naturale ( 238 U al 98,28% , 235 U alio 0,72%) si ottiene^l ??( nat U) ~ 
1,34 . Arricchendo ulteriormente l'uranio ( aumentando la percentuale di ) si puo 
aumentare il valore di r] in modo da avere piu neutroni pronti. 

2. Sono disponibili N(i)r] neutroni veloci che possono dare luogo a delle fissioni. Per tenere 
conto di queste fissioni spontanee dei neutroni veloci si introduce un nuovo fattore di 
fissione veloce e : tenendo conto di questo fattore il numero di neutroni diventa quindi 
N{i)r]£. Per l'uranio naturale e = 1, 03. 

3. Se si grafica la sezione d'urto di assorbimento di questi neutroni si nota che , per energie 
comprese tra 10 e 10 2 eV si trovano delle risonanze : i neutroni di queste energie vengono 
quindi assorbiti e persi. 

4 I dati sperimentali sono 

cr F ( 235 U) ~ 584barn;cr F ( 238 U) ~ Obarn; cr a ( 235 U) ~ 97barn; o- a ( 238 U) ~ 2, 75barn 
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Per tener conto anche di questo effetto si introduce un nuovo fattore p. Se p e la 
probability di evitare la risonanza, dopo il passaggio da questa zona di energia i neutroni 
sono diventati N(i)r/ep. Si noti che la presenza di queste risonanze impone un limite alio 
spessore dei moderatori : dovranno essere progettati in modo che il neutrone si trovi 
fisicamente dentro a queste quando incontra le risonanze. 

4. Si definisce inoltre un fattore / per tener conto della frazione di neutroni termici che 
vengano assorbiti. 

II fattore k^ si scrive quindi come 

koo = rjepf (8.46) 
Si puo dimostrare che, nel caso di reattore finite, la formula corretta e data da 

fc = Ml-fcO(l-*th) (8-47) 
dove si e tenuto conto degli effetti dei neutroni veloci (F) e termici. 



thermal neutron 



ll-fl 
thermal 
captures 
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8.10 Reattori Autofertilizzanti 

Si utilizzano materiali fertilizzanti che non sono fissili spontaneamente ma che lo diventano 
successivamente. Si puo utilizzare , ad esempio, Puranio-238. Si noti la catena di reazioni 

n + 238 U -»• 239 U -»• 239 Np + e~ + V con r ~ 23min (8.48) 
239 Np ->• 238 Pu + e~ + z7 con r ~ 2, 3days (8.49) 

Si noti che 238 Pu e materiale fissile ed ha (TF(Eth) = 742 barn. Un altro fertilizzante che si 
puo utilizzare e il gpTh. 

n + If Th -»• 233 Th -»• |f Pa + e~ + 17 con r ~ 22min (8.50) 
|f Pa -»• 233 U + e" + 17 con ~ 27days (8.51) 

Si noti che 233 U e materiale fissile ed ha o"^(-E t h) = 530 barn. Si noti inoltre che Puranio-233 
non e stabile , ma decade a con una vita media di circa 0,16 milioni di anni. Utilizzando 
questi metodi si costruiscono reattori del tipo FBRche utilizzano come materiale refrigerante 
il Sodio Liquido. 





II secondo disegno rappresenta una miglioria che e stata fatta ai reattori di prima generazione 
: in questi ultimi, infatti, l'acqua che veniva riscaldata nel nucleo era la stessa che azionava le 
turbine! 



8.12. SMALTIMENTO BELLE SCORIE 
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8.12 Smaltimento delle scorie 

Un grosso problema della produzione di energia tramite reazioni nucleari e , come noto, lo 
smaltimento delle scorie. In effetti i frammenti prodotti dalla fissione decadono con vite medie 
che possono arrivare fino a milioni di anni, come mostra il seguente grafico. 
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CAPITOLO 8. FISSIONE NUCLEARE 



Appendice A 

Particelle e loro proprieta 



A.l Quark 

A. 1.1 Classificazioni 



Sapore (Flavour) Massa [Mev/c 2 ] Spin (a) Num.Barionico Q/e S II T3 



u (up) 


1,7^3,3 


1/2 


1/3 


2/3 





+1 


+1/2 


d (down) 


4,1 -=-5, 8 


1/2 


1/3 


-1/3 





+1 


-1/2 


s (strange) 


~ 101 


1/2 


1/3 


-1/3 


-1 


+1 





c (charm) 


~ 1270 


1/2 


1/3 


+2/3 










t (top) 


~ 172000 


1/2 


1/3 


2/3 










b (bottom) 


~ 4190 


1/2 


1/3 


-1/3 











A. 1.2 Stuttura a quark degli adroni 
Nucleoni 



Pioni 



Kaoni 



Barioni S 



Barioni H 



I \p) = \uud) 
1 \n) = \uud) 

\tt + ) = \ud) 



\ir ) = \ud) = \tt + ) 



I \K°) = 


|ds) 




= us) 


'|E+> = 


| itits) 


< |S°) ^ 


uds) 




\dds) 




\dss) 


{!■•'." 


uss) 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



Barione 17 : \Q ) = \sss) 
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A.2 Mesoni 

A. 2.1 Reazioni di produzione e decadimento 



Particella Produzione Decadimento 



7T 


p + n — > p + p + tt~ 


7T - 


P + u , 


u ; 7T -)• e + ^ e 


n 
7T U 


P+p^-p+p + TT 




>7 + 7;7T°— )>7 + e + + e 


7T+ 


p+p— »p + n + 7r + 


TT+ - 


fi + + u fl ; 7T+ ->■ e + + z7 e 


K + 


it- + p ^ n + K + + K- K+ ->• + ^ ; 




f + 7T° ; ->• 7T+ + 7T+ + 7T- 


V 


7T + + p — )• p + 7T + + 7/ 


77— > 7r + + 7T° + 7r +7r + ;?7— >7 + 7;r/— > 3ir° 


P~ 


7r~ + p ->• p + p~ 




P -> 


7T+ -7T° 


n 

P 


7 + p ->■ p + p° 


P°- 




- . p o ^ ^0 + ^0 


1 

P + 


7T + + p — )• p + p + 






7T+ +7T° 


UJ 


7T + +p— ^p + 7T + +LiJ 


uj — > ir + 


+ 7T° + 7T 


+ + 7T _ ; a; — > tt° + 7 


A.2. 2 


Proprieta fisiche 










Particella Stranezza (5) 


Massa [Mev/c 2 ] J n 


T(T 3 ) 




7T+ 


135 


0" 


1 (+1) 




7T _ 


135 


0" 


1(-1) 




7T° 


139 


0- 


1 (0) 




1 


493 


o - 


1/2 (+1/2) 




K° 1 


497 


0" 


1/2 (-1/2) 




7? 


549 


0- 







p+ 


770 


1- 


1(1) 




p° 


770 


\- 


1(0) 




p- 


770 


1- 


1(-1) 




uj 


783 


1- 






A. 3 Barioni 

A. 3.1 Reazioni di produzione e decadimento 



Particella 


Produzione 




Decadimento 


A° 


7T+ + n -> A u + K+ 




> P + 7T 


; A u -> n + tt u 


A++ 


vr+ +p ->• A++ 




A++ - 


->■ p + 7T + 


A+ 


7T+ + n ->• A+ 




A+ - 


► n + 7r + 


A 


ir~ + p ->• A 




A — : 


' p + 7T~ 


A~ 


7r~ + n -)• A~ 




A- - 


> n + 7T~ 


E- 


it- +p ->• + xr 




E- — ! 


► n + 7r~ 


E+ 


if~ + p^ S + + 7T~ 


E+ - 


>P + 7T° 


; E+ ->• n + 7T+ 


e° 






E° -4 


•A + 7 




+ p -> S~ + A + ->■ 




~~ ->■ 


A° + 7T- 








H° — >■ 


A° + 7r° 


n~ 


a - - + p -> n - + k+ + k u n- - 


-> A u + K~ 


; fi" -> 


■ H u + 7T- ; 0" -»• S- + 7T U 



A.4. LEPTONI 
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A. 3. 2 Proprieta fisiche 



Particella 


Stranezza (S) 


Massa [Mev/c 


2 ] J n 


T(T 3 ) 




Protone 





938 


1/2+ 


l/2(+l/2) 




Neutrone 





939 


1/2+ 


l/2(-l/2) 




A++ 





1232 


3/2+ 


3/2(+3/2) 




A+ 





1232 


3/2+ 


3/2(1/2) 




A 





1232 


3/2+ 


3/2(-l/2) 




A- 





1232 


3/2+ 


3/2(-3/2) 




N* 





1440 -T- 1520 


1/2+ 


1/2 




A 


-1 


1116 


1/2+ 







£+ 


-1 


1189,4 


1/2+ 






s° 


-1 


1192,6 


1/2+ 


1(0) 






-1 


1177,5 


1/2+ 


1(-1) 






-2 


1321 


1/2+ 


l/2(-l/2) 




E° 


-2 


1315 


1/2+ 


l/2(+l/2) 




n~ 


-3 


1673 


3/2+ 







s*+ 


-1 


1382 


3/2+ 


1(+1) 




s*° 


-1 


1382 


3/2+ 


l(+0) 






-1 


1387 


3/2+ 


1(4-1) 






-2 


1532 


3/2+ 


l/2(+l/2) 






-2 


1535 


3/2+ 


l/2(-l/2) 




A. 4 Leptoni 












Particella Antiparticella Massa [MeV/c 2 ] Carica 


Spin r [s] 


Decadimento 


e e + 


0,511 


±1 


1/2 








<46eV/ 


c 2 


1/2 








105,66 


±1 


1/2 2,20-10~ 6 


-)• e + v e + Vfj, 


Vfx V fx 


j0,50 





1/2 






T - T + 


1784 


±1 


1/2 3,4 


• 10~ 13 -> 


fJ- + + u T ; ->■ e + v e + v T 


U T V T 


jl6,4 





1/2 
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Appendice B 

Relazioni Fondamentali 



B.l Meccanica Quant istica 

Proprieta degli operatori Riportiamo alcune proprieta fondamentali degli operatori uti- 
lizzati in meccanica quantistica. 

• Prodotto : fgt/j = f{gij)) 

• Commutatore : [/, g] = fg — <?/ / in generale. 

1. [f,[g,h}} + [g,[h,f]] + [h,[f,g}]=0 . 

2- [f,gh] = g[f,h] + [f,g]h. 

3- [fg,h] = f[g,h] + [f,h]g. 

4- [/,<?]* = -[/,<?]• 

5 - I X ' = -f( X ) > P er °g ni f( X )- 

6- if,9 k } = ^- 

7. In particolare [x,f(p)] = ih^ e \p,f(x)] = - ih^. 

8. Se [f,g] = c allora [e A ^, 5 ] = Ace A/ . 

• Se due operatori /, g commutano allora esiste una base di stati ortonormali e com- 
pleti {VVi} tali che fip n = f n ip n e g^ n = g n tyn- Ovvero esistono autostati simultanei 
di / e g. I due operatori sono quindi due quantita fisiche compatibili che possano 
assumere simultaneamente valori ben definiti. 

Operatori di p eq Gli operatori di posizione e di impulso sono definiti rispettivamente da 

• #(<?,*) = qtp(q,t) q^q 

• M?> t) = - ih^(q, 

Indeterminazione di Heisenberg Aq ■ Ap > h/2 , ovvero non si possono misurare con la 
stessa precisione sia la posizione che Pimpulso di una particella (pacchetto d'onda). II 
pacchetto d'onda Gaussiano minimizza questa relazione di indeterminazione. 

Evoluzione di un sistema isolato L'evoluzione di un sistema isolato si ricava dell' equazione 
di Schroedinger indipendente dal tempo 

ih-iP(q,t) = H(q,p,t)^(q,t) (B.l) 

H e definito come operatore Hamiltoniano : H(q,p,t). 



Ill 
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APPENDICE B. RELAZIONI FONDAMENTALI 



Evoluzione di un sistema isolato (non dipendente dal tempo) 

\ihd t tpn{q,t) = H(q,p,t)ip n (q,t) 



ihdtipn = E n il> n (B.2) 

Hlpn = Enlpn 



Quindi 



Equazione di Schroedinger 1 dimensionale 

h 2 d 2 



+ V(x)\^{x) = E^(x) (B.3) 



2m dx 2 

Questa equazione differenziale deve essere risolta 

1. Trovando i valori E n per cui esistono soluzioni (ricerca dello spettro) 

2. Trovando le ip n , autofunzioni relative agli E n ,con le condizioni 

(a) 11^(^)11 = 1 => stati legati , ovvero autovalori E n discreti. 

(b) Altrimenti se ||^(ar)|| > +oo basta richiedere che ip si mantenga limitata 
all'innnito (parte continua dello spettro), ovvero che appartenga all'insieme 
S = {4>(x) : lim| xKoo x N ■ ${x) = 0, ViV}. 

Buche di potenziale Si risolve l'equazione di Schroedinger in tutto lo spazio rispettando le 
condizioni al contorno. 

1. Buca infinitamente alta 

I , fuori dalla buca , , „ k„k 2 ir 2 h 2 9 . 

ip = I k = n-n a => E n = J} — = ~n 2 B.4) 

[sin(fcx) , dentro la buca 2m 2ma z 

4>n{x) = \ 2 sin —x (B.5) 
V a a 

2. Buca di profondita finita : Gli stati legati soddisfano il sistema 

£tan£ - n , n pan J2m\E\ , 
£cot£ = — rj , n dispari ,dove £ = — a e 77 = — a (B.6) 



2h 1 2h 



t2 , 2 _ ma z y 

Si risolve graficamente. 



Oscillatore Armonico L'Hamiltoniana e H = p 2 /2m + mu 2 x 2 /2 , quindi bisogna risolvere 

d 2 ib 2m . _ o 9 , , , 

_^ + _(£_ mw 2 x 2 /2)v , = () 

. Si ricavano gli autovalori 

E n = ^(2n + 1) = wft(n + 1/2) 
. La funzione d'onda dell'n-esimo stato e data da 
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Operatori di creazione e distruzione Si definiscono rispettivamente l'operatore di dis- 
truzione e creazione come 



(B.7) 



Invertendo gli operatori si ha 



(B.8) 



< ^ = yS(« + « t ) 

Applicando questo operatore alio staton -esimo si ha (ipn = ti) 

a\n) = y/n\n — 1) ; <r \n) = Vn + l\n + 1). 

Processo d'urto Processo in cui una particella libera attraversa una zona di potenziale che 
varia con la posizione. La funzione d'onda subisce quindi una modificazione nel passag- 
gio. L'onda incidente Y>£ = Ae lkx , A = m/h\k\ viene divisa in onda riflessa oc e _lfcr e 
onda trasmessa oc e lkx . 



Barriera di potenziale In questa condizione il potenziale e dato da 

V 



, per x < 0, x > a 
V > , per < x < a 



1. Se E > Vo le soluzioni sono 

V» 7 = e ikx + A e~ ikx con k = V^E/h (B.9) 

il>n = B e ik ' x + B' e~ ik ' x con k' = ^2m{E - V )/h (B.10) 

tl>iii = Ce lkx solo onda trasmessa (B.ll) 

Imponendo la continuita si ricava 

^ \_Jyras\_ _ _ _ _ ^ (B .l 2 ) 



tras | 

J~\ Ak 2 k 12 + (k 2 - k 12 ) 2 sm 2 (k'a) 

2 iJ2\2-2iJ, 



p = l^/z| _ (fc^-nW^q 

" \J inc \ 4k 2 k> 2 + (k 2 -k' 2 ) 2 sin 2 k> a " 



(a) D + R = leD^0: diversamente da quanto succede classicamente c'e la 
possibility che la particella attraversi la barriera (effetto tunnel). 

(b) Per E » V , D — >1,R — ► 0. 

(c) Per ^j2m(E - V )a/h = mr D = 1. 
2. Se E < Vo le soluzioni sono 

^/ = e fa + ie- fa (B.14) 
V> 77 = B e~ KX + 5' e KX con k = v^m^o - E)/h (B.15) 
tf; in = Ce ikx (B.16) 
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Quindi 



D 
R 



4fc 2 K 2 

4k 2 K 2 + (k 2 + k 2 ) 2 sinh 2 na 
(k 2 + k 2 ) 2 sinh 2 na 



Ak 2 K 2 + (k 2 + k 2 ) 2 sinh 2 kcl 

(a) D > sempre. 

(b) Se V -»• oo o a -»• oo allora L> ~ e - 2 V /2m (^-- B ) a / n . 

(c) Lo spostamento di coordinate provoca la comparsa di una fase. 
Gradino di potenziale II potenziale in questo caso si esprime come 



(B.17) 
(B.18) 



V 



se x < 
V > se x > 

1. Se -E < Vb le soluzioni sono date da 

tPj = e ikx + Ae~ ikx con k = V2mE/h 
$ u = Ce~ KX con k = ^2m(V - E)/H 

In questo caso R = 1 e D = 0. 

2. Se -E > Vb si hanno le soluzioni 

fa = e ite + Ae" fa con k = s/2mE/h 
= Ce ik ' x con k' = ^2m(E - V )/h 

Ponendo le condizioni al contorno si ricava 

Akk' 



D = 
R = 



(k + k') 2 
(k - k') 2 



(B.19) 



(B.20) 
(B.21) 



(B.22) 
(B.23) 

(B.24) 
(B.25) 



(k + k') 2 

Buca di potenziale 5 In questo caso V(x) = —g8(x) con g > 0. Le condizioni al 



contorno da porre sono 

fv(o„) = v(o+)^V(o) 
\v , (o + )-V'(0-) = -^(o) 

1. Sc E < (spettro discreto) 

^(x) = Jk (#(x) e kx + 0(x) e- kx ^j con k = 



(B.26) 



-2mE 



h 2 



2. Se E > (spettro continuo) 



^( x ) = t?(_ x ) 



-ik'x 



+ 0(rr) 



-iJfc': 



con fc' = v / 2mE/h 2 . Le ampiezze sono legate dalla matrice di transizione 



1 + ia 



ia 



la 



ia/ V£ 



(B.27) 
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Barriera di potenziale 5 E' simile al caso precedente, ma stavolta g < 0. Non ci sono 
stati legati e lo stato di diffusione generale e espresso dalla stessa ip con i coefficienti 
del caso precedente. 



D 
R 



1 

a 2 
1 + a 2 



(B.28) 
(B.29) 



Momento angolare Si defmisce il momento angolare classicamente ma utilizzando la for- 
mula quantistica per Pimpulso. 



L = f x p = — \h(f x V) 



(B.30) 



Separando le componenti si ha 



' L x = -\h{yd z - zd y ) 
L y = — \h{zd x — xd z ) 
L z = -ih(xd y - yd x ) 



(B.31) 



Ovvero 



Li = e ijk Xjp k dove e ijk = < 



1 se (ijk) = (123), (231), (312) 
-1 se (ijk) = (321), (132), (213) 
Oaltrimenti 



(B.32) 



Commutatori Valgono i seguenti commutatori 







'[Li,L 2 ] 


= itiL3 


Li, Lj] — \heij k L k — 


^> < 


[L 2 ,L 3 ] 


= \hL\ 






AL 3 ,L 1 ] 


= ihL 2 



(B.33) 



• Le componenti del momento angolare sono operatori Hermitiani. 

([L i ,x j ]=[he ijk x k ^ B34 ^ 
\[Li,pj] = ihe jik p k 

• Se si defmisce I? = l? x + L 2 + L 2 allora [L 2 , Li] = per i = 1, 2, 3 

• L 2 ed L3 possono avere simultanaemente misure definite 

Momento angolare orbitale In generale il momento angolare totale si indica con J ed e 
composto dallo spin , S e dal momento angolare orbitale, L. Per J valgono le stesse 
proprieta gia dimostrate per L. 

Operatori di salita e discesa per L Si possono definire gli operatori di salita e discesa 
seguenti. 

[ Li = L\ + iL 2 
r r • j (B.35) 
L_ = L\ — 1L2 
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Valgono le seguenti proprieta , di immediata verifica. 
'[L+,L_] = 2L 3 

< [L 3 , L + ] = L + ; L 2 = L + L_ + h\ - L 3 = L_L + + L 2 3 + L 3 (B.36) 
JL 3 ,L_] = -L_ 

Matrici di Pauli Per le particelle a spin 1/2 , j = 1/2 si possono definire le matrici di Pauli 
, che rappresentano l'insieme degli elementi di matrice di J«. 

;) v) -°o 

Le matrici di Pauli godono delle seguenti proprieta 

• [<Ti/2, crj/2] = ie ijk a k /2 

• *? = I 

• 0i( jj = -Oj(Ji = ie ijk a k (i / j) 
Spinori Vengono definiti due spinori di base 

Q = |l/2, 1/2) ee I t) ; (J) = |l/2, -1/2) ee I I) (B.38) 

attraverso i quali si possono scrivere tutti gli elementi dello spazio , indicati con 

Composizione dei momenti angolari Dato un sistema formato da due particelle con stati 
IV'iMV^) la funzione d'onda del sistema e data da (considerando stati puri) = 
IV'i)!^)- I momenti angolari si sommano 

J = Ji + J2 ; [J-i, Jj] = itijkJk 

Per studiare gli autovalori si possono scegliere due basi 

1. La base in cui sono diagonali J 2 , J 2 , Jiz, foz- 

• Si scelgono i vettori di base \ji,j 2 , mi, m 2 ) 

• Valgono le equazioni agli autovalori 



Jf\ji,j2,mi,m 2 ) = + l)h 2 \ji,j 2 ,mi,m 2 ) 
Jiz\ji,32,mi,m 2 ) = mih\ji,j 2 ,mi,m 2 ) 



(B.40) 



e analogamente per jf e J 2z . 
2. La base in cui sono diagonali J 2 ,J 2 ,J 2 , J z 
• Si scelgono i vettori di base \ ji, j 2 , J,M) 
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Valgono le equazioni agli autovalori 



J 2 \h,j 2 , J, M) = ji(h + i)^ 2 |ji,j 2 , J, M) 

jf |ji , j2 ,J,M)= 32 (j 2 + 1 ) h 2 1 ji , 32 , J, M) 
J 2 \juj 2 , J, M) = J(J + l)H 2 \j 1 ,j 2 , J, M) 
Jz \ji , 32 ,J,M) = Mh\h , 32 , J, M) 



(B.41) 



Segue 



Autovalori M di J 2 Poiche J 2 = J\ z + J2 Z => M = m\ + rri2 ■ Quindi all'autovalore 
mi sono associati (2ji + 1) valori , mentre ad 771,2 corrispondono (2j2 + 1) autovalori 
: in totale (2ji + l)(2j 2 + 1) 

Autovalori Jdi J 2 II valore massimo dei due momenti angolare e tale che J = 31+ 32 
'■ ji e 32 assumono il valore della proiezione del momento angolare ttt-i, 7772. Quindi 
\ji,32, mi-, m 2) < — > \ji,32, J = ji + J2, M = mi + m.2). H valore minimo possibile 
di J e |ji — J2I , ovvero J puo assumere i valori 

J = 3i+ 32,3i + 32 ~ 1, \ji ~ h\ 

Coefncienti di Clebsch-Gordan Sono necessari per passare tra le due basi del momento 
angolare. Per ottenerli basta sviluppare gli stati di una base in termine dell'altra base e 
viceversa. 

|ii,j2,"7i,777 2 ) = j 2 , J, M)(j!, j 2 , J, M\j!, j 2 ,m 1 ,m 2 } (B.42) 

J,M 

\ji,32,J,M)= ^ \ji,j2,m 1 ,m 2 )(ji,j2,m 1 ,m 2 \ji,j2,J,M) (B.43) 



mi ) ni2 



Vale ovviamente 



(ji,j2,J,M\j 1 ,j2,m 1 ,m 2 ) = (ji,j2,m l ,m 2 \ji,j 2 ,J,My 



(B.44) 



Spin Rappresenta il momento angolare intrinseco delle particelle ( non e legato alia massa 
delle particelle). Per le particelle che possiedono spin non nullo la descrizione dello stato 
attraverso la funzione d'onda deve poter determinare la probability che lo spin della 
particella abbia direzione definita nello spazio. 



ip(r) = ip(r, s) = 



( ^i(r) \ 



V/Wi(r)/ 



(B.45) 



Spin 1/2 Gli operatori si scrivono come 



S 2 = h 2 s(s + l) = h 2 3/4 
Si = (Jih/2 



S± = S x ± iSy =^ 



1 




;S. 





1 



(B.46) 
(B.47) 
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Valgono le proprieta 



a 



e ia ' CT = cos | a| + i r=T o : sin|a| (B.48) 



a 



n ■ s = -n- a (B.49) 



Problemi 3-dimensionali Si risolvono introducendo le nuove variabili 

J r = r i — r*2 

L mi+m 2 

L'Hamiltoniana si trasforma quindi in 



(B.50) 



H = -—^- rV| - ^V r 2 + V(r) (B.51) 

2(mi +m 2 ) K 2fi ' W V ; 

Dunque la funzione d'onda si divide in 

^ = $(f)0(f) 

Simmetria Centrale Nel caso di simmetria centrale V{r) = V(r). Si cercano quindi soluzioni 
della forma Y>(r) = R(r)Y^ m ($ , (p) , composte da parte radiale e parte angolare con 
armoniche sferiche. Per la parte radiale si deve risolvere l'equazione 



che si risolve facilmente introducendo x(r) = R(r)r . Si deve risolvere il sistema 
'& + { 2 w(E-V(r))-m^} X (r)=0 

/ °°dr| X | 2 = l (B-53) 

,x(o) = o 

Dunque ogni stato stazionario di un sistema 3D a simmetria centrale e identificato dalla 
terna di numeri quantici (n,£,m). 

Particella libera Si consideri il caso V = 0, Vr. Si ottiene l'equazione 

R" + 2 -R' + (k 2 - ^J^) R = ; k 2 = 2mE/h 2 (B.54) 



Se £ = allora sono soluzioni le ONDE PIANE 

J R = regolare in 



r 

[R = A'^r 11 singolare in 
Se £ / si pone Rg = r £ rj£ e si risolve 

Sono soluzioni le ONDE SFERICHE 

(R k: e = 2kj e (kr) 
[Qk,e = 2kn £ (kr) 

dove si sono definite le funzioni di bessel sferiche je(kr) e ne(kr). 



(B.55) 



(B.56) 
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Funzioni di Bessel-Hankel Si definiscono le funzioni di Bessel sferiche 



Per x — > , 



Per x — > oo 



'^W = (-iF(^)T 



a* (2£-l)!! 



. , N 1 (£ + 1)tt , 1 . (£ + 1)tt 
je{x) ~ - cos x ; Hi) ~ - sin x — 



(B.57) 



x 2 l/ x 2 

Sono utili anche le seguenti deflnizioni di funzioni di Hankel sferiche 

(hf\x) =je{x) + in e (x) 
\hf\x) =je{x) - in e (x) 

Gli andamenti asintotici si ottengono utilizzando le funzioni di Bessel. 



(B.58) 



Sviluppo in onde parziali Le soluzioni per momenti angolari diversi si possono correlare 
tra di loro : un'onda piana si puo sempre sviluppare in termine di onde sferiche utiliz- 
zando 

oo 

e ikz = e ifcrcos^ = {21 + 1) i^(fcr)P/(cos<? (B.59) 
e=o 

Buca 3D Prendiamo il potenziale a simmetria centrale 

I U altrimenti 

Si risolve l'equazione di Schroedinger , trovando le soluzioni 

• Per r < a 

Rf *) = Aje (kr) con k ee M ^ 2 + V ° ] > (B.61) 

• Per r > a 



R<f^ = Bh { p (i«r) con k' = W 2 ™ E = i« (B.62) 

Analogamente a quanto fatto per la buca unidimensionale si risolve l'equazione grafica- 
mente trovando che 

• Per \j2mV^a 2 j% < tt/2 non ci sono stati legati 



Per tt/2 < ^2mVoa 2 /h < 3ir/2 c'e un solo stato legato 
Per 2>tt/2 < ^/2mVoa 2 /h < 57r/2 ci sono 2 livelli discreti 
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Atomo di Idrogeno L'Hamiltoniana delPatomo di idrogeno si pud scrivere come 

ft 2 e 2 

H = V 2 con m ~ m e (B.63) 

L'equazione da risolvere e dunque 

R „ + lg_fil±ll R+ ^( E+ f) R = (B . 64) 
r r z n z \ r J 

• Gli stati legati sono possibili solo con E < . 

• La soluzione dell'equazione con il metodo delle serie di potenze porta alia condizione 
di quantizzazione dell'energia 

'^W^™^^^ 1 ' 2 ' 3 -> (B - 65) 

• L n-esimo livello e Efco^ 21 + 1 = n ^ volte de 

genere. 

• La soluzione per l'equazione radiale e data da 
ove 



2 _ 3 /2 (n-l-1)! 

^=-^ y {n+m * 

e -k 2 /^ 1 sono i polinomi associati di Laguerre. 
Per i primi stati 

• R 1:0 (r) = 2r B 3/2 e" r / r s 

Dal modello quantistico si ricava che il moto dell'elettrone nel nucleo e non relativistico 
(v ~ c/137) e che l'energia di ionizzazione e di circa !4eV 



Appendice C 

Equazioni di Lane Emden 

Viene presentato un articolo in cui le equazioni di Lane Emden vengono risolte analitica- 
mente. 
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23 Lecture 23: The Lane-Emden Equation 

"Science is facts; just as houses are made of stones, so is science made of facts; but a pile of stones 
is not a house and a collection of facts is not necessarily science." 

Henri Poincare 

The Big Picture: Today we discuss the Lane-Emden equation, which describes polytropes in 
hydrostatic equilibrium as simple models of a star. We also derive the Chandrasekhar limit for the 
formation of a black hole. 

The Lane-Emden Equation 

Last time we introduced the polytropes as a family of equations of state for gas in hydrostatic 
equilibrium. They are given by the equation of state in which the pressure is given as a power-law 
in density: 

P = np i , (462) 

where k and 7 are constants. The Lane-Emden equation combines the above equation of state for 
polytropes and the equation of hydrostatic equilibrium 



dP , GM(r) 



If we solve for the equation above for M(r) 

w , . r 2 dP dM Id fr 2 dP\ , t AS 

M U = -0lF => " "g^U*)' (464) 

and compare it to what we obtain from considering the spherical shell in hydrostatic equilibrium 

dM = 4irr 2 pdr =>■ = 4nr 2 p, (465) 

dr 

we obtain 



dM 


1 d 


(r 2 


dP 


dr 


~Gdr ' 


VP 


dr 




1 d 


(r 2 


dP 




r 2 dr 


VP 


~d7 



= -AirGp. (466) 

After inserting the polytropic equation of state [eq. (462)], the equation above becomes 

1 d ( r 2 _ 7 _idp 

After defining quantities 



^P"~ ~T = -**Gp. (467) 
r z dr \ p dr ' 



p = xe n , 

7 = — , (468) 
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the eq. (467) becomes 



1 d 

r 2 dr 



nr'n + l , l/wlU/n d{W 



X9 n n 



(X9 n y 



Tl + 1 l-n 

, „ ftA n 



dr 
1 d 



-4mG\0 r 



r2 d9 
r 2 dr V dr 



(469) 



We now make this equation dimensionless by introducing a radial variable £ 



a 



a 



II ■ 1 l-n 
ftA n , 



4ttG 



(470) 



to finally obtain the Lane-Emden equation for polytropes in hydrostatic equilibrium: 



1 



d 



a 



H) 2 d(aO 



d6 



1 d f 2 d9 



(471) 



This is a second order ordinary differential equation, which means that it requires two boundary 
conditions in order to be well-defined: 

1. Define the central density p c = A. Then 

p = \9 n => e{0) = 1. (472) 

— pg = —p c g = 0, because g c = (there is no mass inside zero radius). 



2. At r = 0, £ 



Therefore, 



dP _ 7 _!dp d0 
dr dr ^ d£ 



d£ 



0. 



(473) 



Analytic Solutions of the Lane-Emden Equation 

The Lane-Emden equation can be analytically solved only for a few special, integer values of 
the index n: 0, 1 and 5. For all other values of n, we must resort to numerical solutions. However, 
it is beneficial from both pedagogical and intuitive standpoint to derive these analytical solutions, 
which is what we do next. 

Analytic solution for n=0. 

After substituting n = into the Lane-Emden equation [eq. (471)], we obtain 




(474) 
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But, using the boundary conditions, we obtain 

(16 d6 1 1 - 

= =^> c 1= ^7 = -o^ => ^ = -7:^ + 02 

«4 5=0 «? 3 6 

=► 0(0) = 1 pa = l => O = l"k 2 - (475) 

6 

From the equation above, we see that this configuration has a boundary at £ = V6, where 9q — > 0. 
Analytic solution for n=l. 

After substituting n = 1 into the Lane-Emden equation [eq. (471)], we obtain 



Introduce the variable \ 



X(0=tf(0 (477) 



Then 



« «ve; ' (478) 

and the Lane-Emden equation in eq. (476) becomes 

=> %T = -f =► X" = "X => X" + X = 0. (479) 
This is a harmonic oscillator with general solutions 

X(0 = Asin£ + Bcos£, (480) 

or, in terms of # = x/£ 

, . . sin £ cos £ 
0(£) = A-^ + i?-^, (481) 

After imposing the first boundary condition, the general solution is obtained: 

cos £ 

9(0) = 1 => B = 0, because lim — - — = oo 

e^o £ 

sin£ 

A = 1, because lim — — = 1. 
=► = (482) 



The second boundary condition M 
rule 



?=q = is explicitly satisfied, because, after applying I/Hospital's 



£cos£-sin£ -£sin£ + cos£ - cos£ 1 

hm = hm —2 = — hm sm£ = 0, (483) 

S^o £ 2 s->o 2£ 2^0 s v 7 

as required. From the eq. (482) above, we see that this configuration is has a boundary at £ = tt, 
where 9\ — > 0. 
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Analytic solutions of the Lane-Emden equation 



Q. 




Figure 41: Analytic solutions for the Lane-Emden equation with n = 0, 1, 5. 
Analytic solution for n=5. 

The solution of Lane-Emden equation with n = 5 is analytically tractable, yet quite complicated 
to integrate. The solution is 

6 5 (0 = J_ . (484) 
This configuration is unbounded: £ G [0,oo), and lim^oo^s = 0. 

[For explicit derivation, see S. Chandrasekhar's An Introduction to the Study of Stellar Structure 
(University of Chicago Press, Chicago, 1939), p. 93-94] 



The Chandrasekhar Mass Limit 

Consider a star which has, through gravitational contraction, become so dense that it is sup- 
ported by a completely degenerate, extreme relativistic electron gas (i.e, p > 10 7 g cm -3 ). The 
pressure in terms of the density is obtained by combining the eq. (436) 



he ( 3 



1/3 



and 



P= -\-\ n 4 / 3 

. 7T , 



P 



(485) 
(486) 
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to obtain 



P 



he ( 3 



1/3 



4/3 



m H p e , 

(6.63 x 10" 27 erg s) (3 x 10 10 2p) 1/3 



8 



1 



'1.67 x 10 



-24 ^ 4 /3 



4/3 



P 



< 10 15 f — 

Me 



4/3 



(487) 



which is an equation of state for a polytrope with 7 = 4/3 and k 
of the index n = is n = 3. 



1.24xl0 15 

-4/3 



Corresponding value 



The mass corresponding to this polytropic configuration can be computed as follows: 



M 3 



p(r)d r = 4tt 



Xp(r)r dr = An 



\9 4 {aifd{ai) 



AirXa 



AttXc 



d 



e 



,d6 



di V d£ 



,d0 



(488) 



where we have used the Lane-Emden equation in eq. (471). The constant A is defined in eq. (470), 
and for n = 3 is 



a 



Tl + 1 l^n. 

re A n 



Acr 



A 



4ttG 

K ^2 

^r X 3 



a 



K -2 
7TU 



1 3/2 



K 13/2 



(489) 



The term in brackets can be evaluated numerically (Table 4.2 of Astrophysics I: Stars by Bowers 
&, Deeming) to about 2.02, so the total mass is 

1.24xl0 15 l 3 / 2 ,, /9 

1.24 xlO 15 r 2.02 



Itt 



-4/3 
Me 



tt(6.67 x 10" 8 
1.16 x 10 34 1.16 x 10 34 



2.02 = 4tt 

M 



tt(6.67 x 10-8) 



t'4 



M 3 



5.81 



1.99 x 10 33 



(490) 



Let us now compute p e for a star with relativistic matter degeneracy. In such a star, it is convenient 
to define the matter density, due essentially to the ions, as p = m#/i e n e . Also, let us consider 
contribution from hydrogen (subscript H), helium (He) and elements with atomic weight greater 
then 4 (Z). Then, from the definition in eq. (439), we have 



m H n H 2^m He n He 



A + -Y + -Z. 
2 2 



m, 



2 mz m Z n Z = PH + 2 PHe A p Z 



4 p 2A p 



(491) 
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Also, conservation of mass imposes that 

X + Y + Z = 1 => Z = 1-X -Y (492) 

so 

-I = X + 1 -Y + 1 -(l-X-Y)= 1 -X + 1 -= 1 -±^ 

Jx e 2 2 V ' 2 2 2 

- S = ^ - "-IT*- 

The stars that are undergoing extreme relativistic degeneracy of matter are highly evolved (near 
the end of their life-cycle) , which means that it is reasonable to assume that most of their hydrogen 
fuel has been burned up, so 

X w fi e sa 2. (494) 

Finally, we combine this result with the eq. (490) to obtain the Chandrasekhar mass limit: 

5 81 5 81 

M Ch = —M Q = - r M Q => M Ch = 1.45M . (495) 

When a star runs out of fuel, it will explode into a supernova or a helium flash (see Fig. 16). The 
Schwarzschild mass limit implies that star remnants with mass M > Mch cannot be supported by 
electron degeneracy and therefore will collapse further into a neutron star or a black hole. 
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